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Under salt water for an hour this Type K Potentiometer arrow) was returned 


to the factory for repair. Restored to its original accuracy, it is again in use. 


Water-Soaked Potentiometer 
Restored to Full Accuracy 


When the flood subsided after 
the New England hurricane ot 
1938, the men at the Marine 
Biological Laboratory, Woods 
Hole, Mass., began to take stock 
of the damage. ‘To their disap- 
pointment, they found that one 
of their ‘Type K Potentiometers, 
used to calibrate ammeters, volt- 
meters and wattmeters, had been 
under salt water for an hour. 
Corrosion soon began to appear, 
and the potentiometer seemed a 
total loss. 


After a thorough inspection, 
however, some resistors were re- 
placed, and salt deposits thor- 
oughly removed from the top 
plate and case. ‘The recondi- 
tioned instrument is again in use, 
at its original accuracy. 


EARLY TYPE K TURNS UP 


In addition to the calibrating 
potentiometer, there’s another in- 
teresting Type K at the Labora- 
tory ... serial number 10141 
purchased in 1908 by a power 


Where it has been 
since, and how it got to Woods 
Hole, is a lost story. 


company. 


But even 
after 32 years, this Type K is 
still in first-class mechanical con- 
dition. And electrically, the ini- 
tial stability of its resistors has 
It stands as 
a tribute to the skill and care of 
the instrument makers of another 


increased with age. 


veneration. 


Some of these men are still 
working at L&N with modern 
materials and, with this unique 
background of experience, today’s 
instruments can be expected to 
serve even longer. 

CONDENSED CATALOG 

Type K Potentiometers are de- 
scribed in Catalog E-50B(3), 
sent on request. 





LEEDS 


110 RESISTANCE STANDARDS 
GUARD L&N ACCURACY 


Stock Instruments Used in Making 
And Calibrating Other Equipment 


L&N’s working standards of 
resistance are checked against 
other L&N resistors used as pri- 
Working and 


primary standards are _ identical 


mary standards. 


to stock instruments except for 
these periodic records of resist- 
ance, 

There are 92 working stand- 
ards of resistance at L&N... 
18 primary standards. Working 
units are checked against primary 
units every 4 months. Primary 
units are certified at the National 
Bureau of Standards every year. 
Resistance of all is very constant. 


In practically every case the vari- 





Checking working standards in a constant 
(-+-o.1 C) temperature chamber at the L&N 
plant, Phila., Pa. Switches are turned and 
standards moved by rods projecting through 
the chamber. 

ation from year to year is well 
within the limit of error of the 
comparison. 

Used in calibrating every L&N 
instrument, these standards are 
typical of the constancy of L&N 
equipment in general. Many other 
L&N instruments are also used as 
primary and working standards. 
All are described in Catalog E 
which will be sent on request. 


LEEDS & NORTHRUP COMPANY, 4978 STENTON AVE., PHILA., FA. 


& NORTHRUP 


Measuring Instruments - Telemeters - Automatic Controls - Heat Treating Furnaces 
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recise Measurements of Frequency 
with 


G-R Measuring Equipment 











. 
[ ROM THE G-R PRIMARY STANDARD OF FREQUENCY, 
fundamental frequencies of 50, 11, 10, g, 1 and o.1 kilocycles 
are obtained. By means of harmonic frequency multipliers in the 
primary standard, a large number of harmonics of each of these fre- 
quencies are available, also, at the output terminals of the standard. 
For the precise measurement of frequencies in terms of these stand- 
ard frequencies, suitable auxiliary and interpolation equipment is 
needed. The G-R Frequency Measuring Equipment assembly, with 
| the Primary Standard of Frequency, makes possible the accurate, di 
| rect, precision measurement of any frequency from a few cycles up 
to 25 megacycles, or higher if external receivers are used. 
he units comprising the measuring equipment are as follows: 
Heterodyne Frequency Meter covering the range of 100 ke to 
5,000 ke in sixteen steps. Sufhcient harmonic output is obtained 
to produce beats in a high-frequency receiver operating at 30 Mc 
or higher. For many measurements, suthcient accuracy is ob- 
tained through use of direct-reading finder and interpolation 
dials on the heterodyne frequency meter. 
Interpolation Oscillator which is a direct-reading, linear-scalk 
audio-frequency oscillator covering frequencies between o and 
5,000 ke. It is used to measure the audio-frequency differences 
between the unknown frequency (or a submultiple thereof) 
and a standard 10o-kc harmonic. When the difference is very 
small, use is made of the 9-11 ke standard harmonics to avoid 
measurements with very low audio frequencies. 
Heterodyne Detector with plug-in coils covering the range 
from 25 ke to 25 Mec. This detector is used to obtain beats be- 
tween the standard and the unknown radio frequencies. 
Regenerative Selective Amplifier which is used to select any 
multiple of 1 ke between 1 and 10. This amplifier is particu- 
larly useful when the cathode-ray Comparison Oscilloscope is 
used in calibrations in the upper audio-frequency and lower 
radio-frequency ranges. 
Comparison Oscilloscope with 100-cycle and 1,000-cycle smooth- 
ing filters, networks for obtaining circular sweeps at these fre- 
quencies, and switches for connecting units of the frequency 
standard and measuring equipment to the oscilloscope. 
| The assembly is housed in a steel cabinet and includes a built-in 
loud-speaker. Permanent shielded connecting cables between the 
measuring equipment and the Primary Standard and for all inter- 
connections between units of the measuring assembly are furnished. 
All connections for actual measurements are made by suitable switch- 
ing on the centralized coupling panel, through which the standard 
and unknown sources and the various measuring instruments are 
connected. Suitable power supplies for a-c operation from 115-230 
volts and so—60 cycle mains are individually built in each instrument. 
| GENERAL RADIO COMPANY, Cambridge, Massachusetts. 
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ISTORIANS looking back a hundred years 

from now will be able to make a much 
more objective evaluation of the trying times 
which we are passing through than we are able 
to do now. Yet one can be fairly sure that, on 
the red side of the ledger sheet, will appear the 
large loss of life (unfortunately from the ranks 
of those most physically fit), the permanent 
injuries and diseases of those engaged in combat 
and those enslaved in concentration camps, the 
stifling of intellectual freedom, and finally the 
destruction of great masterpieces of art and 
architecture. The destruction of battleships, 
airplanes, and other .implements of war will 
probably not be an important item since these 
would become obsolescent in any case. Even 
the prodigious waste of labor in manufacturing 
these implements hardly constitutes a net loss 
because it is difficult to prove that busy men 
obtain less satisfaction out of life than those who 
have a great deal of leisure time. It may easily 
be that even the great destruction of buildings 
now going on will be placed on the profit side 
of the ledger because in 
the place where this de- 


and its results still stand in peace as well as in war. 

Professor J. C. Hunsaker brought this point 
out very forcefully in his paper at the Christmas 
meeting of the American Association for the 
Advancement of Science on progress in the aero- 
nautical industry. He discussed the possibility 
of applying much of the research now going on 
to the improvement of commercial airplanes. 
For instance, the use of the spectroscope in the 
study of the nature of combustion has led to 
greater power-over-weight ratios for engines that 
will aid planes used for non-military purposes as 
much as those in military service. The principles 
needed to design planes that can take off with 
heavy bombing loads may be applied with 
only small modification in peace times as well 
as in war times. Our increase in knowledge 
of aerodynamics will lead to improved wing 
sections with a greater “‘lift to resistance” ratio 
and thus lead to improved civil planes. Per- 
haps even such controversies as air cooling vs. 
liquid cooling will be settled by the research 
forced upon us by war conditions. 

The gist of this en- 
tire discussion can be 





struction now. stands 
new streets, modernand 
more efficient buildings 
will grow up. There is : 
: ; D. EVANS 
one item that is sure to 
re profit side; th: 
eon the pr rfit side ‘that Masinae. @ Chane 
isthe research which the 
intense competition of G. FAILLA 
war time has forced upon 
us. Research like time tt- 


self always goes forward 





For April 
Special Issue on Applied Nuclear Physics 


{pplications of Nuclear Physics, by ROBLEY 
Separation and Use of Stable Isotopes, by 
Biological Effects of Ionizing Radiations, by 


{lso abstracts of papers presented at the 
Conference on Applied Nuclear Physics 


summed up in a few 
words. The truly useful 
result of all our defense 
effort is the additional 
knowledge gained from 
the increased research 
program. Let us make 
sure that thisimportant 
fact is not forgotten by 
those setting up budgets 
for national defense. 
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Theory of the Plastic Properties of Solids. IT 


BY FREDERICK SEITZ 


Randal Morgan Laboratory of Physics, University of Pennsylvania, Philadelphia, Pennsylvania 


AND 
T. A. READ* 


Westinghouse Research Laboratories, East Pittsburgh, Pennsylvania 


Part A, Single Crystals 
2. THe THEORY oO! 


N the 


Stipe (CONTINUED) 


first installment of this 


outlined the basic facts concerning the plastic 


properties of single crystals and proceeded with 
a discussion of the theories of slip. It was found 
that the theory of dislocations provided the most 
plausible mechanism of the slip process suggested 
thus far, and it was pointed out that the experi- 
ments on internal friction in single crystals can 


be given a good qualitative interpretation in 





Fic. 19. Stress-strain curves for aluminum obtained in 
compression and tension. (After Taylor The experi 
mental values coincide for both types of deformation 


The ordinate is stress and the abscissa strain. 


terms of this mechanism. We shall now proceed 


to test this theory further by 


considering its 
application to other topics surrounding the slip 


| rocess. 


e. Work hardening 

The critical shearing stress for a given speci- 
men is strongly dependent upon the amount of 
previous slip that has taken place within it. 
This hardening influence of cold work is shown 
strikingly by conventional stress-strain diagrams 


1 The first installment appeared in the previous issue of 


this journal (February, 1941). 
* Westinghouse Research Fellow. 


170 


series, we 


in which the critical shearing stress is plotted as 
a function of the strain. Curves of this type are 
usually obtained in practice by placing a single 
crystal under tension or compression and ob- 
serving the stress required to produce plastic 
flow for a wide range of extension or compression, 
suitable correction being made for the change in 
dimensions of the specimen and orientation of 
the slip planes in relating the applied load to 
the actual shearing stress. Stress-strain curves 
of this tvpe are shown in Fig. 19 for aluminum? 
and in Fig. 20 for sodium chloride* and silvet 
chloride. 

features of this 


In order to appreciate the 


change let us consider the general behavior of a 


of MM2 











sS-—- 
Fic. 20. Stress-strain curves for sodium chloride upper 
curve) and silver chloride (lower curve). (After Stepanow 


Che ordinate is stress; the abscissa is strain. 


2(;. I. Taylor, Proc. Roy. Soc. 145, 388 (1934 
A. W. Stepanow, Physik. Zeits. Sowjetunion 8, 25 
1935). 
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typical specimen during an extension test. 
Starting with the virgin material, which we shall 
assume to be well annealed, we find that there 
is a region of stress in the immediate vicinity of 
the origin of its stress-strain curve in which the 
material behaves elastically. This range extends 
up to values of the order of 107 dynes/cm*, as 
we have seen previously, or over about one- 
fiftieth of the entire range of ordinate shown in 
Figs. 20 and 21. If the crystal is extended beyond 
its Hooke’s law range to a point such as A in 
Fig. 21, it does not return to its original length 
when the stress is relieved, but retains practically 
all of its extension. Generally speaking, there is 
a new elastic stress-strain region about the new 
length. This region usually extends over a larger 
range of stress than in the virgin material, as is 
shown in the figure, for plastic flow does not 
begin until the stress approximately reaches the 
value A. Similarly, if the material is further 
extended to B by additional slip, it will possess 
a new elastic range about the new length for a 
still wider region of stress. Thus cold work 


STRESS 





STRAIN 





Fic. 21. Behavior of a ervstal during straining. In the 
\ irgin state, its elastic range covers a region ot stress near 
the origin too small to be shown; after plastic flow to A, 
obtained by stressing, the crystal retains an amount of 
strain corresponding to O’, when the stress is relieved. 
Che range of stress in which Hooke’s law is valid is O’— A 
Similarly, after extension to B, the strain O” is retained 
and the new range of elasticity corresponds to O''—B. 
Lines such as O’—A and O”—B have a slight forward 
slope corresponding to Hooke’s law. 


increases the range in which the material is 
elastic. Within practical limits, the actual values 
of the elastic constants of pure materials do not 
seem to be altered as a result of extension. 
Taylor? has pointed out that in many simple 
metals the part of the stress-strain diagram 
nearest the origin can be fitted with a parabolic 
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Fic. 22. Stress-strain curves in zine at various tempera 
tures. Equation (1) does not appear to be valid in any of 
these cases. (After Boas and Schmid 


curve corresponding to the equation 
S*=ae, 


where S is the strain, o is the stress and a is a 
constant. The case of aluminum, shown in Fig. 
19 is an excellent example of this type. There 
appear to be exceptions to this rule, however, 
such as the case of zinc illustrated* in Fig. 22, 
in which the stress strain curves are practically 
linear over a range of temperature. Figure 23 
shows that the parabolic type of curve is valid 
in aluminum over a wide range of temperature,’ 
so that the difference between the curve for the 
two metals seems to be due to some underlying 
difference in the materials. 

It is interesting to note that more extensive 
straining can be produced by compression than 
by extension since the tensile strength for rupture 
is ultimately reached in a tension experiment. 
Thus compression strains of fifty have been 
attained by Taylor® in copper, whereas exten- 
sions of much more than ten are rare. Work 
hardening continues to the rupture point in 
tension tests; however, the critical shearing 
stress reaches a maximum value in compression 
tests after approximately tenfold compression 
has occurred and then remains constant. Thus 
about the same degree of work hardening is 
achieved for both kinds of strain in single 


crvstals. 


*W. Boas and E. Schmid, Zeits. f. Physik 61, 767 
(1930); 64, 845 (1930). 

5\W. Boas and E. Schmid, Zeits. f. Physik 71, 703 
1931). 

6G. I. Taylor and W. Farren, Proc. Roy. Soc. 111, 529 
1926); 116, 16, 39 (1927). 





An interpretation of work hardening on the 


basis of dislocation theory was first given by 


Taylor,? who pointed out that dislocations are 
surrounded by a stress field and hence should 


exert forces on one another. A straightforward 





investigation on the basis of elasticity theory 
. T ] 
| 
t + + + 
| | 














+ > TS ee 
= — 
ee 600° —— 
. 20 40 60 80 00 i20 40 160 80 
ENLONGATION IN % 
Fic. 23. Stress strain diagrams for aluminum at various 


temperatures. All of these are least 


the origin 


parabolic, at neal 


shows that at distances not too near the center 
of a long line dislocation the stress components 
fall off as the inverse first power of the distance 
from the line of dislocation. The angular varia- 
tion of the stress field’ is shown in Figs. 24a, b, c 
and d. The first of these illustrates the shearing 
stress in planes parallel to that in which the 
dislocation travels, the second shows the shearing 
stress in planes orthogonal to these and parallel 
to the axis of the dislocation line, whereas the 
third and fourth illustrate, respectively, the 

’ The mathematical equations for the stresses about a 


dislocation in an infinite isotropic medium are, in 
of the coordinate axes of Fig. 25, 


terms 


y 
a g = ( = ? 
a 
Xx 
g Oyz € 9 
Ty 


where ¢=G)X/z, in which d is the relative slip distance neat 
the center of the dislocation and G is the rigidity modulus 


of the solid 





tensile stresses in the same planes. The angular 
variable in all of these polar plots is the azimuthal 
angle in the plane normal to the line of the 
dislocation (Fig. 25), that is, the plane shown in 
Figs. 8, 9, 10 and 12. Figures 24a and 24b are 
identical because of the relation o,,=0¢,, between 
the shear components. Thus the shearing stresses 
are largest in the forward direction lying in the 
actual slip plane, decrease to zero in the direction 
normal to this plane, and reverse their sign in 
the 
stress in planes parallel to the slip plane is zero 


the next two quadrants. Similarly tensile 
in the slip plane, has a maximum in one direction ‘ 
normal to the slip plane and a minimum in the 
other, and so forth. Thus a dislocation lying in 
any direction from the given dislocation except 


pe 4 
| 2 
| GC; | 0, 
| *¥ y* 
| 
| | 
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(c) (d) 
Fic. 24. Polar plots of the angular variation of stresses 


about a dislocation. The polar angle lies in a plane normal 
to the line of the dislocation and is oriented in the mannet 
shown in Fig. 25; that is, the compressed part of the dis 
location is above the slip plane and the extended part is 
below. (a) is the shearing stress in planes parallel to the 
slip plane, whereas (b) is the shearing stress in planes 
orthogonal to these and parallel to the axis of the dis 
location line. (c) and (d) show, respectively, the tensil 
stresses in the same planes. 


along the 90° and 180° azimuths is subject to a 
shearing stress. 

For the present we shall confine the discussion 
to dislocations traveling in parallel planes and 
extending in the same directions. Simple con- 
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the 
repel one another and unlike 
ones attract. Hence unlike dislocations traveling 


siderations show that like dislocations, in 


sense of Fig. 9, 


in the same planes and having parallel lines 
would tend to get as close together as possible, 
that is to be on the same normal plane. For 
this reason, it is possible to form stable lattice- 
like arrays of positive and negative dislocations 
analogous to the lattices of positive and negative 
charges occurring in ionic crystals such as rock- 
salt (Fig. 26). Taylor investigated the properties 
of such lattices of dislocations under external 
shearing stresses and showed that even if the 
isolated dislocations would move under zero 
stress in an undeformed lattice, the lattice of 
dislocations will not shear unless the 


stress 
exceeds a critical value dependent upon the 
density of dislocations. 

It is easy to see from dimensional reasoning 
that this critical shearing stress will increase as 
the square root of the density of dislocations in 
a plane normal to their length. As we mentioned 
above, the shearing stress between two disloca- 
tions drops off as the inverse first power of the 
distance between them, whereas the density .V 
in the plane cutting across the dislocation lines 
drops off as the inverse second power of the 


average distance between neighboring disloca- 
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Fic. 25. The coordinate system chosen in Fig. 24. 


tions in the lattice. Thus the shearing stress ¢ 
required to the 


should vary as the square root of the density: 


overcome interaction forces 


o-=cN}, (1) 
where ¢ is a constant. Taylor has shown that ¢ 


VOLUME 12, MARCH, 1941 


has the value 
c= aGi, 


where G is the rigidity modulus of the solid, \ is 
the distance of slip produced by passage of one 


+ + + + + 

fe) fe O fe) O, 
¥ * + + salielae PLANES 
O fe) O fe) O : 

+ + + oo + 
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Fic. 26. A stable array ol dislocations. \fter Taylor.) 
lhe negative dislocations are indicated by circles, the slip 
planes by dashes. 


dislocation through a crystal, and a is a constant 
near unity depending on the precise arrangement 
of dislocations. 

With these facts in mind, Taylor suggested 
the following picture of the hardening process: 

1. Dislocations of opposite types are generated 
at block equal 
numbers, and move into the crystal blocks under 
the action of the stress. 


boundaries in approximately 


through the 
blocks, but become stuck, perhaps as a result of 


2. They do not pass entirely 
their mutual interaction, and form a more or 
less regular lattice structure. 

3. As the number of dislocations increases, the 
shearing stress required to make thetn move 
becomes higher and higher and the material 
hardens in accordance with Eq. (1). 

Now if the are 
uniformly distributed throughout all blocks, it is 
clear that the strain produced as a result of 


dislocations more or less 


their generation and motion will be proportional 
to their density. Thus, according to Taylor's 
picture, the critical shearing stress should vary 
as the square root of the strain S. Using more 
quantitative reasoning, Taylor has derived the 
following relation: 

a =kG(vAS/L)}, (2) 
where k is a constant of the order unity related 
to c in (1), and ZL is the linear dimension of a 
block. 

As we saw above, this type of relation is in 
excellent agreement with the type of stress-strain 
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curve found in certain metals, such as aluminum, 


chloride and_ silver 


chloride. In fact, by fitting Eq. (2) to the curve 


and in the salts sodium 
for aluminum, in the manner shown in Fig. 27, 
Taylor has found that the best value of ZL is 
about 5.3-10~' em, which, in order of magnitude, 
is the value we should expect from other evidence 
concerning the size of mosaic blocks. 

2) 


\ close analysis of the derivation of Eq. ( 
shows that an essential condition for this para- 
bolic type of relation is that the relative spacing 
of dislocations in the directions parallel and 
normal to the slip plane remain fixed during the 
slip process. It turns out that o would increase 
more rapidly with S if the spacing normal to the 
than that 
parallel to the slip plane, and, in fact, a linear 
connection 


slip plane decreased more rapidly 


between stress and strain could 


readily be obtained by assuming the proper 
relation between the two spacings. Thus cases 
such as zinc, in which the parabolic law is not 
valid (Fig. 22), do not necessarily imply that 
Taylor’s basic concept of the hardening process 
is incorrect. 

In any event, Taylor's theory gives rise to a 
number of questions on which it seems worth 
while to speculate. These are as follows: 

1. Is the length of a fully formed dislocation 
also of the dimensions of a block? If so, there 
are two problems: First, that of explaining what 
occurs at the boundary of two blocks in a region 
where a dislocation ends; and second, that of 
explaining the fact that slip bands often run 
completely around a specimen as though slip 
took place through a plane. 

If neighboring blocks are joined together 
nearly as firmly as parts of individual blocks, as 
we must assume at present, it follows that the 
atoms at the end of a dislocation that terminates 
at a block boundary will encounter resistance to 
motion from the atoms in the neighboring block. 
Moreover, one block will exert shearing stresses 
on its neighbors if a dislocation moves in it and 
not in the others, these stresses being such as to 
oppose the motion. The second type of stress 
obviously will be greatly reduced if dislocations 
are present in the neighboring blocks and move 
more or less simultaneously. A partial reduction 


will occur even if dislocations do not run con- 
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tinuously from one block to another; however, 
it will be greatest if they do. 

The resistance encountered by the end atoms 
during the motion of a dislocation provides the 
most serious obstacle to the assumption that 
fully block 


boundaries. It is a matter of observation® that 


formed dislocations terminate at 
extensive shearing strain can occur in a fraction 3 
of a second for stresses of the same magnitude . 
as those shown in Figs. 19 and 22. According to 
Taylor’s picture of the shear process, this implies 
that a dislocation can move 10~ cm in a fraction 
of a second. Now the time 7 required for an 


atom impeded by barrier to jump from one 





equilibrium position to another is determined in 
order of magnitude by the equation. 


1 
t=—exp (€/k7), 
Vv 


where v~10" sec.~! is the oscillation frequency 





of the atom and € is the activation energy 
required to surmount the barrier. Thus the 
4 
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Fic. 27. Comparison of the observed stress-strain curve 
of aluminum and that obtained from Taylor’s theory with 
L=5.3-10-* cm. The ordinate is expressed in units of 
dynes/cm?. 


minimum time 7° required to jump 10' atomic 
distances is 


1 =10~* exp (€/R7) sec. 
If this is to be less than 1 at room temperature, } 


it is necessary that e be less than 0.5 ev. We 
may conclude either that fully grown dislocations 
*M. Manjoine and A. Nadai, Proc. A. S. T. M. (sum 


mer meeting 1940), have shown for example, that normal 
stress-strain curves may be obtained for tron in 107° sec. 
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run through many blocks, or that they terminate 
at block boundaries and are impeded at the end 
points by very small barriers. 

The fact that slip planes often appear to run 
throughout the crystal specimen is_ readily 
explained if dislocations pass through the entire 
specimen. However, the actual width of slip 
bands is not known because of the finite resolving 
power of micré scopes. Thus the observed results 
would be equally well explained if the disloca- 
tions were much shorter and were distributed 
throughout a lamellar region a hundred or so 
atoms broad parallel to the slip plane. 

2. Why are slip bands separated by finite 
distances? This problem is important not only 
because we are interested in explaining the origin 
of slip bands, but also because the density of 
dislocations in the regions where slip occurs 
should be higher than the mean values used in 
deriving Eq. (2) if slip is highly localized in 
separate planar regions. 

The facts concerning our knowledge of slip 
bands seem to be as follows: 

(a) Slip bands are spaced by distances of the 
order of 1u at ordinary temperatures in highly 
worked seems to be 


material. This conclusion 


agreed upon by a number of 


independent 
investigators.® It also seems to be established as 
a general rule that this spacing increases with 
increasing temperature. 

(b) The spacing between bands is random. 
This point does not seem to be generally accepted 
by all investigators. For example, Straumanis'® 
has claimed that the spacings in zine are integer 
multiples of a constant ; however, extensive and 
apparently trustworthy results of Gough and 
Cox!" indicate that the spacing is random in 
silver. 

The fact that definite slip bands occur at all 
that 


generated with equal ease at all points in the 


seems to imply dislocations cannot be 


specimen. From our present viewpoint, it 1s only 
natural to conclude that the regions of easiest 


nucleation are those where the most. stress 


‘See for example the discussion by E. Schmid and W. 
Boas in Kristallplastizitét (Springer, Berlin, 1935); C. F. 
Elam, The Distortion of Metal Crystals (Oxford University 
Press, 1935). 

10M. Straumanis, Zeits. f. Krist. 83, 29 

H. J. Gough and H. L. Cox, J. 
1931). 


1932). 
Inst. Metals 14, 71 
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magnification occurs. Now we pointed out in 
Part c of this section that only very large stress 
magnifications of the order of 10° 
the nucleation 


independent of temperature. There is no reason 


would be 
sufficient to make process of 
at present for concluding that magnifications of 
this size ordinarily occur; however, even much 
smaller concentrations would have an important 
influence on the rate of nucleation by lowering 
the barrier for formation to a 
height. 


fraction of its 


Two reasonable connections between the slip 
bands and weak spots occur to one: 

(a) In the first place it is possible that each 
slip band is formed by dislocations generated at 
a single spot. In this event, fully formed dis- 
locations would necessarily extend continuously 
through many blocks, and, contrary to Taylor's 
hypothesis, would move through more than one 
block when slip occurs. It is not necessary. to 
assume that all dislocations generated by a single 
weak spot lie in the same atomic plane, for it is 
possible that the region of stress magnification 
extends over a number of planes. However, it is 
necessary to assume that dislocations of opposite 
sign are generated at the weak spot and move in 
opposite directions in a static experiment. This 
implies, in turn, that a lattice of positive and 
negative dislocations of the type postulated by 
that 
hardening would require an 


Taylor cannot form so the problem of 


entirely different 
solution. 

One of the advantages of this picture is that 
it provides a ready explanation of the fact that 
the glide planes do not come closer than 1y in 
spite of the fact that larger and larger stresses 
must be employed to produce slip as work 
hardening proceeds, for it may be assumed that 
all of the weak spots are employed in the fully 
worked solid. 

(b) A number of spots of varving degrees of 
weakness occur on each block and are distributed 
more or less at random throughout the specimen. 
Slip nuclei may be generated at all of these 
weak spots, but a slip band occurs only if a 
number of very weak spots happen to lie very 
nearly in the same plane parallel to the slip plane. 
In this event, slip occurs initially in many parts 


of this planar region causing a concentration of 


175 





stress in the remainder, which, in turn, aids in 


the generation and motion of dislocations in 
these regions. In this picture, the dislocations in 
a given slip plane can be generated at many 
points so that both positive and negative kinds 
will pass in neighboring planes and interfere in 
the manner postulated by Taylor, and Taylor’s 
theory of hardening can be applied. 

In either of these pictures the strain is not 
uniform throughout the specimen, for practically 
no slip occurs in the regions between bands. If 
d is the unknown thickness of a band and D is 
the average spacing between bands, the actual 
strain S, within the slip band is related to the 
average strain S, obtained from direct measure- 
ment of the change in shape of the specimen, 
by the equation 


S,=DS/d. 


Now the S appearing in Eq. (2) should be the 
actual strain in the slip band since this quantity 
determines the density of dislocations in these 
regions, according to Taylor’s picture. Thus we 
seem obliged to rewrite Eq. (2) in the form 


a =kG(ADS/Ld)'. 


D is the order of 1p, but d, being below the 


resolving power of optical microscopes, pre- 


sumably can vary from 0O.1n to an atomic 
distance. In the first of these extremes the 
estimate of L is still about 10-* cm, but in the 


that dislocations 


move through many blocks if Taylor’s theory 


second it must be assumed 
is valid. 

It is possible, of course, that submicroscopic 
changes take place within the regions between 
the obvious slip bands and make the arguments 
leading to the revision of Eq. (2) invalid. In 
this connection, it would be interesting to know 
whether or not the hardness in regions between 
slip bands, as measured by an_ indentation 
method such as the Brinell or Vickers methods, 
increases during slip. If so, it might be concluded 
that submicroscopic slip occurs in these regions. 
that 
strain curves for zinc and aluminum change in 


3. Figures 22 and 23 show the stress 


going from room liquid-air 


temperature, the initial slopes becoming more 


temperature to 


steep. It is remarkable, however, that extensions 


of the order of a hundred percent may be 
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obtained at low temperatures, for we should not 
expect temperature fluctuations to be sufficient 
to allow the dislocations to elongate in a reason- 
able length of time if our picture of their forma- 
tion is correct. Even more remarkable is the fact 
that similar elongations have been produced in 
cadmium at liquid hydrogen temperatures (20°K) 
in a few seconds.‘ Figure 28 shows the manner 
in which the initial slope of the stress-strain 
curve for cadmium varies with temperature. A 
similar curve has been obtained for zinc. It is 
interesting to note that this slope increases very 
rapidly just below room temperature and then 
flattens again below liquid-air temperature. This 
result seems to be relatively insensitive to the 
rate at which the strain occurs and indicates 
that the low temperature mechanism of slip 
contains features essentially different from those 
of the high temperature mechanism. The temper- 
ature changes occurring in cubic crystals such 
as aluminum and copper do not seem to have 
been investigated as thoroughly as those for the 
hexagonal metals. 

If our previous picture of the formation of 
dislocations is correct, we apparently must 
conclude that the slip at very low temperatures 
already formed. In 


is caused by dislocations 


X SLOW EXTENSION ~ 20 MIN 
© RAPID EXTENSION ~~ E 
ii 
. 
x 
= 
a 400 } 
i 
L i 
Fic. 28. The initial slope of the stress-strain curves for 


cadmium as a function of temperature. The ordinate is 
expressed in relative units. (After Boas and Schmid 


this event, the work hardening occurring at low 
temperatures originates in the fact that disloca- 
tions become stuck at lattice imperfections. For 
plausibility, these conclusions evidently require 
that 
temperature remain in spite of annealing. The 
fact 
temperature gives a measurable internal friction 


some of the dislocations formed at room 


that slight handling of crystals at room 
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presumably shows that this assumption is not 
unreasonable. As we shall see in the next part 
of this section, independent estimates indicate 
that there are about 10 dislocations per atomic 
plane in a block in the work-hardened state. 
If one ten-thousandth of this number are present 
in the annealed state, an extension of 100 would 
be attained if each dislocation moved 10? cm. 
Thus we apparently must picture the slip 
process at low temperatures as occurring by the 
motion of a comparatively small number of 
dislocations through large distances. This con- 
clusion evidently favors the viewpoint that block 
boundaries do not seriously limit the motion 
and extension of dislocations, and hence favors 
the possibility that fully formed dislocations 
extend through many blocks. 

The possibility that dislocations are formed 
and are set in motion entirely as a result of the 
action of applied stresses must not be left out 
of mind. A picture of this type evidently is not 
easy to justify on the basis of the model of a 
dislocation we have been employing, but the 
conclusions of the preceding paragraph would 
bear considerable revision if justification were 
found. 

It is interesting to note in passing that the 
existence of a stress field extending for large 
distances about an isolated dislocation in the 
interior of a solid implies that the energy of such 
a dislocation is not all localized in the few atoms 
at its center, or that part is distributed through- 
out large regions of the solid. It is not difficult 
to show™ that a slip nucleus or a line dislocation 
near a surface do not have a similar extended 
stress field so that their energies are localized in 
the atoms near the center of the dislocated 
regions. It follows that as an isolated dislocation 
moves from the surface to the interior the elastic 
energy of the system gradually increases, the 
additional energy being supplied by the applied 
stresses. When many dislocations are present 
and form a stable lattice of the type discussed 
by Taylor, their mutual interaction energy is 
related to a decrease in the extended stress 
rather than to a decrease of the energy of the 

2 This subject is being made the topic of an extensive 


investigation by J. S. Koehler, Bull. Am. Phys. Soc. 15, 
21A (1940). 
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comparatively small number of atoms near the 
center of dislocations. 


f. The energy increase associated with work 
hardening 

It is found experimentally that the internal 

energy of solids is raised during work hardening. 

The most careful measurements of this energy 

have been made by Taylor and Quinney," who 


TABLE III. Values of the maximum energy stored during work 


hardening. (The values are expressed in cal./gram. 


\LUMINUM COPPER IRON NICKEI BRASS 


1.1 0.5 ‘2 0.78 0.49 


found that the maximum increase attainable 
is of the order of magnitude 1 cal./gram for a 
number of metals. This energy may be measured 
in several ways, the most reliable results being 
obtained by measuring the difference between 
work input and heat output during straining 
and by measuring the heat evolved when a 
specimen is heated to temperatures sufficiently 
high to remove the effects of work hardening. 
It is found that the stored energy increases with 
the amount of cold work during the initial 
stages of straining and reaches a saturation value 
when the critical shearing stress approaches a 
constant value, as occurs in compression. This 
limiting value is attained most readily by torsion 
and compression, the values reached during 
tension tests with polycrystals being considerably 
less, presumably because fracture occurs before 
the same degree of strain is reached. The 
circular points in Fig. 29 show the percentage 
ratio of the stored energy to the work done on 
specimens of copper during torsion tests, the 
abscissa being proportional to the torsion strain. 
The triangular points are similar results obtained 
from extension tests. Values of the maximum 
stored energy are given in Table III for several 
metals. 

We naturally must look for the source of this 
increased internal energy in the increased imper- 
fection of the crystals. The two most likely 


8 (4. I. Taylor and H. Quinney, Proc. Roy. Soc. 143, 307 
(1934); 163, 157 (1937). Earlier work in tension was 
carried out by G. I. Taylor and W. Farren, Proc. Roy. Soc. 
107, 422 (1925). 





possibilities are the distortion of blocks, which 


is responsible for broadening of the x-ray 
diffraction lines of crystals, as we saw in the 
Part b of this section, and the number of disloca- 


tions. Cogliotti and Sachs," and later Boas,'® 
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Fic. 29. The percentage ratio of the stored energy to the 
mechanical work done on copper specimens during cold 
work, (After Taylor and Quinney). The triangles refer to 
tension tests and the circles to torsion tests in which much 
greater straining is possible. The abscissa is proportional 
to the strain. The falling portion of the curve implies the 
stored energy reaches a saturation value. 


pointed out that the distortion of the blocks is 
not sufficiently large to account for more than a 
very small fraction of the stored energy. Thus 

Dehlinger and Kochendérfer'® have found that 
the percentage distortion of lattice constant in 
the most highly cold-worked metals is of the 
order of 0.1, which corresponds to an elastic 
energy of about 10~* cal./gram. Apparently we 
must conclude that the stored energy is to be 
associated with dislocations. This assumption is 
reasonable since the number of dislocations 
should increase with cold work according to 
Taylor's theory of hardening. 

. If we make this assumption we may readily 
estimate the density of dislocated atoms in an 
approximate fashion. We saw previously that 
it is reasonable to suppose that the energy of a 
dislocation per atomic distance along its length 

“VV. Cogliotti and G. Sachs, Zeits. f. Physik 74, 647 
(1932). 

'% W. Boas, Zeits. f. Krist. 96, 214 (1937); 97, 354 (1937). 
Other papers on this topic are: G. R. Stibitz, Phys. Rev. 
52, 619 (1937); F. E. Haworth, Phys. Rev. 52, 613 (1937); 


G. W. Brindley and P. Ridley, Proc. Rov. Soc. 50, 501 
(1938). 


This work is surveyed in a review article by U. 


Dehlinger and A. Kochendérfer, Zeits. f. Metall. 31, 231 
(1939). 
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is about 1 ev, or 10-" erg. Since the stored 
energy is of the order of magnitude 10° ergs per 
mol (one mol contains 6.10”* atoms), it follows 
that about one atom in a thousand is near the 
center of a dislocation in a highly worked 
specimen. This is in good qualitative agreement 
with the number that we would estimate directly 
from Taylor’s picture of the hardening process. 
If an average dislocation moves a distance of 
10* atoms, as is assumed in Taylor’s theory, and 
if the strain at maximum hardening is about 
ten, there must be ten dislocations per atomic 
plane in each block in the hardened crystal. 
Hence about one atom in a thousand must be 
near the center of a dislocation." 

Excellent corroborating evidence for the con- 
clusion that the energy of cold work is to be 
associated with defects extending over atomic 
dimensions has apparently been obtained by 
Hengstenberg and Mark and by Brindley'* and 
his co-workers. As we pointed out in the previous 
installment (Part 6 of this section) the principal 
influence on the x-ray diffraction patterns of 
defects of atomic size is a decrease in the in- 
tensity of the diffracted line because of incoherent 
scattering by the defects. Brindley and _ co- 
workers have found that the intensity of diffrac- 
tion lines of metals actually decreases with 
increased cold working and Boas'® has shown 
that the elastic energies associated with the 
atomic displacements derived from their results 
are of the same order of magnitude as the 
observed stored energies. This does not prove 
that the defects are dislocations, of course, but 
provides additional support for the picture that 
defects of atomic size play a primary role in 
the hardening process. 

An interesting and critical experiment on the 
mechanism of slip at low and high temperatures 


174 fraction of one dislocated atom per thousand 
normal atoms corresponds to a density of dislocation lines 
of about 10" per cm?. W. F. Brown, Phys. Rev. 58, 736 

1940), has derived the same density from an investigation 
of the manner in which the magnetization of highly worked 
materials approaches saturation. The writers are indebted 
to Professor Brown and Dr. R. Smoluchowski for discus- 
sions of this point. 

‘8H. Hengstenberg and H. Mark, Zeits. f. Physik 61, 
435 (1930); Zeits. f. Elektrochemie 37, 524 (1931); G. W 
Brindley and F. W. Spiers, Phil. Mag. 20, 882 (1933); 
G. W. Brindley and P. Ridley, Proc. Roy. Soc. 50, 501 

1938); 51, 432 (1939); Phil. Mag. 21, 778 (1936); Proc. 
Roy. Soc. 52, 117 (1940). 
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would be provided by measurements of the 
energy stored during cold work at very low 
temperatures. If the supposition, made in the 
previous part of this section, that a compara- 
tively small number of dislocations are responsi- 
ble for slip at low temperatures, as compared 
with the situation assumed by Taylor at high 
temperatures, we should expect striking differ- 
ences in the ratio of stored energy to work done 
in low temperature and high temperature cases. 

Simple estimates, based on a comparison of 
the work done during hardening and thermal 
data for metals, show that sufficient heat is 
produced during extensive hardening to melt 
about 5 percent of the material in the cases 
listed in Table III. Thus if straining actually 
occurs in very narrow slip bands, it is conceivable 
that local melting would be produced if the strain 
occurred in a time short compared with that for 
thermal equilibrium to be reached. Unfortu- 
nately, the time for equilibrium appears to be 
too short for melting actually to be observable. 
For example, if the slip bands are 500A wide, 
the time required for thermal equilibrium is of 
the order of 107! second. 


g. Resoftening 


The effect of work hardening may be undone 
by appropriate heat treatment. Apparently this 
resoftening takes place to some extent at all 
temperatures above absolute zero, but it may be 
greatly accelerated by heating. Thus zinc re- 
softens in the course of several days at room 
temperature, whereas the rate is negligibly 
small at liquid-air temperatures. It is common 
in practical work to designate the lowest temper- 
ature at which softening takes place at a_practi- 
cally interesting rate, such as a day or a week, 
the recovery temperature. This is defined to 
within an ambiguity of twenty degrees or so 
depending upon the range of time in which 
recovery is expected. Just above the temperature 
at which resoftening is first observable in a 
practical sense, a change in crystalline texture 
of worked specimens occurs. This effect, which 
is known as recrystallization, will be discussed 
later in connection with polycrystals. Since the 
temperatures at which recrystallization begins 
have been tabulated for a large number of 
materials and since this temperature marks an 
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upper limit to a range of about thirty degrees or 
so in which resoftening can occur without change 
in texture, we list the values'® of the recrystalli- 
zation temperature for the polycrystalline form 
of several metals in Table IV. 


ABLE IV. Recrystallization temperatures of several metals. 
These temperatures mark the upper limit of the range in 
which resoftening can occur without change in texture. The 
values are expressed in degrees centigrade and are only 
approximate. ) 


aluminum 150 nickel 620 
cadmium 50 platinum 470 
copper 200 silver 200 
gold 200 tantalum 1020 
iron 450 tin 0) 
lead 0 tungsten 1210 


magnesium 150 zim 15 


Karnop and Sachs*® have shown in the cases 
of aluminum and copper that the recrystalliza- 
tion temperature is higher in single crystals than 
in polycrystals, for which the values in Table IV 
are valid. Thus recrystallization occurs only at 
about 600°C in single crystals of aluminum. 
Resoftening, however, occurs at far lower temper- 
atures, in the range below the recrystallization 
temperature for highly worked polycrystals given 
in Table IV. It is to be noted that in less severely 
cold-worked specimens recrystallization takes 
place at temperatures higher than those given in 
Table IV. 

Figure 30 shows stress-strain curves for zinc 
obtained at room temperature under two differ- 
ent sets of experimental conditions. In the case 
corresponding to Fig. 30a the specimen was 
given a succession of fifty-percent extensions, 
the load removed for a “rest period” of about 
one minute between the extensions. In the 
second case, the extensions were identical, but 
the rest period was increased to a day. It is 
clear that the first period of rest was not suffi- 
ciently long for very extensive resoftening to 
occur, whereas softening was almost complete 
after the second. The dotted curves in the 
figures show the vield stress calculated from th« 
original critical shearing stress of the specimen. 
The vield stress varies as a result of reorientation 

'° Taken from the article by R. S. Archer, A. S. M. 


Handbook, (1939), p. 198. 


20 R. Karnop and G. Sachs, Zeits. f. Physik 42, 
1927). 
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of the specimen during extension. The observed 


value agrees closely with this calculated one at 
the start of each run in case b. 

The phenomena of resoftening would seem to 
be closely connected with the time dependence 
of the type of internal friction discussed in Part d 
of this section (previous installment), for the 
variations take place in comparable periods of 
time in both cases. On the basis of dislocation 
theory, we are tempted to conclude that above 
the recovery temperature dislocations introduced 
into the crystals during straining diffuse out and 
thereby reduce the solid to its virgin state. 
This exodus should be accompanied by a 
decrease in internal friction, at least in the later 
stages when the number of dislocations has 
become so small that they do not interfere with 
one another appreciably, since there are fewer 
centers to dissipate the vibrational energy. Such 
a decrease actually occurs, as we saw previously 
in the case of zinc. 

It seems worthwhile to speculate on some of 
the details of the process by which dislocations 
diffuse out of crystals. In the first place, it is 
possible that dislocations move about primarily 
as a result of thermal agitation, just as atoms 
diffuse as a result of thermal motion. On the 
other hand it is possible that the motion is 
produced as a result of stress acting on the 
dislocations and that thermal fluctuations play 
a role only insofar as they furnish the activation 
energy required to get the dislocation past 
obstacles that impede very small regions along 
the line of the dislocation. It will be recalled 
that in Part ¢ of this section we inferred that 
temperature plays a similar role in the process 
of elongation of a dislocation. Concerning the 
possible origin of stresses in the second case, we 
noted at the end of Part e of this section that 
the elastic energy of a solid increases as a 
dislocation moves from the surface to the interior. 
This implies that the solid exerts stresses on the 
dislocation of such a sign as to force it to the 
surface.” The force on an isolated dislocation 
may be estimated roughly in the following way. 
If we have a dislocation at a distance x from the 
surface, the requisite boundary conditions may 
be approximately satisfied by assuming that 
there is a dislocation of opposite sign at a 
distance —x on the other side of the surface. 
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This ‘“‘image”’ dislocation exerts a shearing stress 
of magnitude 
Gr/ 24x (3) 


(see footnote 7 for notation) on the given 
dislocation. The quantity (3) is about 10° 
dynes/cm? for a dislocation midway in a block 
and should be sufficient to drive an isolated 
dislocation out of an otherwise perfect crystal 
if our views concerning the ease with which 
dislocations may move is correct. Naturally the 
required stress is larger when there are many 
dislocations present; however, we may visualize 
the process as taking place in such a way that 
dislocations near the surface which are under 
greater stress, according to Eq. (3), move out 
first and the others then follow in succession. 

Equation (3) is valid, of course, in the case 
in which dislocations are produced at block 
boundaries, or in the interior of blocks, in pairs 
of opposite sign which move in opposite direc- 
tions. In this case the image dislocation is 
replaced by a real one. 

Thus estimates based on (3) indicate that 
dislocations will tend to move out of blocks when 
the temperature is sufficiently high to provide 
the activation energy for getting by localized 
obstacles. The assumption that this actually is 
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Fic. 30. Stress-strain curves for zinc at room tempera- 
ture in cases in which the specimen was allowed to rest 
during successive elongations. In case (a) the rest time is 
about a minute and in case (b) it is 24 hours. (After 
Haase and Schmid.) 


the mechanism requires, however, that we accept 
some of the details of Tay lor’s picture concerning 
the distances that dislocations move from the 
places where they are generated and it seems 
worthwhile to consider further the possibility 
that they would wander about as a result of 
thermal agitation even if no stresses were present. 
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In the first place, let us assume that the 
dislocation runs through only a single block. If 
we use the model of a dislocation considered in 
Part c (Figs. 8, 9, 10 and 12), we must conclude 
that one part of a given line dislocation cannot 
readily move ahead of the remainder since a 
large amount of energy would be required to 
break the line into two parts. Now, exclusive of 
any friction the atoms at the end of the disloca- 
tion may encounter when the dislocation moves, 
the activation energy required for motion of a 
line 10‘ atoms long from one equilibrium position 
to another will be a small fraction of an electron 
volt if our assumption concerning the low critical 
shearing stress is correct. Thus the only obstacle 
in the way of rapid thermal diffusion of disloca- 
tions, aside from end friction, is the question of 
whether it is probable that all atoms in the line 
can move in phase so as to jump simultaneously 
from 


one 


equilibrium position to the next. 


Apparently there is no difficulty arising from 
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Fic. 31. The critical shearing stress of a sequence of 


silver-gold alloys as a function of composition. These 
metals are miscible in all proportions. (After Sachs and 


Weerts.) 


this source, for if one atom in the dislocation 
gains sufficient energy to take the entire disloca- 
tion to the next equilibrium position, and is 
moving in the proper direction, it will be able to 
drag the others with it the 


because of close 
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coupling between atoms in the dislocation. In 


this respect a dislocation 10‘ atoms long is 


analogous to a particle undergoing Brownian 
motion in air or in a liquid. The only difference 
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Fic. 32. The critical shearing stress of a sequence of 
zinc-cadmium and zinc-tin alloys. At 250°C somewhat 
more than 1.5 weight percent of cadmium is soluble in 
zinc, whereas at the same temperature less than 0.1 weight 
percent of tin will dissolve in zinc. These results illustrate 
the fact that soluble impurities have a greater influence 
than insoluble ones. (After Rosbaud and Schmid.) 


between the two cases appears to lie in the fact 
that a dislocation requires an activation energy 
for motion between equilibrium positions. 

Naturally, if the barrier encountered by end 
atoms is very large, the rate of diffusion will be 
cut. However, we have already seen from other 
evidence, discussed in Part e, that the assumption 
that dislocations terminate at block boundaries 
requires that the barriers at these points be 
very low. 

If dislocations are very long, running through 
many blocks, their natural rate of diffusion will 
be correspondingly decreased; however, it does 
not seem 


to definite 


conclusions concerning their mobility without 


permissible draw any 
more quantitative information concerning their 


properties. 
h. The hardening influence of impurities 


The of on the 
critical shearing stress of metals is shown in 
Figs. 31, 32 and 33. The first of these?! shows 


influence impurity content 


"1G. Sachs and J. Weerts, Zeits. f. Physik 62, 473 (1930). 
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the variation in critical shearing stress in a 


sequence of silver-gold alloys that were carefully 
annealed before measurements were made. The 
striking features of these results are the fact 
that the hardness of solid solutions is greater 
than that of the pure metals and that the 
variation in hardness is a simple function of 
composition. This increased hardness of alloys 
as compared to pure metals is, of course, partly 
responsible for their wide use. The second of 
these figures” illustrates the relative increase in 
hardness of zinc-cadmium and zinc-tin alloys. 
Cadmium is soluble at 250°C to the extent of 
about 1.5 weight percent in the first of these 
alloys, whereas less than 0.1 weight percent of 
tin is soluble in zinc. This case provides an 
illustration of the general rule that the hardening 
effect of a soluble impurity is greater than that 
of an insoluble one. Figure 33 shows the effect 
of very slight traces of silver impurity on the 
critical shearing stress of extremely pure mercury. 
The results were obtained by Greenland® on 
single crystals of mercury that was carefully 
purified by distillation. It is very interesting to 
note that these results indicate that impurity 
contents as low as 10-7 or 10-* have a measur- 
able effect on the hardness. 

In this connection, it is interesting to note the 
influence on hardness of different atoms present 
in the same atomic percentage. Figure 34 shows** 
the Brinell hardness of a number of copper 
alloys, per atomic percent of the dissolved metal, 
as a function of the relative difference in sizes of 
the dissolved atom and the copper atom. 
Although the Brinell hardness is not a direct 
measure of the critical shearing stress, we may 
assume that the two increase together. Figure 34 
then indicates that dissolved atoms have a 
greater hardening influence the greater the 
fractional difference in the diameters of the 
solvent and solute atoms. An approximately 
linear relation of the type of Fig. 34 is also 
valid in a sequence of lead alloys.” 

The interpretation of the hardening influence 
of impurities or alloying elements on the basis of 
the theory of dislocations is not entirely apparent, 

2 P. Rosbaud and E. Schmid, Zeits. f. Physik 32, 197 
1925). 


22K. M. Greenland, Proc. Roy. Soc. 163, 28 (1937). 
* A. L. Norbury, Trans. Faraday Soc. 19, 586 (1924). 
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particularly in cases such as mercury in which 
very small atomic fractions seem to be effective. 
Since soluble impurities have a greater influence 
than insoluble ones, we may conclude that the 
effect occurs primarily within the volume of the 
crystals rather than at their surface or at block 
boundaries. Moreover, since there is a simple 
relation between hardening and atomic size, it 
seems reasonable to suppose that the effect is 
related to the distortion of the lattice by the 
alloying atoms. Now it is clear that these atoms 
will compress or expand the surrounding lattice 
in a way depending upon whether their size is 
larger or smaller than that of the normal atoms. 
In either case, the impurity atom will be the 
center of a more or less spherically symmetrical 
stress field. This radially symmetric compres- 


sional or tensional field varies with the distance 


~ 


r from the center of the impurity atom in the 
manner”? 


o-=4Gfa®/r’, o9= —4Gfa*/r 4) 
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Fic. 33. The critical shearing stress of single crystals of 
very pure mercury containing small amounts of silver. 
\fter Greenland.) The results were obtained at about 
—60°C, which is only about 23° below the freezing point. 
Work hardening does not occur at this temperature so 
that the specimens were not sensitive to previous handling. 


where ¢, is the stress in the radial direction and 
os is that in the angular direction, G is the shear 
modulus, which is about 10" dynes/cm’, f is the 
relative difference in radii of the solute and 


solvent atoms, and a is the atomic radius, which 


** See for example, A. E. H. Love, The Mathematical 
Theory of Elasticity (Cambridge University Press, 1927), 
142. The present discussion foll ‘losely tl l | 
p. 142. Lhe present discussion tollows closely that employer 
by Mott and Nabarro in connection with precipitation 
hardening, discussed in the next part ol this section. 
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is usually about 1.5-10~-* cm. The shearing stress 
in any plane may be computed from this by 
projecting the tension stress in the plane. In 
passing by this stressed region, a dislocation will 
be acted upon by a shearing stress having an 
azimuthal variation of the type shown in Fig. 35, 
that is, the dislocation will be retarded on one 
side of the 90° azimuth and will be accelerated 
on the other. In either case, it will experience 
greater difficulty in passing the atom, so that 
the critical shearing stress should increase, just 
as in the case in which there are other dislocations 
present to impede the ease of motion. 

This picture is in qualitative agreement with 
the fact that the hardness is a maximum for 
50-50 composition in an ideal substitutional 
alloy, for the amount of inhomogeneity should 
then be a maximum, at least in disordered 
substitutional alloys. It seems to face some 
difficulty, however, in interpreting the hardening 
in very pure mercury, although the following 
argument indicates that such an. interpretation 
may be within the scope of the theory. 
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DIFFERENCE IN SIZE OF ATOMS (RELATIVE) 


Fic. 34. Relation between the relative increase in hard- 
ness and the fractional difference in sizes of the solvent and 
solute atoms for a number of copper alloys containing 0.1 
atom percent of alloying metal. (After Norbury.) The 
abscissa is the difference in atomic volume of solvent and 
solute. This difference is proportional to the difference in 
atomic radii. 


To begin with let us consider the case in 
which the crystal is at absolute zero of tempera- 
ture so that thermal fluctuations do not play a 
role. The radius of the cylinder about a line 


dislocation 10' atoms long in which there is one 
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silver atom when the concentration of silver 
atoms is about 1077 is 50A. Thus at an average 


position each dislocation encounters a shearing 





Fic. 35. The variation in shearing stress with azimuthal 
angle about an impurity atom. The shearing stress is 
measured in a system of parallel planes and at distances 
equally spaced from the impurity atom. 


stress of about 10° dynes cm®* extending over a 
region approximately fifty atoms square. Unless 
the applied stress is sufficient to compensate 
for this stress arising from the impurity atom, 
the dislocation will become stuck. It is interest- 
ing to note that the observed critical shearing 
stress in the relatively pure mercury is of this 
order of magnitude. 

Now the work required to push the dislocation 
past this difficult spot may be shown to be of 
the order of 10-" erg or 0.1 ev. This barrier is 
sufficiently small that we should expect the 
dislocation to diffuse past it readily even at dry 
ice temperatures. It is possible, however, that 
speeds of migration of the type required for 
observation of slip in Greenland’s measuring 
apparatus required that the applied stress be 
sufficient to neutralize the retarding stress 
arising from impurities. 

If dislocations are considerably longer than 
10-* cm, the average distance from an impurity 
atom to a dislocation should be much smaller 
than 50 atom distances and the resisting stress 
should increase. Measurements of the critical 
shearing stress of the mercury crystals at 
temperatures sufficiently low that thermal acti- 
vation is unlikely should provide a direct 
measure of the mean stress encountered by 
dislocations and hence should furnish an indirect 
method of estimating their length. 





One of 


reasoning is 


the consequences of the foregoing 


that the critical shearing stress 
should be proportional to the three-halves power 
of the 


concentration of impurities for small 


concentrations. This may be seen from the fact 
that the square of the radius of the cylinder 
about a dislocation containing one impurity 


atom varies inversely as the concentration of 
the impurity atoms, whereas the stress about 
the impurity atom varies inversely as the cube 


(4). It is 


impossible to tell from Figs. 31 and 32: whether 


of the distance, according to Eq. 


or not a law of this type is valid for very small 


concentrations. Greenland’s results, shown in 


Fig. 33, seem to require a more complicated 
This work, 


however, is open to the objection that silver is 


dependence upon concentration. 


only slightly soluble in mercury and it is not 


proved that the silver atoms are randomly 


distributed in the crystals on which the measure- 
ments were made. 

In addition, it is interesting to note that the 
previous reasoning predicts that the critical 
the 


relative difference in atomic radii of the impurity 


shearing stress should vary linearly with 


and solvent atoms. This is in reasonable accord 
with Fig. 34 and with results for lead described 


by Norbury. 
i. Precipitation hardening 


In many alloy svstems the solubility limits of 
given phases decrease with decreasing tempera- 


ture. Thus copper is soluble in aluminum to the 


extent of 4.1 percent at 500°C, whereas the 
solubility is less than 0.5 percent at room 
temperature. For this reason perfectly homo- 


geneous crystals of aluminum containing several 
percent of copper are stable at elevated tempera- 
tures, but become thermodynamically unstable 
when cooled near room temperature. It is found 
in many such cases that precipitation of the 
dissolved metal proceeds in the course of hours 
this 


or days at that 


room temperature and 


precipitation is accompanied by hardening of 
the alloy. The dynamics of this precipitation 
process are fairly complex and will not be 
reviewed here. 

An analvsis of the theoretical factors entering 


into the problem have been given by Mott and 
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Nabarro.*® They point out that the precipitated 
material usually has an atomic volume different 
from that of the homogeneous alloy, so that the 
lattice becomes strained in the vicinity of the 
precipitate (see Fig. 36). This straining implies 
that there are regions about the precipitate 
where the shearing stress is unfavorable for the 
passage of dislocations, just as in the case of 
impurity atoms discussed in the previous part of 
this section. In fact, Eqs. (4) may be employed 
to describe the stress field about a_ spherical 
particle of precipitated material. Using reasoning 
of this Mott 


that in the case of spherical particles, the critical 


type, and Nabarro have shown 


shearing stress is given roughly by the equation 
o=Gfr, 

where G is the shear modulus for the matrix, f is 

which the linear dimensions of 


the 
compressing action of the matrix, and r is the 


the fraction by 


the precipitated particle is changed by 


ratio of amount of precipitate to matrix. 


j- Multiple gliding 

When slip is taking place in a given set of 
parallel planes in a crystal of high symmetry 
theré usually are symmetrically equivalent planes 
which may be regarded as latent slip planes. 
Thus if slip is occurring in one set of (111) planes 
in a face-centered cubic metal there are three 
latent systems of easiest slip. During the process 
of slip it. frequently happens*’ that, as a result 
of reorientation of the specimen, the shearing 


stress in one of the latent slip planes becomes 


larger than that in the actual slip planes. 
Accumulated evidence for face-centered cubic 
crystals of the pure metals indicates that 


hardening occurs in the latent slip planes at 
about the same rate as in the actual slip planes, 
for gliding will occur in the first svstem just as 
soon as the shearing stress becomes larger than 
in the actual slip planes. On the other hand, 
the (101) 
indicates that the latent slip planes harden more 


evidence for slip in planes of tin 


rapidly than the actual planes, that is, 


the after 


slip 


occurs in latent system only the 
*N. F. Mott and N. F. R. Nabarro, Proc. Phys. Soc. 
52, 86 (1940); N. F. R. Nabarro, Proc. Phys. Soc. 52, 90 
1940); Proc. Roy. Soc. 175, 519 (1940). 
27 The topic of multiple gliding is discussed in the books 


by Elam and by Schmid and Boas (see reference 9). 
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shearing stress in these planes has become much 


larger than that in the actual system. Similar 


results have been found by Elam?’ in single 
crystals of copper-zinc and copper-aluminum 
alloys. 

The relationship between hardening in actual 
and latent slip planes is not very simple when 
viewed in terms of the model of dislocations we 
have been considering. To begin with, let us 
consider a hypothetical cubic case in which the 
planes of easiest slip are three systems equivalent 
to (100), and let us assume that the specimen 
contains a uniform distribution of dislocations 
which can move in the x direction and which 
extend along the y direction, so that when they 
move the slip planes are x—y planes and the slip 
direction is parallel to the x axis. Now because 
of the relation ¢;;=0;; between the components 
that 


encountered by 


of shear stress, it follows the average 


shearing stress one of these 
dislocations as a result of the presence of the 
others is the same as that encountered by a 
dislocation which moves in the direction of the 
z axis and extends along the y axis. Thus, using 
Taylor’s picture of the hardening mechanism, 
we should expect 


hardening to proceed as 


rapidly in the y—z planes for slip in the z 
directions as in the x—y planes for slip in the x 
direction, after slip has started in the second 
system. On the other hand, dislocations which 
and extend in the x 


move in the z direction 


direction should encounter no shearing stress 
because of the presence of the other. types of 
dislocations, and we should not expect the x—2z 
slip planes to harden at the same rate as the 
others. Although this is not a practical case, it 
does indicate that we shall meet difficulties in 
the 
latent planes in the face-centered cubic crystals 


attempting to understand hardening of 


in terms of Taylor’s mechanism and line disloca- 
tions of the type we have been treating. In fact 
it would appear that the observation that all 
planes harden equally is explicable only if the 
distortion induced by the dislocations associated 
with the primary slip process has the same cubic 
symmetry as the lattice; this does not seem to 
be conceivable in terms of the concepts we have 
been using. 

It is possible, of course, that the rule stated 
above for the hardening of latent slip planes in 
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(a) 


(b) 


Fic. 36. Schematic representation of the strain about an 
inclusion of precipitate in a precipitation hardened alloy 
After Mott.) (a) Situation occurring if the atomic volume 
of the precipitate is the same as that of the solvent. 
b) Actual situation at equilibrium in a case in which the 
atomic volume of the precipitate is larger than that of the 
solvent. 


pure face-centered cubic crystals is based on 
insufficient evidence and that all latent planes 
do not harden equally. In this event, additional 
experimental work on the point would be of 
great value. 

A discussion of the cases of pure tin and the 
copper alloys, in which the latent planes harden 
more rapidly than the others, does not seem in 
order until the more simple one of pure cubi 


metals is understood. 


k. Asterism 


Although a 


strained by 


crystal that has been uniformly 


slip remains a single crystal, in 


principle, its x-ray diffraction patterns show 


striking changes that are indicative of internal 
the 
Laue x-ray pattern, which is obtained by sending 


distortion. In particular, it is found that 


a beam of continuous 


through the 
st reaks or 


instead of sharp spots, as in an unstrained single 


radiation 


system, consists of radial bands 
crystal. This streaking is generally called aster 
ism.2> The appearance of asterism implies that 
the crystallographic planes which reflect 


the 


par 


ticular wave-lengths§ of incident beam, 


and thereby give Laue spots, are bent as a 


result of the straining. It is not possible to sa\ 
whether or not this bending is accompanied by 
fragmentation. 

It has been 


pointed out by a number of 


"See tor exan ple, Schmid ind Boas and Elam, rele 


ence 9: W. G. Burgers, Proc. Int. Conf. in Physics 139 
1934): | N. da C. Andrade Pro Phys. Sor -. i 
1940) 
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Fic. 37. Origin of asterism on basis of Burger's picture. 
The material between dislocations is distorted as a result 
of the stresses about dislocations 


workers that the distortion of the Laue pattern 
can be explained by assuming that the bending 
of planes takes place about an axis lying in the 
slip plane and normal to the slip direction, the 
bending varying statistically 


angle of over a 


finite range. However, there is no conclusive 
evidence concerning the size of the bent regions 
or their positions relative to the planes on which 
slip has occurred. 

Barrett and Levenson*® have pointed out that 
cold-worked crystals frequently exhibit narrow 
bands on either side of which the crystallographic 


orientation is slightly different. They have 


27 C._.S,. Barrett, Trans. A. I. M. E. 135, 296 (1939); C.S 
) 


Barrett and L. H. Levenson, Trans. A. 1. M. EF 135, 327 
1939), 
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suggested that the disorientation associated with 


those distortion bands is the seat of asterism. 
Burgers*® has suggested that the warping of 
the lattice is directly connected with the presence 
of dislocations in a manner such as that shown 
in Fig. 37. This hypothesis awaits more extensive 
confirmation. In disagreement with it appears 
to be the fact, observed by Karnop and Sachs*° 
that 
extensive resoftening does not appreciably affect 


in aluminum strained by 15. percent, 
the asterism. 

It was recently shown by Andrade and Chow*! 
that the asterism changes its character when the 
crystal is heated sufficiently, for the continuous 
streaks appear to coagulate into spots extending 
along the line of the streak. This result seems to 
imply the 


material changes its orientation to one of several 


that, as a result of heating, bent 
arrangements bearing a fixed relationship relative 
to the bulk of the ecrvstal. Since the original 
rotated material is presumably the most highly 
distorted in 


the crystal, it is reasonable to 


suppose that it is thermodynamically unstable 
and that the reorientation accompanying heating 
is associated with a transition to a more stable 
state. It is also natural to suggest that the newly 
the nuclei 

Part g), 


but this possibility has not yet been subject to 


oriented material is identical with 


responsible for recrystallization (see 


careful experimental investigation. 


J. M. Burgers, Proc. Phys. Soc. 52, 23 (1940). 
' | N. da C. Andrade and = how, Proc Roy SOC 
175, 290 (1940). 
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Rockefeller Foundation Grant 


rhe following statement is an excerpt from the Minutes 
of the November meeting of the Executive Committee of 
the Governing Board of the Americar Institute of Physics: 

‘The American Institute of Physics has for several years 
been greatly indebted to the Rockefeller Foundation for 
encouragement and financial assistance in developing a 
well-rounded and normally self-sustaining publication pro 
gram in the field of physics. A drastic decline in foreign 
dues and subscription income caused by the war in Europe 
has now threatened to cancel the gains which have been 
made. In response to the Institute’s appeal on this account, 
the Foundation has generously granted the sum of $20,000 
to assist the Institute to adjust the publication program to 
the new circumstances. 

‘Therefore, be it resolved: That the Executive Committe 
express for the Institute its full appreciation of this timely 
and generous assistance and convey to the Rockefeller 


foundation an assurance of the gratitude of physicists.”’ 


New Electron Microscope 


Details of a more effective and simplified model of the 
RCA electron microscope were discussed recently before a 
convention of the Institute of Radio Engineers in New 
York City, by James Hillier of the RCA Research Labora 
tories. A section of the new model is shown on the cover of 
this issue; Dr. V. K. Zworykin, under whose direction the 
instrument was designed and constructed, is shown stand 
ing beside it; while Jame s Hillier, co-deve loper, sits at the 
controls. The electron microscope is capable of magnifying 
objects up to 100,000 times and has therefore far-reaching 
possibilities in virtually every field of medical and in 
dustrial research. RCA’s first electron microscope, an 
nounced early last year, was designed primarily for use by 
electrical experts and physicists; the new model is smallet 
and better adapted for use in general medical and chemical 
laboratories since it may be oper ited by any competent 
laboratory worker. 

Despite its radical design improvements, Mr. Hillier 
said, the new microscope vields nothing in operating 
efficiency to the earlier model. The new instrument, 
housed in a steel cabinet, is completely self-contained and 
is ready for operation when plugged into an ordinar\ 
light socket, according to Mr. Hillier. It is not affected by 
magnetic fields set up by other electrical equipment, as 
was the earlier model, nor is it disturbed by ordinar: 
room vibrations. An important feature is an ingenious 
“submarine lock” arrangement for changing specimens in 
the evacuated chamber of th instrument without ad 
mitting air. A number of viewing windows are provided 
through which several persons can study specimens at 


the same time. 
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The RCA electron microscope operates very much like 
an ordinary light microscope, except that instead of glass 
lenses and light, it utilizes a stream of electrons prope lled 
at high speed through a series of magnetic coils, M1 
Hillier explained. ‘Commercial possibilities for the ele 
tron microscope, even in such an all-important field as 
national defense (where it may be used to improve ma 
terials and processes) led to the development of this 
simpler instrument,” he said. “We are now engaged in 
producing the new microscope for use in leading educa 
tional and commercial research laboratories throughout 


the country 


* 


Book on ‘‘Temperature—Its Measurement and Control’’ 


It was announced last year that a book on ‘‘Temperature 

Its Measurement and Control in Science and Industry” 
Was in preparation to make permanently available the 
valuable papers presented at a symposium on this subject 
sponsored by the American Institute of Physics with the 
cooperation of the National Bureau of Standards and the 
National Research Council, held in New York. This 
monumental volume will be ready for distribution about 
the middle of February. The original papers have been 
carefully edited by authorities in their fields. They have 
also been rearranged, supplemented and exhaustively 
indexed. All phases of temperature are covered in great 
detail. By sper ial arrangement, the book is to be published 
by the Reinhold Publishing Corporation, 330 West 42nd 
Street, New York City, at a list price of $11.00 per copy 
\ valuable Appendix of 25 reference tables is also included 
and this may be purchased separately bound for $1.00 
per copy. With introductory material and indexes, the 
volume runs to 1375 pages. It will be especially valuabl 
for chemists, teachers, physicists, biologists, medical men, 
metallurgists, geologists, electrical, mechanical and _ re 


frigerating engineers, and plant-operating officials. 
* 


Honor Parole 


\ccording to a wireless dispatch to The New York Times, 
Dr. Paul Langevin, formerly Professor of Experimental 
Physics at the College de France, who had been detained 
in the Cherche Midi prison in Paris, has received an honot 
parole. His release was ordered for reasons of health. 
Professor Langevin was deprived of his professorship 


because of his politi al activities. scence. 
* 
Holder of 52 Patents Cited 


Dr. C. Guy Suits, 35-year-old physicist who was named 
assistant to the Director of the General Electric Research 
Laboratory in October, has been chosen as one of the 
outstanding young men of 1940 in the annual selection of 
Durward Howes, editor of America’s Young Men, a 
biographical work. Dr. Suits, a graduate of the University 
of Wisconsin, has 52 patents and credit for creating the 
hottest flames ever made by man. These were electric 


ircs twice as hot as the surface of the sun 





Dublin Institute for Advanced Studies 


The founding of the Dublin Institute for Advanced 
Studies, by an act of the Oireachtas, in October, was 
announced in a recent issue of Nature. For the present it 
consists of two schools, the School of Celtic Studies and 
the School of Theoretical Physics. The Institute will 
provide facilities for advanced studies and research in 
special branches of knowledge, according to the announces 
ment, and for the publication of the results of such studies, 
irrespective of whether or not they have originated from 
the Institute. Moderate fees will be « harged, but a limited 
number of scholarships are available 

Che scope ol the School of Theoretical Physics is de- 
scribed as ‘‘the investigation of the mathematical principles 
of natural philosophy and their application to the sciences 
in which they obtain.”’ The first senior professor appointed 
to the School of Theoretical Physics is Professor E. 


Schrédinger, formerly of the University of Graz. 
* 


University of Michigan Honors Retiring Professor 


Professor Harrison M. Randall, Chairman of the De- 
partment of Physics of the University of Michigan has 
retired after forty vears of service in the University, 
according to a recent note in Science. In his honor East 
Physics Hall is to be known as the Harrison M. Randall 
Laboratory of Physics, the regents of the University have 
announced. Dr. Randall has been appointed Henry Russel 
lecturer for 1940-41. This award goes to the faculty mem- 
ber judged to have attained the highest scholastic distine- 
tion during the year. 


* 
Many-Valued Logics 


Seven years ago a physicist, Dr. F. Zwicky of the Cali 
fornia Institute of Technology, declared that he thought it 
would be easier to explain the atom if he could use a many 
valued logic instead of the usual two-valued logic. By a 
many-valued logic he meant one which denies the principl 
of reductio ad absurdum. This new logic states that, in 
stead of the two possibilities, ‘‘so’’ and “not so,”’ there are 
n possibilities. When is small, there might be a chance of 
saying what the possibilities are. For instance, when n is 
six, the possibilities might be “certainly so,” “probably 
so,”” “possibly so,”’ ‘possibly not so,” “probably not so,”’ 
and ‘“‘certainly not so.’’ Even the most precise and cleat 
cut statement may fall into any one of m possible cases, 
according to the new theory 

Systematic work has been done on the proble m of the 
development of many-valued logics by a number of 
workers, and a survey report of their findings was presented 
by Professor J]. Barkley Rosser of Cornell University at the 
November meeting of the American Mathematical Society 
held at Wayne University, Detroit, Michigan. The prin 
cipal result, according to Professor Rosser’s statements, is 
the discovery that a many-valued logic proves to be much 
more complicated than anyone expected. In a three-valued 


logic, for instance, there are 256 distinct meanings of the 
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simple word, “and”; in a four-valued logic there are 
14,348,907 distinct meanings of the same word. All this 
complication makes things hard for the mathematician, 
but allows more freedom for the physicist. Imagine, pro 
ponents of the theory point out, how many possible kinds 
ol atoms a physicist could make out of electrons and 
protons with 256 kinds of “‘and”’ to choose from, and how 
much more chance he’d have of finding one that worked 
the way it ayght to! 


* 


Sanitizing Service 


‘*Protecto-Ray Bathrooms” is the term applied to a new 
and exclusive service of sanitizing bathrooms in the Hotel 
New Yorker in New York City, according to Frank L. 
\ndrews, the hotel’s president. A special Protecto-Ray 
machine has been developed by the hotel in conjunction 
with the General Electric Lighting Laboratory in Nela 
Park, Cleveland. It is a specially constructed portable unit 
on wheels with ultraviolet ray equipment of high intensity 
which, the announcement states, destroys all bacteria, 
germ life and fungi known as athlete’s foot. These ma- 
chines, fifty of which are in use in the hotel, are operated 
by the regular maids. After each bathroom is cleaned 
upon the departure of a guest, a Protecto-Ray machine is 
wheeled into the bathroom and turned on for a minimum 
of ten minutes, regulated by an automatic timing device. 
The bathroom door is then sealed with Cellophane, 
assuring the next guest that the bathroom is thoroughly 
sanitary. 


* 


Mathematical Biophysics 


The table of contents of the March issue of the Bulletin 
of Mathematical Biophysics, edited by N. Rashevsky and 
published by the University of Chicago Press, includes th 
following articles: 

Some Remarks on the Movement of Chromosomes During Cel! 

Division, by N. RASHEVSkKy; 

On Reinforcement and Interference Between Stimuli, by G. Youn: 

Studies in the Mathematical Biophysics of Discrimination and Con- 
ditioning. I, by H. D. LANDAUL; 

Cellular Forms and Movemerts, by A. S. HouseHOLDE: 


* 


Spectrographic Analysis Conference 


A second Midwest Conference on Spectrographic Analy 
sis will be held in the Kent Chemical Laboratory at the 
University of Chicago on Saturday, March 29, according 
to an announcement received from that Laboratory. Bi 
cause of the growing importance of absorption spectroscopy, 
the announcement states, this year’s conference will be 
held in two divisions, one dealing with emission spectra, 
the other with absorption spectra. The former is being 
organized by W. C. Pierce, of the University of Chicago, 
and the latter by Elmer Miller, of the University of Minne 
sota. Further information and copies of the final program 


may be obtained from either of these men. 
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New Books 


Physical Chemistry of High Polymeric Systems 








H. MArK, Adjunct Professor of Organic Chemistry 
Polytechnic Institute of Brooklyn, New York, New 
York. Volume 2 of J/ligh Polymers. Pp. 345+vii, 
Figs. 99, 16234 cm. Interscience Publishers, Inc 

215 | ourth Ave New York 1940, Price $6 50. 

xcept for minor changes and additions, this book is a 
translation of the author’s recently published Ad/gemein 
Grundlagen der hochpolymere Chemie, the first volume of 
Meyer and Mark's /lochpolymere Chemie. Ein Lehr- und 
Handbuch fiir Chemiker und Biologer. 

Physical and physico-chemical methods of value for 
studying relationships between composition and structurs 
ind between structure and properties are, with few excep 
tions, the same for high polymeric systems as for systems 
containing only small molecules. However, because the 
experimental and theoretical complications are generally 
much greater when high polymers are being studied, a 
combination of different methods must frequently be 
employed—and even then a thoroughly satisfactory inter 
pretation of the experimental data is often impossibk 
It, therefore, becomes advisable for anyone studying the 
physical chemistry of systems containing high polymers to 
have as thorough an understanding as possible of thi 
capabilities and limitations of the various conceivably 
applicable research methods. Hence this book was written 

not primarily to summarize the present state of knowl 
edge concerning the physical chemistry of high polymers 
but to provide a sound basis for future progress in this 
field, as well as for the critical evaluation of past work 

lhe author has done a remarkably fine job in combing 
the literature for significant articles, in critically extracting 
from them their essential points, and in presenting these 
briefly, yet clearly and understandably, in their proper 
perspective relative to each other 

The treatment in most of the book is nonmathematical 
ilthough the fundamental equations related to the more 
important phenomena and methods are given and the 
language of the calculus is used fairly extensively in thi 
treatment of such subjects as osmotic phenomena, solu 
bility, viscosity and polymerization kinetics. Many refer 
ences are included to the more important papers and books 
in the fields covered, thus directing the reader to appro 
priate sources for details of the mathematical and experi 
mental methods, further applications, etc. The choice ot 


material for the numerous tables and figures is excellent 
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The first chapter, ‘“The geometry of molecules,’ con- 
tains sections dealing with x-ray and electron diffraction, 
polarizability and dipole moments, and the Kerr effect. 
lhe second chapter, ‘‘The internal motions of molecules,” 
presents mainly some results of band spectrum and 
specific heat studies. The ‘‘Primary and secondary valence”’ 
chapter is a discussion of intra- and intermolecular forces. 
In the chapter on ‘‘Crystal structure, lattice forces and 
chemical structure of the unit cell,’’ crystals are classified, 
first on the basis of the interatomic forces involved, then 
on the basis of their symmetry. ‘Molecular arrangement in 
liquids and the investigation of minute crystals’’ deals 
with x-ray methods for studying such systems. The chapter 
on “The behavior of mixtures’’ presents a condensed 
thermodynamic treatment of mixtures, with especial 
regard to such properties as osmotic pressure and solu- 
bility, in addition to an admirable discussion of the vis 
cosity of liquids and solutions. In the final chapter 
“Kinetic phenomena in high polymer chemistry,” diffusion 
phenomena and the kinetics of the formation and degrada 
tion of long-« hain molecules are considered. 

In conclusion, the reviewer is glad to recommend this 


volume, not onl 


to those whose chief interest is in the 
high polymer field, but also to those interested in corre- 
lating the properties of /ow molecular weight substances 
with their structures. 


M. LL. HuGGins 


Eastman Kodak Compan 


The Story of Superfinish 


ARTHUR M. SWIGeERT, JR. Pp. 672, Figs. 720, 16 233 
em. Lynn Publishing Company, Detroit, 1940. Pric 
$5.00. 

This book describes a recently developed process [ol 
producing commercially fine finishes on machine surfaces 
ol any shape which hitherto could be approached In 
quality only by expensive lapping processes. The new 
process employs a random abrasive action which geo 
metrically develops a smooth true surface by eliminating 
the surface points and projection above the bearing ol 
contact surtace 

I UGENE PERKINS 


University of Pittsburgh 


Handbook of Physics and Chemistry 


24th Edition Pp. 2564, 125183 cm. Chemical 
Rubber Publishing Company, Cleveland, Ohio, 1940. 
Price $3.50 
I:very vear this important handbook becomes bigger 
and better. The physicist, in particular, will be interested 
in the sections on mathematical constants; physical con 
stants in the fields of heat, sound, electricity, light, and 
specific gravity; and in the tables of elements and thei 
properties In the present edition 7 pages are devoted to 
induced radioactivities. 
\ 3l-page index makes most of this material readily 
ivailable. The handbook is well-bound and costs less than 


one-seventh cent per page. 
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An X-Ray Study of the Changes That Occur in Aluminum During the 
Process of Fatiguing 


3 RAYMOND G. SPENCER AND 


ve 


l/bion Colle 


Received November :. 


Carefully prepared test bars of 17S—T aluminum alloy 


have been fatigued; x-ray diffraction patterns of each 


specimen were taken at intervals as fatiguing progressed 
Approximately thirty different test bars have been studied 


Seven series of diffraction patterns, corresponding to 


seven different specimens, are reproduced and discussed in 


this paper. This study shows that (1) fatigue failure of 


17S-T aluminum is always preceded by fatigue damag: 

which is detectable by the use of x-ray diffraction patterns; 
INTRODUCTION 

LARGE percentage of the very large 


A 


of metals consists merely of 


number of experimental studies of fatiguing 
subjecting the 
samples to repeated cvcles of stress and recording 
} the number of cycles required to cause failure. 
, These data vield the so-called S— N curves which 
indicate a more or less accurate endurance limit 
and give, with rather poor accuracy, the number 
of cycles required to cause failure at a given 
stress, but give no information at all about what 
happens within the sample to cause it to fail or 
what physical changes, if any, precede failure. 
At stresses much above the endurance limit the 
number of cycles required to produce failure of 
different specimens of the same material varies 
over wide limits. For example, one can say that 
different specimens of the same material, fatigued 
in exactly 
will fail 


between two and three million cycles of stress. 


the same manner at the same stress, 


after being subjected to somewhere 


Fatiguing of metals might be expected to 


cause any one, or all, of the following reactions: 


a 


(1) rotation or reorientation of some of the 


crystallites, which may or may not be a pre- 


ferred orientation, (2) a warping of some of the 
crystallites which 
condition 


would represent a_ strained 
lattice 


(3) fragmentation of some of the crystallites and 


with consequent distortion, 
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) 


2) nearly all of the 


detectable changes occur in the last 


stages ol tatiguing and consist of rotation of erains, 


introduction of strains in some of the crystallites, and 


fragmentation of some of the crvstallites; (3) these changes 


do not occur when the specimen is fatigued below its 


endurance limit, but always occur 


to failure when 
limit; (4 
than 


prior 


it is fatigued above the endurance failure is 


probably intercrystalline rather intracrystalline in 


character 


+) a gradual weakening and consequent breaking 


of the 


intercrvstalline bonds. Under suitable 


conditions, any one of the above-mentioned re- 
actions, except (4), will produce a corresponding 
and detectable change in the x-ray diffraction 
pattern. If the nature of the specimen in its 
original conditions is such that the diffraction 
lines consist of a row of separate spots, frag- 
mentation would produce a larger number of 
spots, rotation of a crystallite would change the 
position of the corresponding Spot, and lattices 
distortion of a crystallite would cause the corre 
sponding spot to become blurred. kragmentation 
of sufficient extent would cause the spots to 
coalesce. Lattice distortion in a specimen which 
has crystallites small enough to produce uniform, 
solid diffraction lines would cause the lines to 


widen and decrease in intensity. It should be 
pointed out that, while most of the changes that 
one would expect to take place within the metal 
when it is fatigued cause changes in the diffrac- 
tion pattern, it is not always possible to be sure 
which of the possible reactions has affected the 
pattern. For example, if the specimen has small 
grains of the order of magnitude of 10~¢ em, 
either further fragmentation or lattice distortion 
would cause the diffraction lines to broaden and 
it would be difficult to 


distinguish between 


the two. 
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Since any reasonable understanding of the 


causes of fatigue failure awaits a knowledge of 
the changes within the specimen which precede 
failure, it is not surprising to note that several 
investigators have attempted to study these 
changes by using x-ray diffraction methods. 
Charles S. Barrett' has given a critical review 
of the literature up to 1936, including a rather 
complete bibliography. Gough and Wood?~ have 
published a series of papers on x-ray studies of 
both cyclic and static stresses in ferrous metals 
and conclude, among other things, that fatiguing 
of steel below the endurance limit has no 
appreciable effect on the crystalline structure, 
that fatigue failure is preceded by definite and 
detectable changes in crystalline structure, and 
that failure due to static and cyclic stressing is 
preceded by the same changes in structure. 
Wever and Moller® point out the rather obvious 
fact that there is a marked difference in the 
response of the specimen from place to place on 
its surface due to cyclic stressing. They followed 
the practice of etching a small spot on the 
surface of the specimen and assumed that the 
maximum damage would occur at this spot. 
The specimen was locked in the x-ray camera so 
that this etched spot was always irradiated. 
They conclude that they can differentiate be- 
tween safe and unsafe stresses. Wever, Hemple 
and Médller® conclude that fatigue damage can 
be detected by x-ray methods in the case of cold- 
worked steel. Barrett’? has studied both steels 
and aluminum alloys. He states that, in general, 
it is not possible to distinguish between safe and 
unsafe stresses from x-ray diffraction patterns. 
He also states that the changes in the diffraction 
patterns reveal cold work due to fatigue, not 
necessarily fatigue damage, and that cold work 
may take place at stresses both above and below 
the endurance limit. Cold work, of course, does 

' Charles S. Barrett, Trans. A. S. M. E. 25, 1115-1143 
(1937). 

2H. J. Gough and W. A. Wood, Proc. Roy. Soc. A154, 
510-539 (1936 

‘H. J. Gough and W. A. Wood, Proc. Roy. Soc. A165, 
358-371 (1938 

*H. J. Gough and W. A. Wood, Proc. Inst. Mech. Eng., 
London 141, 175-185 (1939 

°F. Wever and H. Modller, Naturwiss. 25, 449-453 
1937). 

®*F. Wever, M. Hemple and H. Miller, Archiv f. Eisen 
huttenwesen 11, 313-318 (1937-38). 

7 Charles S. Barrett, Metals and Allovs 8, 13—21 (1937 
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not necessarily weaken the specimen. One of the 
authors® has shown that x-ray methods can be 
used to distinguish between safe and unsafe 
stresses in malleable iron and to reveal the 
progressive changes that occur in the metal 
when it is stressed above its endurance limit. 
His studies show that these changes do not occur 
when malleable iron is stressed below its en- 
durance limit. Kies and Quick® have worked 
with 25S aluminum alloy and conclude that 
x-ray diffraction patterns do not reveal any 
significant structural changes due to fatigue 
stressing of the alloy. 


MATERIALS USED AND PREPARATION OI 
THE SPECIMENS 


17S-T aluminum alloy, specially prepared bars 
of which were kindly furnished by Mr. D. W. 
Smith of the Aluminum Company of America, 
were used in all of the tests described in this 
paper. The 17S-T ingot was hot rolled to ; of 
an inch diameter. This rod was then cold- 
worked to 0.625 of an inch diameter, straightened 
and heat treated at 950°F. The samples were 
handled in such manner that no internal strains 
could have been introduced into them prior to 
the time that they were put into the fatiguing 
machine. Mr. Smith states that 17S—T aluminum 
alloy has a tensile strength of 61,400, a vield 
strength of 39,000 and an endurance limit at 
500,000,000 cveles of 15,000 pounds per square 
inch. The endurance limit of the specimens as 
we prepared them is about 2000 pounds per 
square inch less than this. 

I-ach specimen was machined to one-half inch 
diameter. A constriction of circular contour was 
then machined near one end. The smallest 
diameter of this constriction was 2 of an inch. 
The last half-dozen cuts of the lathe tool when 
machining the constriction were about 0.0005 of 
an inch in depth. The specimen was then polished 
with fine aloxite cloth. Finally, the specimen was 
etched to remove any crystallites which might 
have been strained during the process of machin- 
ing and polishing, and to remove any finely 
powdered aluminum and abrasive dust which 

* Raymond G. Spencer, Phys. Rev. 55, 991-994 (1939). 


‘J. A. Kies and G. W. Quick, Report No. 659, National 
Advisory Committee for Aeronautics (1939). 
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might have been imbedded in the surface. Speci- 
mens with etched surfaces, of course, have a 
lower endurance limit than polished surfaces. 
However, crystallites which have been broken 
up or strained by machining and polishing must 
be removed from the surface or the results will 
be confusing. A study of surfaces prepared in 
this manner will reveal what changes occur due 
to fatiguing not confused by strains or frag- 
mentation of crystallites which have been pro- 
duced by other processes before fatiguing has 
started. 


THE CAMERA 


A circular camera was used. The film was on 
the circumference and the part of the specimen 
which was irradiated was in the center of the 
circle. The collimating system was so designed 
that a large number of well-defined. spots appear 
in each line. The x-ray beam struck the specimen 
in the middle of the constriction on its minimum 
circumference at a point which was subjected to 
the maximum stress. However, it must be noted 
that the x-ray beam did not necessarily strike 
at the point of maximum strain or where the 
first fatigue crack would start. 

Iron K-alpha x-rays were used. This soft 
radiation penetrated a comparatively short dis- 
tance into the specimen and gave large deviations 
of the diffracted beams. Both of these conditions 
have advantages for this type of work. Each 
diffraction pattern contains lines from the 111, 
100, 110, and 311 sets of planes and the second- 
order lines from the 111 and 100 sets of planes. 
The angle of deviation of the second-order 
diffracted beam from the 100 planes is roughly 
145 degrees. 

The specimens were fatigued in an ordinary 
rotating beam machine in our laboratory. They 
were removed from the machine at intervals as 
fatiguing progressed so that diffraction patterns 
could be made. The camera was so designed 
that a specimen could be removed and replaced 
in the same position. Thus, the same area was 
irradiated each time an exposure was made. 
If, after the specimen was broken, it was found 
that the fracture passed through this area, it was 
put in the camera in a slightly different position 
so that the x-ray beam struck it a small distance 
away from the edge of the fracture. When 
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diffraction patterns were made after a specimen 
was broken, the x-ray beam was never permitted 
to strike the specimen so. that any part of it 
struck a fatigue crack or the edge of a fracture. 
(Except, of course, when we wished to study the 
structure of the fractured surface. The x-ray 
beam was then directed upon the face of the 
fracture.) Under these conditions, we think it is 
logical to assume that the damage revealed by 
the diffraction patterns taken after failure 
occurred prior to failure, rather than to assume 
that the damage was caused by the breaking of 
the specimen. To test this assumption we took 
diffraction patterns at different places on the 
circumference adjacent to the edge of the frac- 
ture of some of the specimens after they had 
broken. We found that in some positions very 
marked damage was revealed, in others only 
slight damage was indicated, and sometimes a 
position was found where no damage at all was 
revealed. For example, Fig. 4(l-) shows a place 
adjacent to the fracture where very marked 
damage had occurred ; Fig. 4(F) shows another 
position on the circumference adjacent to the 
fracture where only slight damage was revealed. 
If the act of rupturing caused the damage indi 
cated, this damage should have occurred at all 
positions adjacent to the fracture. Whenever we 
made a survey of different positions on the 
circumference of a broken specimen, we always 
found one position where very marked damage 
was indicated, and other positions where little or 
no damage was revealed. We assume that the 
fracture which caused the specimen to fail 
started on the side of the specimen where the 
marked damage was indicated. It seems obvious 
that the damage revealed in the diffraction 
patterns occurred prior to failure and was the 
cause of failure. Since, as will be pointed out 
later, failure follows very soon after detectable 
damage occurs, the probability that a specimen 
could be removed from the fatiguing machine at 
the right time, prior to failure, for this damage to 
be evidenced in the diffraction pattern is very 
small. Nevertheless, a few specimens were re- 
moved at the proper time to have damage re- 
vealed prior to failure. The diffraction patterns 
of specimen No. 7 show clearly that detectable 


damage occurred before failure. 
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Fic. 1. X-ray diffraction patterns from specimen No. 1 
Stress 11,000 Ib./sq. in \) before fatigue. (B) after 1000 
cvcles “ 5,000,000 1) 17,000,000 E) 23,000,000 
Fk) 60,000,000 


RESULTS 


Specimen No. 1 was fatigued at a stress of 
11,000 Ib. sq. in., which is approximately 2000 
lb. sq. in. below the endurance limit. Figure 1(A) 
is the diffraction pattern taken before fatiguing 
started, (B) after 1000 cycles, (C) after 5,000,000 
cycles, (D) after 17,000,000 cycles, (FE) after 
23,000,000 cycles, and (F) after 60,000,000 evcles. 
The specimen did not break. The spots in thi 
lines of these patterns are sharp and clear and 
show no changes with increasing fatigue. Even 
after 60,000,000 cycles the diffraction pattern is 
unchanged. 

Specimen No. 2 was fatigued at a stress of 
16,500 Ib. sq. in. Figure 2(A) is the diffraction 
pattern taken before fatiguing started, (B) after 
1000 cycles, (C) after 1,200,000 evcles, (D) after 
7,500,000 cycles, (E) after the sample had 
broken, and (F) at a different place on the 
circumference after the sample had broken. The 
specimen broke after 8,200,000 cycles. Only very 
slight changes occurred in the diffraction patterns 
until sometime after the specimen had undet 
gone 7,500,000 cycles of stress. All the changes 
indicated in (EF) and (F) took place sometime 
during the last 700,000 cvcles. 

Specimen No. 3 was fatigued at a stress of 
22,700 Ib. sq. in. Figure 3(A) is the diffraction 
pattern taken before fatiguing started, (B) afte 
1000 eyeles, (C) after 1,500,000 evcles, (D) after 
2,200,000 cycles, (I) after 2,500,000 cycles, and 
IF’) after 2,800,000 cycles after the specimen had 
broken. (E) indicates that no appreciable damage 
had occurred up to that time, but (F) indicates 
very marked changes. These changes must have 
taken place sometime during the last 300,000 


evel 5. 
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Specimen No. 4 was fatigued at a stress of 
27,500 Ib./sq. in. Figure 4(A) is the diffraction 
pattern taken before fatiguing started, (B) after 
1000 cycles, (C) after 640,000 cycles, (D) after 
744,000 cycles, (IE) after the specimen had 
broken at 760,300 cycles, and (F) at a different 
place on the circumference after the specimen 
had broken. Here again, nearly all of the fatigue 
damage occurred in the last stages of fatiguing, 


in fact, nearly all the damage indicated by (E) 





Fic. 2. X-ray diffraction patterns from specimen No. 2 
Stress 16,500 Ib./sq. in. (A) before fatigue. (B) after 1000 
evcles. (C) 1.200.000. (1D) 7.500.000. (E) and (F) 8,200,000 


took place during the last 16,000 cycles. (Fk) 
indicates much less damage than (IE), which 
shows that different areas of a specimen sub- 
jected to the same stress may respond very 
differently. This shows that x-ray studies of 
fatiguing are probably misleading, if arrange- 
ments are not made to have the x-ray beam 
survey the same area each time a diffraction 
pattern is made. 

Specimen No. 5 was fatigued at a stress of 
30,000 Ib. sq. in. Figure 5(A) is the diffraction 
pattern taken before fatiguing started, (B) after 
1000 cycles, (C) after 440,000 cycles, (D) after 
the specimen had broken at 560,500 cycles, and 
IX) at a different place on the circumference 
after the specimen had broken. In this case, 
marked changes in crystalline structure evi- 
dently took place during the first few cycles of 
fatiguing. (B) taken after only 1000 cvcles 
indicates some fatigue damage. Fatigue damage 
increased steadily until failure occurred. No other 
specimen that we fatigued responded in_ this 
manner. It was first thought that this behavior 
might be due to the relatively high stress at 
which the specimen was run and that the re- 
sponse was due to plastic flow or cold work of 


the surface as Barrett contends. To check this 
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possibility, several other specimens were run at 
the same and higher stresses, but in all of these 
cases practically all of the fatigue damage 
occurred in the last stages of fatiguing. For 
example, see Fig. 6. 

Specimen No. 6 was fatigued at a stress of 
32,000 Ib., sq. in. Figure 6(A) is the diffraction 
pattern taken before fatiguing started, (B) after 
1000 cycles, (C) after 69,500 cycles, (D) after 
153,000 cycles (FE) after the specimen had broken 
at 181,300 cycles, and (F) at a different place on 
the circumference of maximum stress after the 
specimen had broken. Here again, the detectable 
fatigue damage took place in the last stages of 
fatiguing. 

Specimen No. 7 was fatigued at a stress of 
26,000 Ib. sq. in. Figure 7(A) is the diffraction 
pattern taken before fatiguing started, (B) after 
600,000 cycles, (C) after 700,000 cycles, (D 
after the specimen had broken at 725,000 cycles, 
and (FE) is a diffraction pattern taken with the 
x-ray beam incident upon the broken end. 
(C) indicated that definite, detectable damage 
had occurred prior to failure. Since neither (C 
nor (D) indicates the marked damage that is 


often revealed, it is probable that more damage 





Fic. 3. X-ray diffraction patterns from specimen No. 3. 
Stress 22,700 Ib./sq. in. (A) before fatigue. (B) after 1000 
eveles. (C) 1,500,000 1) 2,200,000. (E) 2,500,000. 
F) 2.800.000. 


would have been revealed at some other position 
on the circumference if a survey of the entire 


circumference had been made. 


DISCUSSION AND CONCLUSIONS 


Although diffraction patterns for only seven 
specimens are shown in this paper, at least 
thirty different specimens were studied. The 
fatigue stresses for the seven specimens dis- 
cussed here vary from 2000 Ib./sq. in. below the 
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endurance limit to stresses of more than two and 
one-half times the endurance limit. In all in- 
stances except one, specimen No. 5, nearly all 
of the detectable fatigue damage occurred in the 
last stages of fatiguing. Since no other specimen 
of the thirty studied responded as did No. 5, we 
feel that its behavior is due to some unknown, 
spurious condition. We think that it is per- 
missible, logical and necessary to neglect it in 
drawing general conclusions. It is included here 
only for the sake of scientific completeness. 
Anyone who reads this paper may attach what- 
ever significance to it that he wishes. 

Distortion of various types within a crystallite 
would cause an increase in the size of its corre 
sponding spot in the diffraction pattern. Increase 
in size of some of the spots due to this cause is 
evident. However, some of the large spots which 
we found, we think, cannot be due to distortion. 
If the increase in spot size is due to distortion, 
the density of the spot should decrease, partly 
due to the fact that the diffracted x-rays are 
distributed over a larger area on the film. Some 
of the spots increase in size and, also, in density. 
This is exactly what one would expect if two 
crystallites should rotate so that their spots 
coalesce. These large, dense spots do not occur 
until the specimen is near failure. 

The suddenness with which failure follows 
detectable damage is of considerable interest. 
Detectable fatigue damage indicates that failure 
is imminent. It does not seem reasonable to 
assume that the cyclic stresses cause no damage 
within the metal until failure is almost reached 
and then the incidence of damage sets in. It 
seems to us highly probable that the specimen 
is damaged in some way in the early stages of 


fatiguing, but that this damage is not revealed 





Fic. 4. X-ray diffraction patterns from specimen No. 4 
Stress 27,500 Ib./sq. in. (A) before fatigue. (B) after 1000 
cycles. (C) 640,000. (D) 744,000. (E) and (F) 760,000. 
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in the diffraction patterns. Since the diffraction 


patterns show very slight or no change until 


failure is imminent, this early damage, if it 
occurs, cannot consist of fragmentation of crys- 
tallites or introduction of strains into them. If 
the failure of the specimen is due to inter- 
crystalline rather than intracrystalline failure, 
we would expect the type of response that is 
observed. Any weakening of the intercrystalline 
bonds would not be evident in the diffraction 
patterns. Nothing would be evident until some 
of these bonds gave way and the crystallites 
started to rotate and to slide over one another. 
This would cause the strains and fragmentation 
which are evidenced in the diffraction patterns 
just prior to failure. 

We feel that our studies warrant the following 
Fatigue failure of 17S—T 


is always preceded by 


conclusions: (1) alu- 


minum alloy fatigue 


methods. 


damage which is detectable by 


2) Nearly 


X-ray 


all of the detectable damage occurs 





Fic. 5. X-ray diffraction patterns from specimen No. 5. 
Stress 30,000 Ib./sq. in. (A) before fatigue. (B 1000 
cycles. (C) 440,700. (D) and (F) 560,500 


atter 


in the last stages of fatiguing and consists of 
rotation of some of the grains, introduction of 
strains into some of the crystallites, and frag- 
mentation of some of the crystallites. (3) These 
changes do not occur when the specimen is 
4) Repeated 


cvclic stresses gradually weaken the intercrystal- 


fatigued below its endurance limit. 
line bonds which cause the metal to fail. 


Conclusion No. 1 does not agree with Kies and 


Quick® who state that no significant changes due 
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Fic. 6. X-ray diffraction patterns from specimen No. 6. 
Stress 32,000 lb./sq. in. (A) before fatigue. (B) after 1000 
cycles. (C) 69,500. (D) 153,000. (E) and (F) 181,300. 





Fic. 7. 
Stress 
600,000 cycles. (C 
end. 


X-ray diffraction patterns for specimen No. 7. 
26,006 Ib./sq. in. (A). before fatigue. (B 
712,000. (D) 735,000. (E 


alter 
from broken 


to fatiguing are apparent from their diffraction 
patterns. However, for every one of our speci- 
mens that failed, diffraction patterns revealed 
the fact that marked fatigue damage at some 
position on the circumference of maximum stress 
had preceded failure. Barrett! does not agree 
with No. 3. 


conclusion However, we 


are con- 
vinced that our experimental evidence justifies 
this conclusion. We do not argue that the 


experimental evidence proves conclusion No. 4. 
This conclusion, however, is consistent with the 
facts and it difficult to 


experimental seems 


explain them without it. 
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Use has been made of the radioactive properties of colloidal thorium dioxide, intravenously 
injected into rabbits, to study the functional capacity of the reticulo-endothelial cells for 
phagocytosis. A mathematical analysis of the problem is presented which indicates that when 
the system is not overloaded, the foreign material should be removed from the blood according 
to the relation 

b 

n = nye 2* + > Moe 

«@ 
Here a>? and n is the number of particles remaining in the blood at any time ¢. Experiments 
on 13 rabbits substantiate this prediction. With a standard injection of 1 ce of colloidal thorium 
dioxide per kilogram of body weight, m decreases according to three distinct exponential 
relations. The system is apparently overloaded initially, but after m decreases to about one 
sixth its initial value (at /=0), the fraction of the amount of material present that is removed 
in unit time increases and the rest of the curve follows the predicted relation. The fact that the 
system is initially overloaded was checked by the use of smaller amounts of thorium. In all 
cases more than 99 percent of the thorium was removed from the blood stream within 6 or 8 
hours after the injection. The technique reported here should serve as a basis for testing many 
theories of phagoc \ tosis now largely based on deduc tion. 


HE extensive use of artificially radioactive the radioactivity is apparently insufficient to 

elements such as radio-phosphorus and produce harmful effects. 
radio-sodium as tracers in studies of the chemical When such a colloidal material is injected 
processes involved in plant and animal metabo- into the blood stream, it is removed by cells 
lism is well known. However, the possibility collectively known as the reticulo-endothelial 
that long-lived naturally radioactive elements system. These cells are most numerous in the 
can be used in studying other physiological — liver, spleen, bone marrow and lymph nodes but 
processes such as phagocytosis apparently has in addition are widely distributed in the con- 
not been fully explored. The disadvantages nective tissues of the body. Obviously the 
usually presented against the use of naturally functional activity of such a widely distributed 
radioactive elements are that they do not enter group of cells is difficult to study. It was ap 


into the normal chemical reactions of the body parent, however, that the ability of these cells 


and that the cumulative effect of their radiations to remove colloidal thorium from the circulating 
are harmful to the tissues over a period of time. blood would serve as a functional test of the 
Colloidal thorium is, however, a very satisfactory — reticulo-endothelial system. Therefore, a meas 
substance for a quantitative study of phago- ured amount of thorium dioxide was injected 
cytosis because it is readily taken up by the into the blood stream of rabbits and the decrease 
phagocytic cells and the relative amount in the in the radioactivity of blood samples taken at 
circulating blood, at any time, can be determined — selected intervals was determined. 

from its radioactive properties. Intravenous It is not difficult to estimate the type of 
injections of colloidal thorium dioxide have been — results that should be expected in an experiment 


used in clinical x-ray studies for the past ten of this sort. Since the total mass of particulate 


years, and in the dosage of 1 cc per kilogram 
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matter injected into the blood stream is relatively 
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small and since the dimensions of a colloidal 


particle are extremely small, even in comparison 
to the dimensions of an individual cell, phago- 
cytosis probably takes place almost immediately 
on the contact of a particle with a phagocytic 
cell. Furthermore, the capacity of a cell for 
further phagocytosis is probably not appreciably 
affected by the ingestion of successive particles 
until the mass of material taken up by an 
individual cell is of the order of a few percent 
of the total cell volume. These are not unreason- 
able assumptions since histological examination 
of portions of the liver show that individual cells 
take up literally millions of the colloidal particles 
even when relatively small doses of colloidal 
material are used. This is based on the fact that 
the colloidal particles of thorium dioxide are too 
small (diameter order of 10~* cm) to be observed 
with an ordinary microscope even at the highest 
magnification but agglomerates of the thorium 
dioxide (diameter order of 10~* em) can be easily 
distinguished in the interior of individual cells a 
few hours after an injection. 

Although we have postulated that phago- 
cytosis takes place as soon as a particle makes 
contact with a phagocytic cell, it is evident that 
only a small percentage of the particles present 
actually disappear from the blood stream in a 
given small interval of time, for the trajectories 
of the individual particles are of the same order 
of magnitude as the diameters of the particles 
themselves and most of their motion is due to 
the motion of the supporting fluid, the blood. 
Thus we can assume that the removal of an 
individual particle is purely a matter of chance 
since only those particles are removed that are 
carried into contact with phagocytic cells by the 
blood stream. It is evident, then, that on this 
basis the fraction of the total number of particles 
present that is removed in a given short interval 
of time is constant. 

It may be necessary to postulate also that 
some of the particles are only temporarily 
removed from the blood stream and that a 
fraction of these temporarily removed particles 
is continually being reintroduced into the blood 
stream. The possibility of temporary removal is 
a valid assumption since it has been demon- 
strated that some of the reticulo-endothelial cells 
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of the liver and spleen, after they have phago- 
cyted particulate material, leave their fixed 
abode on the walls of the sinusoids and enter the 
circulating blood. The majority of such circu- 
lating phagocytes are arrested in the capillaries 
of the lung but doubtless a smaller number 
reach the peripheral blood. Furthermore, because 
of the very small size of the colloidal particles 
and the fact that the blood has a finite viscosity 
and its velocity is very small adjacent to the 
cell boundaries, some of the thorium dioxide 
may be temporarily stored in pericapillary fluid 
spaces and on the surface of the general endo- 
thelium. That thorium dioxide so removed from 
the blood may later re-enter the circulation is 
distinctly probable. 

Let us say, then, that a represents the fraction 
of the amount of thorium dioxide present that is 
disappearing permanently in unit time and that 
8 represents the fraction of the amount present 
that is disappearing temporarily in unit time. 
To take account of the reappearance of these 
latter particles in the blood stream, we can let 
y represent the fraction of the amount tempo- 
rarily removed that is reappearing in the blood 
stream in unit time. Then, if we assume that 
initially mo particles of thorium dioxide are 
injected into the blood and that at any time, f, 
there are n particles still present in the blood 
while m,; are permanently phagocytized and mn, 
are temporarily removed from the circulating 
blood, we can represent the variation of m with 
time by two first-order differential equations. 


dn,/dt=an 
dny dt=Bn—yno. 


Here the first equation simply states that 
since a@ is the fraction of the amount of thorium 
dioxide present that is disappearing permanently 
in unit time, an is the instantaneous rate at 
which particles are permanently phagocytized. 
Similarily the second equation states that the 
instantaneous rate at which particles are col- 
lecting in regions which will subsequently re- 
introduce the thorium into the blood stream is 
given by the difference of the rate at which 
thorium dioxide particles are accumulating in 
these regions, Bn, and the rate at which these 
temporarily removed particles are reappearing 
in the blood stream, ys. As in the case of the 


JOURNAL OF APPLIED PHYSICS 








-- 
en 











first equation, the rate at which particles are 
disappearing into a temporary abode is propor- 
tional to the number of particles available for 
phagocytosis, m, the proportionality constant 
being simply 8. Similarly the rate at which 
thorium dioxide particles are reappearing in the 
blood stream is proportional to the number of 
particles which have been temporarily removed 
m2, the proportionality constant being y. By 
combining the above equations to eliminate n, 
we can get two second-order differential equa- 
tions for the variation of nm; and ne with time. 


d*n, dn, 

+(a+B+y) +ayn;=ayno, 
di” dt 
d*n» dn» 

+(a+~6+y) +ayn.=0. 
dt? dt 


If a, 8, and y are true constants (a reasonable 
assumption under the conditions of the problem) 
the solution of these 
ing value for n. 


equations yields the follow- 


n=No—NnN1—Ne 


or 
n=[A+} |noe—*'—[LA —} |mce—"*, 
where 
a+tpBp-—y7 
°  la+b+4)*—4ey)! 
h='(a+6+y)—}((a+8+y7)2?—4ay)! 
an 


s(a+B+y7)+3((a+8+y)?—4ay)’. 


Now if we assume that y is small in comparison 
to a and 8, we can easily simplify this solution 
to give the following value for n, 


b 
n= nye **+-—-mnye*'_ 
2a 
where 
a=at+ sp 
b=ay (a+B) 
and 


b<Ka. 


Thus, following an injection of a particulate 
substance into the blood stream, we might expect 
the number of particles present in the blood to 
decrease with time according to a rapidly varying 
exponential function initially when the first term 
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is most important and then according to a less 
rapidly varying exponential function when the 
first term has dropped to a small value and the 
second term is the predominant one. 

This analysis, of course, only can apply when 
the number of free particles in the blood is small 
enough so that the number of particles making 
contact with a single phagocytic cell at one time 
is a small number—perhaps one. If the amount 
of colloidal material in the blood is large or if 
the distribution of the colloidal material through- 
out must be 
immediately after an injection, the chance of an 


the blood is non-uniform, as it 


individual particle being removed in a given 
interval of time will be reduced and the rate at 
which particles are disappearing at any time 
will decrease. This condition will persist until 
enough phagocytosis has taken place to reduce 
the number of free particles in the blood to a 
value small enough to allow phagocytosis of all 
particles making contact with phagocytic cells. 
That is, the variation of n with time would only 
follow the above equation after a time interval 
long enough to allow for the decrease of n to 
the required value. Between the times /=0 and 
t=t,, when this condition is reached, n would 
be given very closely by 


n=n,'e~*, 


where ¢ is a smaller number than a in absolute 
value. 

These relations were verified by a series of 
tests on 13 rabbits, as follows: Seven rabbits 
were given a standard intravenous injection of a 
colloidal suspension of 25 thorium 
dioxide in 20-percent dextrin (Thorotrast) in 
the amount of 1 cc per kilogram of body weight. 
Two rabbits received per kilo of the 
25-percent Thorotrast and two received 1 cc 
per kilo of 6 percent Thorotrast for a net dose 
of one-fourth the standard amount. Finally one 
rabbit kilo of 25 percent 
Thorotrast and one received 2 cc per kilo of 6 
percent Thorotrast for a net dose of one-half the 
standard amount. The six percent Thorotrast 


percent 


0.25 cc 


received 0.5 ce per 


was used to determine whether a relative increase 
in the amount of dextrin for a given amount of 
thorium had 
difference in 


any effect on phagocytosis. No 
the results could be determined 


when equivalent amounts of thorium were used. 
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Fic. 1. Ionization chamber and linear pulse amplifier for 
testing alpha-ray activity. 


Following the injections, a series of 2 cc blood 
samples were removed from each animal at 
various intervals ranging from 15 minutes to 10 
hours. The blood samples were centrifuged and 
the red cells were discarded after being washed 
twice with normal saline solution and recentri- 
fuged. Previous tests had shown that the amount 
of radioactive material discarded with the red 
cells was small in comparison to the total amount 
present and that discarding the red cells greatly 
simplified the reduction of the remaining ma- 
terial. The washings were added to the blood 
plasma and the total evaporated to dryness and 
charred. The residue was then transferred to a 
small ignition crucible and ignited over a Meker 
burner to remove the organic matter. After 
allowing the crucibles to cool, they were soaked 
overnight in separate jars each containing 400 cc 
of distilled water to remove the soluble salts 
still present. Finally the samples were again 
dried and then tested for their alpha-ray activity 
by means of an ionization chamber and linear 
pulse amplifier as indicated in the block diagram 
of Fig. 1. (The resistance R actually used was 
simply the leakage resistance from the grid of 
the first Western Electric 259B tube in the 
amplifier.) By following this procedure we were 
able to get initial counting rates of approximately 
130 to 150 alpha-particles per minute with the 
largest dosage, while normal blood samples gave 
a counting rate, including the background in the 
counting circuit, of less than 0.3 count per 
minute. 

Figures 2, 3 and 4 indicate the type of results 
obtained when the experimental counting rates 
were corrected for the counting rate obtained 
with normal blood samples. Figure 2 shows 
clearly the speed with which the circulating 
blood is freed of the particulate thorium dioxide 
after the injection of 1 cc per kilo of Thorotrast 
into a typical rabbit. About 85 percent of the 
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particulate matter is removed in the first four 


hours after the injection and more than 99 
percent is removed in five and one-half hours. 
The solid curve in this figure is obtained from 
the exponential lines of Fig. 3. Although the 
break in the curve is not particularly evident 
from the experimental data plotted in Fig. 2, 
it is distinctly visible when the data is plotted on 
semi-logarithmic paper as it is in Fig. 3. The 
close correspondence between the experimental 
points and the exponential relations is very 
striking. In all cases the deviations of experi- 
mental points from the exponential lines can 
be attributed to experimental error. It is clearly 
evident from Fig. 3 that the number of particles 
remaining in the blood decreases with time 
according to three distinct laws. This has been 
verified in all tests. 

Figure 4, which summarizes all the tests made 
with the standard injection of 1 ce per kilogram 
of body weight, shows that the counting rate 
decreases along one exponential curve until it 
reaches the neighborhood of 20 to 25 counts 
per minute but beyond this point which may take 
anywhere from 190 to 410 minutes, depending 
on the animal, the fraction of the amount 
present that is disappearing from the blood 
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Fic. 2. Graph indicating the speed with which the 
circulating blood is freed of the particulate thorium 
dioxide after the injection of 1 cc per kilo of Thorotrast 
into a typical rabbit. 
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Fic. 3. Graph showing how the number of particles 
remaining in the blood decreases with time 


stream in unit time increases by a factor of 4 or 
5 and the counting rate decreases along a steeper 
exponential curve. We can deduce from this 
that when the counting rate drops to approxi- 
mately 23 counts per minute, the number of 
particles present in the blood stream has de- 
creased to such a value that the number of 
particles striking an individual phagocytic cell 
in any given time interval is no longer more than 
the cell can ingest. Finally, there is a second 
break in the curve when the counting rate is 
somewhat less than one per minute and the 
fraction of the number of particles present that 
are disappearing from the blood in unit time 
decreases again. We can assume that from here 
on the number of thorium dioxide particles 
reintroduced into the blood stream has begun 
to be important in comparison with the number 
of particles still circulating in the blood. That is, 
the slope of the exponential curve as predicted 
by the equations above has been reduced by the 
reintroduction of thorium dioxide into the 
circulating blood. In any case, practically all of 
the thorium is removed in about 8 hours. With 
smaller doses (see Fig. 3) the thorium is removed 
in a correspondingly shorter time, for although 
the lower portions of the curves are all of 
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Fic. 4. Graph summarizing all tests made with 1 cc per kilo. 


essentially the same shape the initial amount of 
thorium introduced is much less and conse- 
quently the time necessary to reach a counting 
rate in the neighborhood of 20 per minute is 
considerably reduced. Because of the differences 
between individual rabbits, even of the same 
weight, the dose of 1 cc per kilo cannot be 
expected to lead to an exact comparison between 
the animals. However, the results show that the 
only apparent difference between the curves 
obtained with various rabbits is in the initial 
chance for removal. The first portion of the 
curves in Fig. 4 varied slightly in slope with 
individual rabbits but the other portions were 
all very nearly alike. Finally we can conclude 
that the experimental results verify, at least 
qualitatively, the assumptions made as to the 
mechanism of phagocytosis. Apparently, when 
the blood stream is not overloaded, the removal 
of particulate matter of colloidal size from the 
circulating blood depends solely on the chance 
that a particle is carried into contact with a 
phagocytic cell. When the blood stream contains 
more than the equivalent of about one-fifth or 
one-sixth of a cubic centimeter of Thorotrast 
per kilogram of body weight, however, the 
system can be considered to be overloaded and 


the fraction of the amount of particulate matter 
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present that is removed in unit time is con- 
siderably reduced. 

Many functions have been assigned to the 
group of cells primarily concerned in_ this 
experiment. It is generally agreed, however, that 
one important function of the reticulo-endo- 
thelial system is the removal of any foreign 
material (including bacteria) from the circu- 
lating blood. In this manner the cells are an 
important factor in the defense mechanism of 
the body. Obviously quantitative measurements 
of the functional capacity of this widely dis- 
tributed group of cells under varying conditions 
are important. Colloidal thorium dioxide should 
prove a good agent for such studies since it is 
specifically removed by the reticulo-endothelial 


cells and since its radioactive properties make it 
possible to determine easily the relative amount 
remaining in the circulating blood. The technique 
and the normal data recorded in this experiment 
should thus serve as a basis for testing in a 
more objective way this function of the reticulo- 
endothelial system. 

We are indebted to Messrs. R. L. Guest, B. 
D. Lohmiller and B. E. Lambrecht who assisted 
in the preparation of the blood samples and in 
collecting the data. We also wish to acknowledge 
the assistance of the Wisconsin Alumni Research 
Foundation which through grants-in-aid made 
it possible to build the linear amplifier and 
obtain the animals and other supplies for this 
research. 
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pH measurements on sludge of particulate carbonaceous materials and ash were made with 
a vacuum tube potentiometer amplifier. Dust, precipitated out of flue gas of industrial fuel 
burning plants through the use of the series impinger, gives PH values of sludge in the range 
2.5 to 5. Sludge of ash of coal used in these tests gives pH values in the range 7 to 11.5, de- 
pending upon the seam and locality in which the coal is mined. The results of this study may 
aid the “fingerprinting’”’ of dust, responsible, in part, for atmospheric pollution in American 


cities. 


I. INTRODUCTION 


URING the past fifty years, many cities of 
the middle west, where bituminous coal is 
the main fuel burnt, have been trying to reduce 
pollution of the atmosphere. Until recently, 
efforts were directed mainly to a reduction of 
visible smoke, the source of which can be readily 
determined. An observer, stationed on the roof 
of a tall building can watch the performance of 
hundreds of stacks, and notify violators of the 
smoke ordinance. 
Recently, however, the attention of municipal 
authorities has been directed to contamination of 


* Paper presented at the Washington, D. C. meeting 
of the American Physical Society, April 25, 1940. 

+t Now with Automatic Electric Company (Division of 
Assoc. Tel. & Tel. Corp.), Chicago, Illinois. 
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the air by fuel burning plants that appear to be 
operating with clean chimneys. A reduction in the 
amount of visible smoke followed the intro- 
duction and acceptance by fuel users of the 
forced draft stoker, but the problem of excessive, 
tar-free particulate carbonaceous materials and 
ash often being discharged into and subsequently 
settling out of the atmosphere now confronts city 
authorities. There is a need for a positive method 
of linking carbonaceous materials and ash (either 
entrained in or settling out of the atmosphere) 
with the source. Simply stated there is a need 
for a method of “‘fingerprinting”’ dust. 

This paper presents work done to date by the 
authors in attempting to find some easily 
measurable property of the chimney dust, which 
is also characteristic of the fuel burnt. The 
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success of Wiegand! in classifying types of gas 
carbon according to pH determined for the sludge 
suggested to the authors that similar measure- 
ments could be readily carried out on atmospheric 
contaminants. Wiegand found, for example, that 
“Peerless” or printing ink colloidal gas carbon, 
which is best milled in drying oils, has a sludge 
pH of 2.6 (—logio [47 ]= 2.6). Deactivated rubber 
carbon sludge has a pH of 10.6, and this form of 
carbon is milled into rubber, which usually forms 
stable sols in alkaline media. 


II. COMPOSITION oF CoAL AsH 


Fieldner and his collaborators? have found 
significant differences in the chemical composition 
and softening point for ash obtained from fuel 
mined in various coal regions of the United 
States. The quantities of MgO, CaO, Na.O and 
KO are variable in amount, and hence the 
alkaline properties of ash sludge ought to be 
more or less pronounced for different samples of 
coal. For an ash sample obtained from coal mined 
in the No. 3 Pocahontas bed in West Virginia, 
having a softening temperature of 2730°F there 
was 2.1 percent CaO, 0.9 percent MgO, 1.0 
percent NasO and 0.4 percent K.O present. For 
an ash sample obtained from coal mined in the 
No. 6 Franklin County, Illinois bed, having a 
softening temperature of 2320°F, there was 10.1 


‘e POTENTIOMETER AMPLIFIER 
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Fic. 1. Vacuum tube potentiometer amplifier 
for pH measurements. 


1'W. B. Wiegand, Ind. Eng. Chem. 29, 953 (1937 
Columtian Colloidal Carbons (New York, 1938), p. 54. 

2A. C. Fieldner, A. E. Hall and A. L. Feild, Fusthility 
of Coal Ash and Determination of the Softening Temperature, 
U. S. Bur. Mines Bull. 129 (1918); W. A. Selvig and 
A. C. Fieldner, Fusthility of Ash from Coals of the United 
States, ibid. 209 (1922). 


VOLUME 12, MARCH, 1941 











TEST NO.9 TEST NO.II 
11 SPREADER STOKER a SPREADER STOKER 
COAL FROM COAL FROM 
WESTERN KENTUCKY STAUNTON, ILL. | 
10 ® RUN NO! 10 
© RUN NO2 > RUN NO.2 $ 
® RUN NO3 : RUN NO.3 
9+ ® RUN NO4 9+ RUN NO4 
@ RUN NO.5 @ RUN NOS 
® RUN NO6 ® RUN NO6 
w gt LAUNDRY Ww gi PAINT MFG. CO 
8 3 
2 # 
” 
3 $7 t 
x a = T 3 
a a 
a 6h 4 a 6L 4 
io at 2 2 2 
2| Sx 8 «| 5 8 
x a 
5+ On S+ 2 
« z< > ra 
Ww a =4 
Z = $3 =x 
4+ & 4- 23 wor 
> } z 
a 
me 
3* 3 








Fic. 2. pH values obtained for tests No. 9 and 11. 


percent CaO, 1.6 percent MgO, 0.7 percent Na,O 
and 0.8 percent CaO present. 


II]. MrTHop FOR MEASURING PH OF SLUDGE 


The circuit diagram of a pH potentiometer 
amplifier of the type used in this work is shown in 
Fig. 1. Dr. P. Klopsteg of the Central Scientific 
Company sponsors this type of instrument, and 
it has a wide variety of extremely useful appli- 
cations. Buffers effective at PH 3, 4, 7, 8 and 10 
were prepared according to the tables of Clark 
and Lubs and were used to zero the instrument 
from time to time at a known pH value. The 
potential developed is balanced by the potenti- 
ometer circuit, which is arranged to take into 
account reversal of polarity. Balance is indicated 
by the vacuum tube amplifier and output meter 
arrangement. The instrument can be adjusted to 
compensate for changes in room temperature. 

The sludge is prepared by boiling 0.28 gram of 
solid material in 10 ce of distilled water in a 
covered 100-cc beaker for fifteen minutes. The 
suspension is permitted to cool and settle, and 
about five cc of supernatent liquid is decanted. 
The sludge remaining is transferred to the five-cc 
beaker into which the glass and calomel electrodes 
are inserted. Before this procedure was fixed 
upon, a sample of Illinois coal ash was disinte- 
grated in distilled water and the pH of the sludge 
determined at intervals during the boiling. The 
pH for sludge cooled to room temperature was 
12.0, 12.1, 12.15 and 12.05 at 10, 100, 130 and 160 
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minutes, respectively, after the start of boiling. 
The material was comminuted during the boiling 
with a stirring rod. 


IV. COLLECTION OF SAMPLES 


The main project, of which this work is an 
outgrowth, involved a quantitative study of dust 
loading in flue gases from various types of in- 
dustrial fuel burning equipment, and results of 
the survey have already been published in 
part.** Dust was precipitated out of the flue gas 
in a series impinger, which appeared to be 
suitable for the work because dust, rather than 
tarry smoke, is the material collected. In five 
hours of operation a sample in excess of 10 grams 
may be collected, sampling being performed at 
flue gas velocity. Four to six runs were usually 
made in each plant. The impinger liquid was 
poured through a 325-mesh screen in some cases 
to separate + and — sizes. The dust was 
collected in Alundum filtering crucibles. The 
amount of soluble material was determined by 
evaporating to dryness a 100-cc portion of the 
filtrate. A sample of coal used as fuel in each test 
was analyzed. When a dust collector was present 
in the plant, samples of dust retained in it during 
each run were collected. An attempt was made 
to obtain samples of dust settled outside of the 
boiler plant for comparison with that collected in 
the source. 


V. pH MEASUREMENTS 


Measurements on pH of sludge of impinger 
dust, ash of impinger dust and ash of the coal 
used were made for each run of a test on flue gas 
dust loading in the survey. Figure 2 gives results 
for tests No. 9 and 11 on a laundry and paint 
manufacturing company boiler plant, respec- 
tively. Western Kentucky coal was used in the 
former test and coal from Staunton, Illinois, in 
the latter test. For run No. 6 of test No. 9, the 
coal was changed, with a consequent reduction in 
dust loading, and there remains some doubt as to 
its source. Impinger dust PH values in the range 
3-4.5 indicate that a variable amount of SO; or 


§“‘Low cost flue dust sampler,”” Power 84, 70 (April, 
1940). 

‘Osborne Monnett, “Air pollution in Chicago,”’ paper 
read before the Milwaukee meeting of the American 
Society of Mechanical Engineers, June, 1940. 
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;. 3. pH values for sludge. 


SO. may be adsorbed on the carbonaceous dust 
from the flue gas or liquid, depending upon the 
amount of sulphur in the coal used, the degree of 
subdivision of the dust or the action of the 
impinger in absorbing acid gases. Sludge of ash of 
impinger dust is nearly neutral because of the 
solvent action of the acid liquid in the impinger 
on the original dust entering it. The difference in 
sludge of ash of coal from western Kentucky and 
from Staunton, Illinois, is probably caused by a 
higher content of alkali and alkaline earth ele- 
ments in Illinois coal. 

All of the data obtained in the pH tests was 
plotted as in Fig. 2, but space does not permit all 
of the diagrams to be presented here. Averages 
were made for all runs of each test, and the data 
are shown in Fig. 3. In addition to pH tests on 
sludge of impinger dust, its ash, prepared without 
fusion, and that of the coal from which it is 
derived, tests were made on coal ash treated with 
0.001N and 0.1N HCl, followed by ignition in 
order to find out why the impinger dust ash 
nearly always gives a lower pH value than the 
sludge of ash of the coal from which it is derived. 
In two cases, tests 18 and 20, samples were 
removed from a dust collector associated with the 
stack of the boiler plant. In the first case, a fresh 
water washer is used and ih the second case a 
multiple baffle Bubar type collector is employed. 
Sludge of coal ash averaged 11.1 and 7.8 and 
sludge of ash of dust from the collector averaged 
10.8 and 8.0, respectively, in these two cases. 
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The ash was unfused in both these tests, and 
hence the agreement is fairly good. The pH of 
sludge of dust from the collector for test No. 20 
is 3.9. This value was found inconvenient for 
plotting in Fig. 3. 

From the results of Fig. 3, it appears to be 
fairly certain that the —325-mesh dust has a 
sludge pH which is a few tenths of a unit lower 
than for the +325-mesh dust, probably because 
of greater surface area. A difference is not marked 
for pH of sludge of ash of this test. Intrinsic coal 
ash consists of mineral matter, much of which is 
present in the cell walls of plants prior to 
carbonization. Burning of the coarse particles, 
without fusion, brings about a finer state of 
division. In Fig. 3 S designates a spreader and C 
a chain grate type stoker. O designates an oil 
burner. /7 designates a combination with one 
boiler having forced draft and the second natural 
draft, both being hand fired. 

Dust from an oil burner obtained in test No. 
15, which isa mixture of lampblack and petroleum 
coke particles, gives a PH of 3.0 for the impinger 
dust sludge and slightly less for the —325-mesh 
fraction, which was separated in some runs. The 
ash was neutral. In test No. 17, where pulverized 
petroleum coke was used as fuel the pH for the 
impinger dust fell between 4 and 5, and that for 
the coke ash near to 11.3 (not shown in Fig. 3). 

Significant differences in PH of sludge of ash of 
coal obtained from various localities are also 
shown in Fig. 3. Illinois coal gives the most 
alkaline sludge PH values, ranging from 9.2 to 
11.5 units. The No. 6 seam is worked in many 
localities in Illinois. For the Maumee, Indiana, 
coal samples, where the No. 5 seam is probably 
worked, the pH for the coal ash sludge was much 
lower, falling near to 7.8. The fact that the ash 
remained unfused in the pulverized coal burning 
plant, probably indicates the absence of appreci- 
able quantities of alkali or alkaline earth ele- 
ments in the ash of the coal used in test No. 20. 
The sludge of ash of western Kentucky coal gave 
a PH value near to 8 and that of eastern Kentucky 
coal near to 7.3. If mixtures of coal from various 
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seams were sampled, the pH of sludge of ash 
would fall between the PH values of the sludge of 
either constituent. 

An attempt was made to correlate PH meas- 
urements on sludge of dust and its ash, collected 
by a Lake View housewife complainant, with 
similar determinations on combustion chamber 
dust and ash obtained from a nearby laundry 
equipped with a forced draft spreader type 
stoker, probably burning coal from northern 
Illinois. Values were 7.4 and 11.2 for the first two 
samples, and 11.4 and 11.7 for the last two 
samples, respectively. The similarity in PH values 
for sludge of the dust ash and of the ash of the 
coal used in the laundry was an aid in establishing 
a possible source of neighborhood atmospheric 
pollution. 


VI. SUMMARY 


The use of the PH potentiometer amplifier in 
studies on atmospheric pollution is discussed and 
PH values are presented for sludge of dust and 
ash discharged in the atmosphere by modern 
industrial plants burning various types of mid- 
western coal as well as oil and_ pulverized 
petroleum coke. It is found that the ash of the 
dust obtained with the series impinger gives a 
sludge pH value which is lower than that for the 
ash of coal from which it is derived. However, a 
close correspondence is found between the pH of 
sludge of ash of dust precipitated in other types 
of dust collectors with the pH of sludge of the ash 
of coal from which the dust is derived. PH values 
for sludge of ash of coal from eight different coal 
mines are presented, and it is shown that such 
values for ash of Illinois coal are appreciably 
higher than for Kentucky coal. The use of pH 
measurements to aid in “fingerprinting’’ dust is 
discussed. 

Assistance in the preparation of this report 
was supplied by Official W.P.A. Project 30051, 
sponsored by the Chicago Department of Smoke 
Inspection and Abatement, and co-sponsored by 
the Lewis Institute. 





Some Useful X-Ray Data 


EUGENE W. PIKE 
Andrews and Perillo, Inc., Long Island City, New York 
Received September 23, 1940) 


The measured mass absorption coefficients of Allen and Wrede for H, Be, C, N, O, F, Na, Al, 
Si, Cl and Br have been fitted to the ‘“‘universal law” 


u/p=art"+c/(1+0.048/)h), 


by least squares, and the values of a, nm, and c tabulated. The energy distribution in the con- 


tinuous x-ray spectrum in the range from 20 to 70 kv, as measured by Kirkpatrick and by Ulrey, 


has been recomputed, and presented graphically in a form which permits easy interpolation 


Che recomputation of Ulrey’s data is described in some detail. 


HE widening industrial use of x-rays in the 

range 0.2A<A<1.0A for the inspection of 
packaged foods and similar articles brings an 
increasing need for more complete and better 
organized tabular data on mass absorption 
coefficients and energy distributions in this range. 
The data presented below were assembled by the 
writer for his own use in problems of this kind, 
and have been extremely useful. 


I. Mass ABSORPTION COEFFICIENTS 


The most extensive recent measurements of 
mass absorption coefficients are those of Allen,' 
and those of Kustner and his school, from whom 
Wrede? gives the most complete data. These two 
sets of data agree very well, and after excluding a 
few points obviously widely in error, they were 
combined. For interpolation purposes, the so- 
called ‘‘universal formula” 


u/p=ar"+c/(1+0.048/2r) (1) 


was used, fitted to the observed data by the 
method of least squares. Table I(A) gives the 
constants so obtained, and Table I(B) some 
interpolated values. The elements F and Si, 
marked by parentheses, have not been measured, 
and the constants were determined by inter- 
polation on large accurate charts. 

The constants so derived differ in several cases 
from those derived by Wrede ;? a careful graphical 
check against the original data proves his values 
to be in.error. 

1S. J. M. Allen, Handbook of Chemistry and Physics, 


22nd edition (Chemical Rubber Co., 1938), p. 1502. 
2W. Wrede, Ann. d. Physik 36, 684 (1939). 
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The single value for H in (CHg), is taken 
from Taylor.’ 


Il. THe ENERGY DISTRIBUTION IN THE CON- 
TINUOUS X-RAY SPECTRUM 


The direct experimental measurements of this 
function are limited to the early and incomplete 
work of Kirkpatrick,‘ Dauvillier,> and Ulrey.® 
The former two agree very well over the range 
where Dauvillier’s inadequate corrections for 
second-order reflection from the crystal do not 
prejudice his figures. Kirkpatrick’s corrections 
are complete and accurate ; his final result is ‘‘the 
intensity of x-radiation inside the bulb” of a 
tungsten target Coolidge tube. The target face 
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> E. G. Taylor, Phys. Rev. 20, 709 (1920). 

‘P. Kirkpatrick, Phys. Rev. 22, 37 (1923). 

> A. Dauvillier, Ann. de physique 13, 49 (1920). 
®D. C. Ulrey, Phys. Rev. 11, 401 (1918). 
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TABLE I(A). 
ELEMENT Be ( N O KF Na 
a 0.38 1.11 1.91 2.93 4.33 8.60 
n 3.11 3.00 3.09 3.09 3.07 3.08 
( .188 195 .20 .20 .20 .20 
Br Br 
ELEMENT Al (Si S Cl 4<0.917 X>0.917 
a 13.9 16.9 25.7 30.6 201 26.6 
n 3.01 3.00 2.99 2.96 2.71 2.80 
‘ .20 21 ae 23 
Paste 1(B) 
H H 
r He) CHa)» ( N 0 I Na \l Si Ss Cl Br 
0.28 0.39 0.195 0.207 0.228 0.258 0.341 0.471 0.55 0.77 0.903 6.35 
32 39 .210 .228 .260 305 A431 619 74 1.06 1.252 9.12 
36 39 .228 255 301 365 .546 813 .97 1.42 1.69 12.9 
40 39 .249 .286 351 439 .690 1.058 1.27 1.87 2.24 16.7 
44 39 .274 325 A412 529 868 1.357 1.63 2.43 2.91 21.45 
A8 A3 .303 371 A84 .637 1.078 1.71 2.06 3.09 3.72 27.4 
ae A3 338 425 570 765 1.32 2.42 2.56 2.87 4.64 34.1 
56 43 378 A488 .671 914 1.62 2.60 3.16 4.77 $72 42.4 
.60 A3 424 561 787 1.086 1.97 3.17 3.84 5.81 6.97 50.2 
64 A3 476 .645 921 1.28 2.36 3.81 4.62 7.00 8.40 59.8 
.68 A3 535 741 1.074 ) 2.81 4.54 5.51 8.35 10.0 70.5 
By 44 .26 .601 849 ee 1.77 3.31 5.35 6.56 9.87 11.8 82.4 
76 A4 .675 .970 1.44 2.05 3.88 6.27 7.62 11.54 13.8 95.4 
80 44 157 1.106 1.66 2.37 4.52 7.29 8.85 13.4 16.0 110 
.90 AS .999 1.407 2.30 3.32 6.41 10.31 iZ.3 19.0 22.6 151 
917 159 
21 
1.00 A5 1.303 2.10 3.12 4.52 8.79 14.1 17.1 26.0 30.8 26.6 


was inclined at 45° to the electron beam, and the 
measurements were made in the direction of 
specular reflection of the beam from the target 
face. 

Ulrey® used almost exactly the same experi- 
mental arrangement, but measured only the 
ionization current as a function of voltage and 
spectrometer setting. His description of the 
method used is so complete, however, that the 
necessary corrections may be computed, using a 
providential overlap with Kirkpatrick's work to 
determine the final parameters. 

The plot of Jy, the intensity, against \ for the 
various voltages used is shown in the lower graph 
of Fig. 2. The regularities of this system of curves 
Here the 
function (A7J)) / (A27S)) max is plotted in percent asa 


are shown more strikingly in Fig. 1. 


function of (Ao/A), where Xo is the short wave- 
length limit at the voltage in question. Con- 
sidering the complexity of the processes involved, 
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it is unlikely that there is any theoretical sig- 
nificance to these regularities, but they do 
facilitate interpolation. If the straight portions of 
the curves are prolonged, their intercepts with 
the axis of (Ao/A) form an almost linear pro- 
gression with voltage, as is shown in the upper 
curve of Fig. 2. 

The values of J, as a function of voltage and 
wave-length are given in Table II. 


Ill. THe REDUCTION OF ULREY's DATA 


Ulrey measured simply the ionization current 
as a function of the spectrometer setting for 
several voltages. There are four major corrections 
to be applied to these data: (A) The efficiency of 
the ionization chamber. (B) The reflectivity of 
the crystal. (C) The absorption in materials such 
as the wall of the x-ray tube. (D) Second-order 
reflections from the crystal. These will be con- 
sidered in order. 
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A. The efficiency of the ionization chamber 


Ulrey’s ionization chamber was 75 cm long and 
7.5 cm in diameter, and contained air at atmos- 
pheric pressure, saturated with the vapor of ethyl 
bromide at 20°C. This amounts to 1.61 mg of Br, 
0.95 mg of N, 0.49 mg of C, 0.24 mg of O, 0.09 mg 
of H, and 0.013 mg of A, per ce of gas. The x-ray 
beam was 1.8 cm high, and passed axially 
through the chamber. 

Of the radiation entering the chamber, one 
part goes through, without being absorbed, while 
another part is converted into fluorescent bro- 
mine K radiation. Of this last, a portion escapes 
to the walls. The remainder of the original 
energy is transferred to photoelectrons which 
ionize the gas. Following Binks,’ we assume that 
the energy expended per ion pair is essentially 
constant in this region, so that the ionization is a 
true measure of the total energy of the photo- 
electrons. 


At any wave-length A, the volume absorption 


coefficient 8 of the gas in the chamber can be 


computed from the known mass absorption 


7W. Binks, Reports on Progress in Physics (The Physical 
Society, London, 1936), Vol. III, p. 387. 
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coefficients in Table I. The portion of the original 
beam passing through the chamber unabsorbed 
is exp (—758). The portion lost as fluorescent 
radiation can be estimated as follows: 


In any element of volume, the energy absorbed 
from the beam in a wave-length interval dd will 
be BJ,dvdx. Of this energy, a fraction (8;/8) will 
be expended in producing K ionization in bro- 
mine. (8, can be computed as the difference 
between the absorption coefficient of bromine, 
and the LZ absorption extrapolated from the 
formula for (A>0.917) in Table I.) When a 
bromine atom is ionized in the K level by a 
quantum of wave-length X, a fraction (A/1.03) of 
the energy of the photon is converted into 
bromine K radiation, the remainder goes into the 
primary photoelectron, or into soft radiation 
which does not reach the walls. Finally, only a 
fraction of this K radiation escapes to the walls. 
A numerical integration gives this fraction as 58 
percent. 

Combining all these, the efficiency of the 
chamber is 


£=(1—exp (—758))(1—0.58A8,./1.038). (2) 


This ranges from 0.48 at 0.28A to 0.74 at 0.48A, 
and back to 0.56 at 0.90A. It rises abruptly to 
0.97 at the bromine K edge at 0.917A, and this 
discontinuity shows in Ulrey’s original data, 


although it is smoothed away in his curves. 


TABLE II. 
\ 70.7 KV 51.4 KV 40 KV 30 KV 25 KV 20 K\ 
0.28 7.7 
32 13.8 1.9 
36 17.1 Se 
AO 41 19.5 11.0 
A4 20.3 13.2 aoe 
A8 21.2 14.5 5.2 
50 40.5 23 
52 21.6 15. 6.9 1.26 
56 15 8.0 3.26 
.60 35.5 21 ) 8.3 4.4 
64 8.3 5.0 0.57 
68 8.0 5.1 1.44 
70 28 17 
12 1.6 5.1 1.96 
76 7.0 5.07 2.14 
80 21 12 6.7 4.8 2.24 
84 4.6 2.2 
88 4.15 2.1 
.90 14 8.5 
1.00 8.7 5.6 
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B. The reflectivity of the crystal 

Ulrey states that his crystal was a cleaved 
surface of calcite, ‘“‘the best of some twenty 
pairs.’” Davis and Stempel*® have measured the 
reflectivity of calcite in the range in question, and 
find that the reflectivity of an 
cleaved surface is practically independent of 


“excellent” 


wave-length. 
C. The absorption in the path 


The absorption in Ulrey’s arrangement is con- 
centrated in the glass wall of the x-ray tube, of 


* B. Davis and D. Stempel, Phys. Rev. 17, 608 (1921). 





unknown thickness. However, by requiring that 
Ulrey’s curve join smoothly with Kirkpatrick's in 
the region where they overlap (50,000 volts, 
0.4<’<0.48A), this absorption can be deter- 
mined very accurately. It is equivalent to 0.21 
g/cm? (about 1 mm) of Pyrex glass, a very 
reasonable value. 


D. Second-order reflections 


This correction is large and cannot be com- 
puted. For this reason, the usefulness of Ulrey’s 
curves is limited to the region Ay<A<2Xp, 
where second-order reflections do not occur. 





Cathode Spot Initiation on a Mercury Pool by Means of an External Grid 


Mark A. TOWNSEND 
Lamp Development Laboratory, General Electric Company, Nela Park, Cleveland, Ohio 


(Received January 2, 1941) 


Experiments relating to mercury-pool stroboscope tubes with external starting bands are 


described. Of particular interest is the unidirectional current that is found to flow in the plate 


circuit when grounded and its relationship to the formation of a cathode spot on the mercury 


cathode. Cathode-ray oscillograms are exhibited showing the periodic nature of this current. 


The potential of a disconnected electrode is negative with respect to the cathode and has a value 


of several thousand volts. Curves of the volt-ampere characteristics of an anode are given as a 


function of tube dimensions, gas pressure and other factors. 


INTRODUCTION 


HE mercury-pool cathode, although used 
for many years as a practical and reliable 
source of electrons, is still only partially under- 
stood. Always associated with the use of a 
mercury-pool cathode is the problem of initiating 
the cathode spot. Well-known methods of ac- 
complishing this are: (1) separating mercury 
contacts through which a current is flowing, (2) 
applying a surge of current to a high resistance 
rod immersed in the mercury surface (the 
Ignitron), or, (3) applying a surge of voltage to a 
conductor near the mercury surface but elec- 
trically insulated from it.'~* 
This third method listed above is found to be 
satisfactory and convenient for initiating the 


'K. J. Germeshausen, Phys. Rev. 55, 228 (1939). 
2 Y. Watanabe, H. Kasahara and Y. Nakamura, Elektro- 
tech. J. 2, 180 (1938). 
>L. J. Buttolph and D. W. Dana, Rev. Sci. Inst. 4, 206 
1933). 
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cathode spot in the mercury-pool stroboscope 
tube. The stroboscope tube is usually in the form 
of a glass cylinder from }” to 2” in diameter and 
from 4” to 12” in length. A pool of mercury at the 
lower end serves as the cathode and a nickel or 
iron plate at the upper end serves as the anode. A 
conductor, such as a wire screen, is placed around 
the tube in contact with the outside of the glass 


a) 


ae a 


Fic. 1. Experimental types of mercury-pool 
stroboscope tubes. 
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Fic. 2. Control circuit for a mercury-pool stroboscope tube. 


opposite the mercury meniscus, forming an 
external grid. Figure 1(a) shows a typical form of 
this tube. 

It is possible to start a cathode spot capable of 
carrying a high current discharge through the 
main body of the tube by supplying in the proper 
manner a small amount of energy to the external 
grid. The starting time and duration of the dis- 
charge may be controlled to within a _ few 
microseconds. The light from this discharge may 
be used as an accurately timed source of high 
intensity light for both stroboscopic work and 
high speed motion pictures.‘ In addition, the 
accurate control provided by the stroboscope 
tube makes it useful as a timing element in other 
circuits.® ® 

This paper describes experiments on strobo- 
scope tubes. Under certain conditions a unidi- 
rectional current flows in the anode circuit when 
there is no e.m.f. source therein. A correlation is 
established between this current and the starting 
of the tube. The conditions both of the strobo- 
scope tube and of the circuit, under which this 
unidirectional current is produced, are studied in 
detail. 


Normal tube operation 


A control circuit commonly used with the 
stroboscope tube is shown in Fig. 2. Condenser 
C, is charged to about 350 volts d.c. but is pre- 
vented from discharging by the thyratron tube 7. 
Condenser C; is charged to about 1000 volts d.c. 


*H. E. Edgerton, Elec. Eng. 54, 149 (1935). 

5H. E. Edgerton, K. J. Germeshausen and H. E. Greer, 
J. App. Phys. 8, 2 (1937). 

®*T. S. Gray and W. B. Nottingham, Rev. Sci. Inst. 8, 
65 (1937). 
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which is below the steady-state breakdown 
voltage of the mercury-pool tube. A_ small 
tripping voltage applied to terminals 1—2 allows 
the thyratron to become conducting and dis- 
charges condenser C, through the primary of the 
induction coil. This causes a high voltage oscil- 
lation in the secondary of the induction coil 
which is applied to the external grid of the 
stroboscope tube. The high voltage grid excita- 
tion causes the formation of a cathode spot on the 
mercury pool and the discharge of condenser C; 
through the stroboscope tube. It is this discharge 
of Cs; which produces the stroboscopic light. 


The phenomenon of unidirectional current 


In the circuit of Fig. 2, if the anode supply 
voltage is disconnected by opening switch S; and 
if the grid excitation is continuously applied, an 
unexpected result is obtained: The condenser C; 
acquires a powerful charge. If this condenser is 
shunted by a microammeter, a unidirectional 
current of several microamperes is observed in 
the anode circuit. The direction of flow is that 
which would be produced. by electrons striking 
the anode. A faint glow appears along the length 
of the tube accompanied by streamers of light 
extending from the mercury meniscus a short 
distance along the surface of the glass. The 
mercury surface is greatly agitated. The same 
effects are produced if a 60-cycle sine wave of the 
order of magnitude of 5 kv is substituted for the 
induction coil grid excitation. 

A somewhat clearer picture of the nature of 
this current is given by the oscillograph. In the 
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Fic. 3. Circuit for obtaining oscillograms of the 
unidirectional current. 
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circuit of Fig. 3 the voltage drop of the current 
flowing through resistor R, is impressed directly 
on the vertical deflecting plates of the oscillo- 
graph. The wave shape of the grid voltage is 
superimposed on this by means of an electrostatic 
pick-up obtained by simply bringing the lead 
wire close to the high voltage grid. Figure 4(a) 
shows the form of oscillogram obtained with a 
60-cycle grid voltage. Figure 4(b) shows the 
same thing except that the oscillatory induction 
coil grid voltage is used. Both figures show that 
the unidirectional current consists of a series of 
short, high current discharges occurring at defi- 
nite portions of the grid voltage wave. These 
figures represent many superimposed traces on 
the oscillograph screen. The current varies from 
cycle to cycle and, although the pattern has the 
general shape shown here, the details continually 
change. Notable features of the oscillograms are 
that the current flow is in the same direction 
during both half-cycles of the grid voltage wave 
and that the larger currents flow during the 
negative half-cycle. Measurements with a cali- 
brated oscillograph indicate that the average 
height of the current impulses corresponds to 0.06 
umpere for a typical condition. 


Correlation of the unidirectional current with 
cathode-spot formation 


The multi-anode tube of Fig. 1(c) was con- 
nected in such a manner that the unidirectional 
current due to grid excitation could be studied on 
one anode (No. 3) while an are could be estab- 
lished to a second symmetrically placed anode 
(No. 6). The circuit is shown in Fig. 5. The 
following observations were made: (1) It was 
necessary to apply a grid voltage of at least 6.8 kv 
before an appreciable current could be observed 
in the circuit of anode 3, with no voltage applied 
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Fic. 4. Oscillograms of unidirectional current. (A) with 
60-cycle sine-wave excitation. (B) with oscillatory induc- 
tion coil excitation. 
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Fic. 5. Circuit for studying the correlation of cathode spot 
formation with unidirectional current production. 








to anode 6. With 6.8 kv grid excitation an arc 
would form to anode 6 when the line switch was 
closed. (2) When the grid voltage was reduced 
slightly below 6.8 kv, no current could be ob- 
served in the circuit of anode 3. No arc would 
form to anode 6 when voltage was applied to it. 
(3) With an are to anode 6 and with a grid 
excitation of 6.8 kv continuously applied, current 
did not flow in the circuit of anode 3. If the are to 
anode 6 was suddenly interrupted, a current 
began to flow to anode 3 after a time delay of a 
few seconds. If the arc to anode 6 was interrupted 
and the circuit immediately reclosed, the same 
time delay of a few seconds occurred before an 
are would restrike. 

These observations indicate that the are will 
not strike unless conditions are so adjusted that 
application of the grid voltage produces a 
unidirectional current to a non-energized anode. 

The mechanism involved in the formation of a 
cathode spot in this tube is probably similar to 
that of the Cooper Hewitt lamp described by 
Kenty:’ “The electrostatic field between the 
external starting band and the mercury meniscus 
causes field emission of electrons. The electrons 
are projected up the tube—causing ionization of 
the vapor in the tube. Asa result of this ionization 
a_ self-maintaining discharge becomes estab- 
lished.”” In the Cooper Hewitt lamp the anode is 
connected to the starting band and, during 
starting, the cathode is forced to a high negative 
potential relative to the anode and starting band. 
Kenty’ observed that, even with the anode 
insulated, electrons were sometimes projected up 
the tube with high energies. In the case of the 
stroboscope tube described here the anode is only 
a few hundred volts positive with respect to the 
cathode during starting. Nevertheless, experi- 


7C. Kenty, J. App. Phys. 9, 705 (1938). 
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Fic. 6. Unidirectional current as a function of grid voltage 
for different values of retarding voltage. 


ments show that electrons reach the anode with 
energies corresponding to several thousand volts. 
The remaining portion of this paper describes 
experiments which were performed to determine 
the factors influencing the flow of the unidirec- 
tional current to a non-energized electrode. 


Effect of circuit constants on the unidirectional 
current 


The curves of Fig. 6 show the unidirectional 
current as a function of grid excitation for several 
fixed values of retarding voltage. The current 
was measured by means of the circuit of Fig. 7. 
The microammeter was protected from any large 
alternating current by the shunt condenser and, 
with switch S in position 1, the unidirectional 
current was unopposed. With switch S in position 
2, the d.c. power supply was connected in 
opposition to the flow of current. The magnitude 
of this d.c. retarding voltage was controlled by 
varying the voltage input to the rectifier circuit. 
The 60-cycle grid voltage was controlled by 
means of a variable transformer connected to the 
primary of the high voltage grid transformer. 

We see from Fig. 6 that as the grid voltage is 
increased from zero, no current flows until a 
critical value is reached. Past this point the 
current increases rapidly with the grid voltage. 
The current decreases with an increase in re- 
tarding voltage. This decrease would be expected 
to occur since all the electrons going toward the 
anode will be repelled if their energy is below 
that consumed by the field between the cathode 
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and‘anode. The curves of Fig. 8 show this vari- 
ation of unidirectional current with retarding 
voltage for several fixed values of grid voltage. 
These curves represent the energy distribution of 
the electrons producing the current. 

If the anode is left floating when the grid 
excitation is applied it will acquire a voltage 
corresponding to the maximum energy of the 
electrons striking it. This open circuit anode 
voltage is shown in Fig. 9 as a function of the 
peak value of grid voltage. It is seen that the 
anode voltage is always below the peak value of 
the grid voltage although in some cases this 
difference is quite small. In making these 
measurements of anode voltage a null type of 
electrostatic voltmeter was used. 

Of a similar nature to the voltage acquired by 
the anode is the charge acquired by the tube 
walls. The presence of the wall charges can be 
demonstrated by wrapping a conductor around 
the outside of the tube and connecting it to the 
ground through a galvanometer. When the grid 
voltage is applied, a current flows in the galva- 
nometer circuit until the walls become charged. 


Effect of gas pressure on the unidirectional 
current 


At normal temperatures a mercury-pool cath- 
ode tube always contains mercury vapor at a few 
microns pressure. In a sealed-off tube there is also 
usually present a trace of other gases which come 
from the tube walls and internal parts. In making 
a study of the effects of gases on the flow of the 
unidirectional current it would be desirable to use 


various pure gases for separate experiments. 
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Fic. 7. Circuit for measuring unidirectional current as a 
function of grid voltage and retarding voltage: 
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Fic. 8. Unidirectional current as a function of retarding 
voltage for various values of grid voltage. 


However, in approximating conditions met in 
practical tubes it is convenient to use air. The 
method used was to connect the tube under test 
to a vacuum system and admit air at various 
pressures, maintaining other conditions as nearly 
constant as possible. 

Some typical curves are given in Fig. 10 show- 
ing the current at a fixed value of grid voltage as 
a function of gas pressure. The circuit of Fig. 7 
was used with switch S in position 1 and witha 
60-cycle sine wave of grid voltage. Although the 
tubes used to obtain the data of Fig. 10 differed 
widely in construction, the same general trend is 
observable in all cases. 

In addition to causing the greatly increased 
currents shown in curves 1, 2 and 3 of Fig. 10, the 
presence of the gas affected the value of retarding 
voltage required to stop the flow of current. 
Curve 4 of Fig. 10 shows this “‘stopping voltage”’ 
as a function of gas pressure. As the gas pressure 
is increased the stopping voltage decreases until 
at high pressures the current in some experi- 
mental tubes reverses even with no retarding 
voltage applied. In this range of pressure large 
currents flow with comparatively few volts ap- 
plied to the anode and complicated effects are 
encountered which are not directly connected 
with the phenomenon of unidirectional current 
with which this paper is concerned. It is signifi- 
cant, however, that a faint glow was always 
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Fic. 9. Open circuit anode voltage as a function of grid 
voltage for several different stroboscope tubes. 


observable along the entire length of the tube in 
all experiments where the current was produced. 

In order to test the effect of a very low value of 
mercury vapor pressure, experiments were made 
using dry ice around the tube. It was found that 
current could be produced in a standard type of 
sealed-off stroboscope tube even when the mer- 
cury cathode was solidified by placing it in dry 
ice. However, when the tube was connected to a 
vacuum system, well evacuated and entirely 
surrounded by dry ice, no current was produced. 
In this latter case the permissible grid voltage 
was limited by the fact that moisture collected 
on the outside of the tube. When air was admitted 
to the tube surrounded by dry ice the unidi- 
rectional current was again observed. This cur- 
rent increased with the gas pressure as shown by 
curve 1 of Fig. 11. The fact that currents could be 
produced with the mercury cathode solidified led 
to experiments with other types of metallic 
cathodes. Curve 2 of Fig. 11 was taken with a 
tube in which the mercury cathode was replaced 
by a nickel wire screen similar to the external 
grid except mounted on the inside of the tube. 
These curves are not strictly comparable since 
curve 1 is for a grid voltage of 6.0 kv while 
curve 2 is for a grid voltage of 8.6 kv. In spite of 
this difference there appears to be considerable 
similarity between the results with a mercury 
cathode and with a nickel wire screen cathode, 
especially at the higher gas pressures. 
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Fic. 10. Unidirectional current and stopping voltage as a 
function of gas pressure. Curves 1 and 2 show current for 
anodes 4 and 5 of the tube of Fig. 1(c). Curve 3 shows 
current for the tube of Fig. 1(a). Curve 4 shows stopping 
voltage for the tube of Fig. 1(a). The mercury vapor 
pressure was approximately 1 micron, corresponding to 
rpom temperature 


Effect of tube geometry on the unidirectional 
current 


The anode area was found to have negligible 
effect on the magnitude of the current when other 
conditions were kept constant. One experimental 
method for testing the effect of anode area was to 
vary the amount of mercury in the anode end of a 
“U"’-shaped tube (Fig. 1(b)). Changing the area 
from 0.02 cm? to 3.0 cm? made no appreciable 
difference in the amount of direct current ob- 
tained. Because this U-shaped tube was relatively 
long and narrow, experiments were made with 
the tube shown in Fig. 1(c). This tube was 
designed to reduce the effect of highly charged 
walls in the neighborhood of the anode and to 
provide a control of anode area by having several 
anodes of widely differing area arranged at a 
constant distance from the cathode. Using this 


TABLe I. 
(CURRENT IN ANODE AREA 
MICROAMPERES cM? ANopE No, 
1.2 0.2 5 
1.2 1.26 4 
1.2 5.1 3 
1.4 5.1 6 
1.6 5.3 6+5 
1.6 6.76 6+4 
1.7 10.2 6+3 
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Fic. 11. Unidirectional current as a function of gas 
pressure. Curve 1 for a frozen mercury-pool cathode. 
Curves 2 and 3 for a nickel screen cathode; curve 2 with 
the current unopposed, curve 3 with a 17-megohm resistor 
in series with the anode. Except for curve 1, the mercury 
vapor pressure was approximately 1 micron. 


tube with good vacuum conditions, the current 
was practically independent of the anode area as 
is shown by Table I. 

It is noticed that the current obtained with 
two anodes in parallel is increased over that 
obtained for one alone by about the same 
amount regardless of the total areas involved. 
The current differences shown in Table I may be 
attributed to field effects produced by open 
circuited anodes near those to which the current 
is being measured and are slight compared to the 
range of anode areas covered. 

The cathode-anode distance is apparently not 
an important factor in determining the flow of 
unidirectional current. Tests made with several 
different tubes showed only very slight changes 
in current for changes in anode distance of as 
much as 25 to 1. Table II shows some data 
obtained with the tube of Fig. 1(b) using different 
anode distances and keeping other conditions 
constant. 

Tests made with varying thickness of the tube 
wall at the mercury meniscus indicate that the 
current decreases as the glass thickness is in- 
creased. Table III shows some data obtained with 
tubes identical except for the different wall 
thicknesses at the mercury meniscus. 
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In any tube containing a liquid mercury pool, 
the production of direct current is always ac- 
companied by a strong agitation of the mercury 
meniscus. Tests were made to determine the 
effect of different lengths of excited meniscus by 
varying the length of the external grid so as to 
surround only a part of the tube circumference. 
It was found that decreasing the length of the 
grid by as much as one-half did not appreciably 
affect the current flowing. Cutting down the 
length to one-fourth the circumference reduced 
the current by about one-half. It is necessary 
that some portion of the external grid be opposite 
the mercury meniscus. With the grid entirely 
above or entirely below the mercury edge no 


direct current is produced. The maximum current 


TABLE III. 
Tuse No 1 2 3 
WALL THICKNESS 
INCHES 0.01 0.046 0.092 


GRID VOLTS, KV CURRENT, MICROAMPERES 


4.8 0.2 

6.0 0.5 0.01 

7.8 rupture 0.2 

9.6 of wall 0.4 0.1 
10.8 0.6 0? 


is produced when the center of the external grid 
is placed opposite the mercury meniscus. 

This work was done at the Massachusetts 
Institute of Technology in 1937 in partial ful- 
fillment of the requirements for a Master's 
degree. Earlier work on the starting of strobo- 
scope tubes was done in 1936 by Mr. Charles 
Mueller. Work subsequent to that reported here 
was done in 1939 by Mr. Frederick Gemmill and 
Mr. Louis Hutchins. All these theses are on file in 
the library at the Massachusetts Institute of 
Technology. The author wishes to acknowledge 
gratefully the help and encouragement of Pro- 
fessor H. E. Edgerton who suggested the thesis 
and made many valuable contributions. 





The Effect of X-Rays on the Breakdown Strength and Flashover Voltage 
of Certain Dielectrics 


Eric A. WALKER* 
Hooper Laboratories, Tufts College, Massachusetts 


(Received April 17, 1940) 


Qualitative measurements are made on the reduction of breakdown strength of certain solid, 


liquid, and gaseous insulations. 


he observed reduction in breakdown strength is influenced by 


the gap length as well as the quantity of the x-rays. Exposing a solid to air boundary to x-rays 
increases the flashover voltage. A theory for this effect is offered. 


OBJECT 


T is a well-known fact that x-rays will reduce 
the breakdown strength of insulation although 
little or no data have been presented to show the 
magnitude of this reduction.' It is believed that 
electrical breakdown of insulation is caused by 
three different mechanisms. First, there is the 
: Now at the University of Connecticut, Storrs, Con- 
necticut. 


1F, W. Peek, High Voltage Engineering (McGraw-Hill, 
New York, 1929), 
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well-known “ionization by collision’ phenomena 
as explained by Townsend,’ and extended by 
others.* The second method of electrical break 
down is the pyro-electric phenomenon as offered 
by Wagner and Steinmetz.' This theory states 
that local heating reduces the resistance of a 

2 1.S. Townsend, Electricity in Gases (Oxford University 
Press, 1914). : 

} J. Slepian and R. C. Mason, “Electric Discharges in 
Gases,”’ Elec. Eng. 53, 511 (1934). 


4K. W. Wagner, ‘“‘Physical nature of electrical break 
down in solid dielectrics,’’ Trans. A.I.E.E. 41, 288 (1922). 
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small region in a heterogeneous solid which in 
turn increases the heat generated. Thus a cycle is 
originated which progresses until the resistance 
approaches zero and failure is complete. The 
third cause of insulation failure is found in 
chemical changes which are caused by weathering, 
oxidation or ionic bombardment. The times re- 
quired for these three effects to produce break- 
down differ widely. The collision mechanism 
produces electrical failure in a small fraction of a 
second, the thermal mechanism takes several 
minutes because of the thermal capacity of the 
insulations. The last cause of failure, however, 
takes hours and sometimes years to produce a 
noticeable result. Furthermore, if an agent which 
tends to lower the breakdown strength by 
ionization is used, the insulation will rapidly 
recover its original properties as soon as the 
agent is removed because the ions will rapidly 
recombine. If local heating is the cause of break- 
down it, too, will allow the recovery of the original 
dielectric strength unless the insulation has been 
chemically changed. On the other hand, any 
chemical change has caused an irreparable im- 
pairment of dielectric strength. 

It is believed that the illumination by x-rays 
contributes to lowering the breakdown strength 
only because of the copious supply of free ions 
which are produced in the insulation under 
electrical stress. Certainly no permanent lowering 














Fic. 1, Diagramfof circuit used for tests. 
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ic. 2. Breakdown of air in the presence of x-rays. 


of the breakdown strength was produced, for the 
reduction was only noticeable if the x-rays and 
breakdown voltage were applied simultaneously. 
The heat generated in the insulation by the 
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‘1G. 3. Percent difference of breakdown voltage in 
presence of x-rays. 
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Fic. 4. Breakdown of oil. 


x-rays is negligible. Therefore, it is reasonable to 
suppose that the effects to be described are 
caused exclusively by ionization phenomena. 


METHOD 


Figure 1 shows the arrangement used for most 
tests. The transformer voltage was determined 
by reading the peak voltage on the tertiary 
winding. Such procedure is preferable to reading 
the r.m.s. voltage because it eliminates the error 
caused by a nonsinusoidal wave shape. 

The gaps S; and S, may be changed for differ- 
ent tests, and for tests on oil they may be 
immersed in a suitable tank of oil. 

The x-ray tube was placed about 3 feet from 
the sphere gap and the distance was kept con- 
stant throughout the investigation. The x-ray 
tube was the conventional Coolidge type for 30 
ma current. The x-ray transformer will supply 
80,000 volts peak and the testing transformer is 
rated at 200 kv r.m.s. Usually five readings of 
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1G. 5. Breakdown of solid dielectric (yellow 
varnished cambric). 


each point were taken both with and without the 
application of x-rays. 


RESULTS: BREAKDOWN TESTS 


Figure 2 shows the breakdown voltage between 
two 6.25-centimeter spheres as a function of the 
separation of the gap with and without x-rays. It 
will be noticed that the reduction in breakdown 
voltage is relatively constant. This means the 
percentage reduction becomes less marked as the 
separation of the gaps is made larger. Such a 
result is probably to be expected because, as the 
gap separation becomes larger, there is a greater 
quantity of air and, therefore, more free ions in 
the electric field which are available to initiate 
further ionization. Curve A, Fig. 3, shows the 
percentage decrease in breakdown strength of air 
for this gap as a function of the separation 
between the spheres. 

Figure 4 shows the effect of x-rays on the 
breakdown strength of insulating oil between two 


6.25-centimeter spheres. In this case the differ- 
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ence between the voltage necessary to cause 
electrical breakdown in the presence of x-rays 
and the voltage when they are not present 
increases as the gap is lengthened. In fact, the 
percentage change also increases as is shown by 
curve B of Fig. 3. 

Figure 5 shows the results obtained on yellow 
varnished cambric, a material used extensively 
for transformer insulation. Again there is a de- 
crease in breakdown strength but the decrease is 
approximately proportional to the breakdown 
voltage so that the percentage decrease is ap- 
shown in Fig. 3, 


proximately constant as is 


curve C. 


FLASHOVER TESTS 


The electrode arrangement for flashover tests 
is shown in Fig. 6. Various diameters d and height 
h were used. The results differed somewhat but 
the results shown in Table I for a glass tube 1} 
inches high with a diameter of 1} inches is 
typical. 

This shows a percentage increase of 7.5 percent 
in flashover voltage when the surface is illumi- 
nated with x-rays. The explanation of this 
phenomenon is not difficult. Consider a point P 


as shown on the surface of the glass cylinder in 


TABLE I. 


Wirnour X-Ray With X-Ray 
36.5 38.7 
30.8 38.4 
36.6 40.1 
36.4 40.5 
36.4 40.3 
\verage 36.5 39.6 
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Fig. 6. This point is a point on an impedance 


the the lower 
electrode. This impedance is composed of con- 


which extends from upper to 
ductance caused by dust and moisture on the 
surface of the glass and also capacitance coupling 
such as between the point P and P» just above, 
or P, just below. Now if the conductance and 
capacitance are the only impedances and if they 
are evenly distributed, it is evident that the 
voltage gradient along the surface will be uni- 
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Fic. 7. Voltage distribution on insulator with stray 
capacitance to both electrodes. 


form. However, they are not the only impedances. 
Point P has capacitance coupling through the air 
and glass to both the upper and lower electrodes 
as is indicated in Fig. 6. It has been shown that 
the voltage curve for such an arrangement is non- 
uniform.® A typical distribution curve is shown in 
Fig. 7. 

The gradient is obviously not uniform. How- 
ever, it can also be shown that any decrease in 
the resistance along the surface will cause the 
curve to approach more nearly the ideal voltage 
If the flashover 
gradient is fixed this effect will increase flashover 


distribution shown in Fig. 7. 


voltage. Such a mechanism is thought to be the 
reason an increase was noted in the table. 
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Electromagnetic Waves in Transformer Coils Treated by Maxwell’s Equations 
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(Received November 4, 1940) 


A rigorous solution of Maxwell's equations is derived for waves in a thin single-layer coil, a 


concentric ground electrode being at a short distance from the coil. Simple transformer coils 


resemble such an ideal to a considerable degree. Any electromagnetic state of the coil propa 
gates slowly in axial direction with such a velocity that the phase velocity along the helical 


conductor equals the velocity of light. The surge impedance of the coil is constant for low fre- 
quency, but decreases with higher frequency, the boundary wave-length being determined by 


the insulating distances. 


1. INTRODUCTION 


HE behavior of electric waves and oscilla- 
tions in tightly wound coils is a problem 
which has been discussed for a long time in the 
field of communication circuits! as well as in the 
field of power transmission.*~* The investigators 
agree on the character and on the quantitative 
effect caused by the internal capacitance between 
adjacent elements. There are various opinions, 
however, as to the complete influence of neigh- 
boring another. The mutual 
capacitive and inductive interaction of each and 
every turn in a compact coil is still an open 
question and merits further consideration. 
Particularly for 


turns on one 


the 
various conceptions concerning this interaction 
have led to very different conclusions as to the 
frequency of the natural oscillations of the 
winding and 


transformer windings, 


the traveling velocity of surge 
waves. Some authors believe that the core flux 
in the steel sheets common to all turns plays a 
major part.®® Others think that the transverse 
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lagiger Spulen,”’ Ann. d. Physik 43, 749 (1914). 
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Sprungwellen,’’ Elektrotech.. u. Maschinenbau 32, 729 
(1914); Elektrische Schaltvorgénge (Berlin, 1923/33), p. 
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3K. W. Wagner, ‘‘Das Eindringen einer elektromag- 
netischen Welle in eine Spule mit Windungskapazitat,” 
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suchung der Einwirkung von Wanderwellen-Schwingungen 
auf Transformatoren-Wicklungen,”’ Archiv f. Elektro- 
technik 5, 383 (1917). 
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flux between the elements of the winding is of 
primary importance.” * Still others feel that the 
main agent is the axial flux between the primary 
and secondary coils.® !° 

With situation, it may be 
interesting to show that a rigorous solution of 


regard to this 
such a problem is possible by using Maxwell's 
equations. We have only: to stipulate certain 
assumptions which enable us to find suitable 
boundary conditions, assumptions, however, 
which are fulfilled to a considerable degree in 
many actual transformers. 


2. SPATIAL ASSUMPTIONS AND SYMMETRY 
CONDITIONS 


Figure 1 shows the cross section of a usual 
transformer with concentric primary and second- 
ary windings, the resistances of which may be 
neglected. We concentrate our attention on the 
primary high voltage winding which is exposed 
to traveling waves impinging from the incoming 
line. 

We assume first that the secondary winding 
of the transformer can be considered as short- 
circuited and grounded throughout. This as- 
sumption is valid to a good approximation for the 
rapid the 


transformer is in actual operation, since its 


oscillations here considered when 


secondary low voltage external network then 
7L. F. Blume and A. Boyajian, ‘‘Abnormal voltages 
within transformers,” Trans. A.IJ.E.E. 38, 577 (1919). 

SL. V. Bewley, Travelling Waves on Transmission 
Systems” (New York, 1933), p. 207 ff. 

* W. Steidinger, ‘‘Das elektromagnetische Verhalten der 
einlagigen Zylinderspule,”” Archiv f. Elektrotechnik 13, 
237 (1924), 

10 R, Riidenberg, ‘‘Performance of travelling waves in 
coils and windings,’ Trans, A.I.E.E. 59, 1031 (1940). 
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Fic. 1. Spatial position of windings in a transformer. 


forms a very low impedance. The low voltage 


winding thus acts as a screening grounded 
metallic surface with respect to any waves within 
the high voltage coils and therefore no core flux 
can be produced within the tron. 

Second, we assume that the radial thickness 
of the uniformly wound primary coil and, also, 
and 


secondary coils are small in comparison with the 


the insulating distance between primary 


diameter of the coils. If the interspace is rela- 
tively small, which is true in most actual designs, 
the problem is reduced from a cylindrical to a 
plane problem. If the primary coil is very thin, 
any considerable internal capacitance between 
adjacent turns within the primary winding is 
excluded. 

Third, we consider the primary coil as axially 
long compared with the insulating distance to 
the secondary. For most transformer windings 
of the concentric tube type, this condition is 
closely approximated. For shorter coils, as in 
the the 
be considered separately, in 


shell-type transformers, influence of 


finite ends may 
accordance with well-known methods. 

With this idealization we have removed from 
the problem the influence of any direct turn-to- 
turn capacitance, the effect of which is well 


known, while all other mutual interactions of the 
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turns, Capacitive as well as inductive, about which 
investigators are in disagreement, are retained 
in a pure form. 

Two further properties of the coils which are 
present in every actual transformer will be taken 
into consideration: the winding has a helical 
structure, as shown in Fig. 2(a), and it is sur- 
rounded by an exterior electrode represented in 
Fig. 2(a) as a grounded, concentric cylinder. 
With an oil transformer this outer electrode 
consists actually of the metallic tank; with an 
open air-insulated design it consists of the remote 
surface of the earth. In accordance with well- 
known rules, even an irregular cylindrical ex- 
terior electrode can be reduced to a circular one 
and we shall see, moreover, that its influence 
for the most part is slight. 

In accordance with our second assumption, 
a coil with a large winding diameter may be cut 
along an axial line and unrolled from a circular 
form to a plane, as indicated in Fig. 2(b), without 
introducing any considerable error. Fig. 2(c) 
then shows the insulating or dielectric spaces 
transverse to the winding in which we now 
must find the components of the electromagnetic 
field expressed in rectangular coordinates. 

We measure the x coordinate in the direction 
of the former circumference of the coil, y in the 
direction of the axis, and z perpendicular to the 
thin layer of current, which now flows in a 
system of parallel wires having an axial separa- 
tion s equal to the pitch of the helical coil. 

Since the electromagnetic state is equal at 
both edges of the circumference u along the 
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Fic. 2. (a) High voltage coil between core and tank. 
b) Coil, cut open and unrolled to a plane lattice. (c) Plane 
lattice between ground electrodes. 
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entire length of the coil, as in Fig. 2(b), we may 
repeat the representation of the unrolled coil 
periodically in the +x and — x directions, as in 
Fig. 3. By this 
with a 


we cover the entire x—y plane 
lattice of closely spaced conductors, 
the inclined lines, and we can 


now solve our problem mathematically by a 


represented by 


statement valid for any value of x and y within 
the infinite plane. 

If we draw in Fig. 3 along the helical wire of 
the coil a wave with any change in space repre- 
sented by + and — signs, we can repeat this 
wave at the left-hand and right-hand repetitions 
of the pursuing the 
course of such a wave over the real helical wire, 


conductors. Instead of 
we can follow the wave over one turn only of 
the real coil but, in addition, over the left-hand 
and right-hand continuations of this turn in the 
repeated systems of Fig. 3, which are represented 
by dotted lines. If we concentrate now on these 
conductors and look over the entire x—y plane, 
the 
wires is independent of x but depends only on 


we see that the electromagnetic state at 
the axial coordinate y. This axial symmetry of 
the actual electromagnetic field of the coil, as 
shown in Fig. 3 along the chain-dotted line, will 
greatly simplify our considerations. 

In the interspace between the inclined con- 
ductors there may be, of course, some change of 
the electromagnetic state in the x direction, but 
certainly only a small one. We see this im- 
mediately if we either think of our conductors 
as being wide enough to fill the whole space of 
the x—y plane, or if we think of the space be- 
tween the thin wires as being filled up by inter- 
leaved wires, composing a multiple thread helical 
winding in reality. 

Thus for our mathematical development we 
will consider a coil which consists of a thin 
layer, conducting perfectly in the direction of 
the wires and insulating perfectly perpendicular 
to the wires. With this idealization, we can solve 
rigorously the problem of electric oscillations in 
such a coil. We neglect hereby only the change 
of the electromagnetic field within the immediate 
neighborhood of every single wire, which is 
identical with neglecting the capacitive inter- 
action of adjacent However, by 


turns. our 


assumption of constancy in the x direction, we 
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an’ inclined lattice, 
and covering the 


Fic. 3. Development of waves on 
repeated on both sides periodically 
infinite plane. 


the electro- 


magnetic field in the y direction, and therefore 


do not restrict any variation of 
also in the direction of the inclined wires. Hence 
we include by our statement even those oscilla- 
tions which may be shorter in space than a turn 


length. 


3. FUNDAMENTAL EQUATIONS 


If at the wires the electromagnetic field varies 
to a negligible extent in the x direction, it is 
justifiable the more to assume that it remains 
the lattice of 
wires. With the condition, therefore, that 


constant at some distance from 


0/0x=0 (1) 


field within the entire 


space, we have the following Maxwellian equa- 


for every component 
tions in rectangular components corresponding 
to Figs. 2 and 3 and expressed in electromagnetic 
units. The variation of the electric intensity E 
with time ¢ is 


e OF 
=curl // (2a) 
Vg- at 
or 
e OF, OH OlT,,) . 
. - i 
vo Ot Oy Oz 
e OF,  Oll, 
= ( (2) 
Vo Ot Oz 
e OF oll, 
vo? Ot dy 


while the time variation of the magnetic intensity 


Fic. 4. Vector E of 
electrical intensity, 
perpendicular to the 
wire direction w. 








Ww 
Xx 
U 
IT is 
oH 
iv =curl i (3a) 
al 
ol 
OH, dk, dk, 
u - = 
al ay Oz 
OH, dk, 
-—s= (3) 
ol Oz 
oH Ok, 
Mu =- ed 
al oy } 


e and uw denote dielectric constant and magnetic 
permeability within the space considered, and 
vp the velocity of light in vacuum. Furthermore, 
the spatial variation of the magnetic intensity is 


OH, oll, 
div /J=0 or + =(), (4) 


oy Oz 


and that of the electric intensity is 


div E=0 or 
oy Oz 


These relations are valid within the dielectric 
Fig. 2(c). At the the 
unidirectionally conducting coil layer, however, 
the electric intensity must be zero in the con- 


spaces of surfaces of 


ductive direction. From Fig. 4 we see that the 
electric intensity E,, in the direction of the turn 
length w, which is given by circumference u and 
pitch s as 


w= (u?+57)}, (6) 
disappears only if the components in the x and 
y directions are in the proportion 

E,/E,=—s/u. (7) 
Furthermore, since all currents in Fig. 3 are 
flowing in the w direction, they cannot produce 
a magnetic intensity in that direction. Therefore, 


i) 
i) 
Nm 





we have, by the same geometric consideration 
as in Fig. 4, likewise 


H,/H,=—s/u. (8) 


We assume these relations as valid within the 
entire space. 

Equations (1), (7) and (8) state the structural 
conditions of the electromagnetic field of the 
coil and they alone determine one important part 
of the Maxwellian The boundary 
conditions at the entire surfaces of the dielectric 


solution. 


spaces determine another part and will be 


considered later. 
If we eliminate by Eqs. (7) and (8) the x 


components of the electric and magnetic in- 


tensities in Eqs. (2) and (3), we obtain from 
the first group 
s €e 0E, OH, dll, 
= in _ 
a Na 
uve Ot oy Oz | 
e OK, s OH, 
= — - (9) 
vo" Ol u Oz 
« OE. s OH, 
vo Of wu Oy 
and from the second group 
s 0H, OE, OL, 
M ca 7 . a - 
u dal oy Oz 
dH, s dE, 
— yp =—— (10) 
ot u OZ 
OH. sdk, 
— isaac, 
ot u oy ) 
Before proceeding further, we realize from 


eqs. (7) and (8) that the x components, corre- 


Fic. 5. Correlation 
of wire velocity v with 
axial velocity v, of 
traveling waves. 
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the circumferential 
direction of the coil, are very weak if the pitch s 
is only a small fraction of the circumference u, 


sponding to the field in 


as is true for most actual windings. Furthermore, 
the second and the third lines in Eqs. (9) and 
(10) show that the time variations of the electric 
as well as the magnetic field in the y and z 
directions are very small as compared with the 
variation of the field components in_ space, 
determined by the factor s/u. This means that 
the field in the 
static, although rapid time variations in the x 
direction may occur. The smaller the pitch, or 


y and gz directions is almost 


the greater the number of turns per cm coil 
length, the more closely the field approaches the 
static state. 


4. SOLUTION FOR THE FIELD COMPONENTS 


We eliminate now 0/dz from the first two 
lines in Eqs. (9) and (10) and obtain 
e/u Ss\OE, oll. 
a ae 
vy \s us dt oy 
(11) 


II 


e fu s\dll, dk. 
+ : 
vo \s us dt dy 


If we combine the first Eq. (11) with the last 


K’q. (10) we obtain 


€u w\ Ck, PE, 


Vo" sJ ae dy 


follows for J7,. 
If we combine the second Eq. (11) with the last 
Eq. (9), there results for E, and I/, the same 
relation as (12). 


and exactly the same relation 


This well-known wave equation (12) can be 
solved by an arbitrary function of the form 


E,, E., H,, H.=f(y+v,!), (13) 


which shows that the entire electromagnetic 
field is propagated undistorted in the direction y 
of the axis of the coil with a velocity given by 
the constant coefficient on the left-hand side of 


Eq. (12) as 


Vo 5S Wo 
= = | (14) 
(eu(1+u7/s*)) w (eu)? 


no = 


c 
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Fic. 6. Flow of the electric and magnetic lines of force in 
radial and in axial direction of the coil. 


Since this axial velocity depends only on the 
fundamental electromagnetic velocity vo/ (eu)? 
and on the relative pitch or slope of the turns, 
s/w we see that it is a natural velocity inherent 
in every coil given only by the slope but inde- 
pendent of all other dimensions of the coil. It is 
a small fraction of the velocity of light, confirm- 
ing the fact that the change of the phenomena 
in the axial direction is nearly static. The velocity 
is independent of the variation of the field in 
time, or of the exciting frequency. Any shape of 
the waves, sinusoidal as well as irregular, or 
impact waves of any steep front, travel with the 
same axial velocity. The front of every wave is 
propagated throughout in the axial direction 
of the coil. For coils of small pitch the velocity 
may be very low; however, if the pitch s ap- 
proaches the wire length w of a turn, the velocity 
approaches the velocity of light. For stranded 
cables, with finite pitch, forming the limiting 
case of a multiconductor coil, the velocity is still 
somewhat less than that of light; the velocity of 
light is attained only with solid wires. 

Furthermore, we can determine the phase- 
velocity of the waves along the wires. Figure 5 
shows immediately that this is larger in the 
ratio of w/s, and is therefore 

v= (w/s)v,=v0/(en)?. (15) 
Hence we see that every electromagnetic excita- 
tion travels along the wire of a transformer coil 
with the universal velocity of light, as adapted 
to the dielectric and magnetic constants of the 
surrounding material. According to the structure 
of the Maxwellian equations this result is caused 
by the fact that the and magnetic 
intensities vanish in the direction of the wires 


electric 


and that the lines of force are therefore perpen- 
dicular to the wire direction. Since the magneto- 
dielectric constant for oil transformers with 
imbedded solid insulation is of the order of 4, 


NM 
N 
w 





the wave velocity is about half the velocity of 


light in vacuum. Hence the natural frequencies, 
fundamental as well as higher harmonic, can be 
determined easily for any given length of wire 
of the winding. 

After having expressed the dependence of the 
field on the x, y and ¢ coordinates in a form of 
general validity, the variation in the z direc- 
tion remains to be determined. We eliminate 
the time variation 0/dt from the first two 
Eqs. (9) and (10) and obtain 


(1+s*,/ u*)d11,/d2—0H,/dy=0 
(16) 
(1+s°/w)dE,/dz—-dE,/dy=0. 


If we combine these relations with the second 
lines in Eqs. (4) and (5), we obtain for the 
variation of one of the components in z and y 


(w/u)ePH, d2+0H1,/dy*=0, (17) 


and exactly the same equation for each of the 
other three field components, J/7,, E, and F,. 
This relation (17) is similar to the well-known 


/ 


potential equation and therefore we can use any 
suitable solution given by the potential theory, 
except that in comparison with the potential 
problem the z dimension in our problem is 
extended by the ratio of the turn length w to 
the circumference u of the coil. For tightly 
wound coils this ratio is very nearly equal to 1; 
for loosely wound coils it is considerably larger 
and increases to infinity, when approaching the 
limiting case of a straight conductor. 

A simple solution of Eq. (17) can be obtained 
for waves which vary sinusoidally in the y 
direction along the axis of the coil. This variation 
enforces an exponential distribution in the z axis 
perpendicular to the coil surface. We see im- 
mediately that the magnetic intensities 


IT, = H,- cosh Bz: sin a(y+z,/) | 


: (18) 
H,=H,-sinh Bz-cos a(y+v,/) | 
and the electric intensities 
E, = E,-sinh Bz-cos a(y+,/f) | (19) 


E, -_ E»-cosh 8z- sin a(y+v,!) 


satisfy the differential equation (17). One such 
entire group of field components is valid within 
the internal dielectric space of the coil, and an- 
other one within the external space. Both may 
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be distinguished by subscripts i and e where 
necessary, as in Fig. 2(c). 

We have expressed Eqs. (18) and (19) in terms 
of hyperbolic functions, rather than exponential 
functions, in order to make the equations sym- 
metrical in +z and —z. Moreover, in this 
manner we can satisfy immediately the boundary 
conditions H7,=E,=0 at the perfectly con- 
ducting ground electrodes in Fig. 2(c), if we fix 
the two origins z=0 for the internal and the 
external spaces at the respective positions of 
these electrodes. 

By inserting the field components (18) and 
(19) in Eq. (17), we obtain a conditional equation 
between attenuation constant 8 in the z direction 
and wave density a@ in the y direction as 


B=(u/w)a, (20) 


valid for both the internal and external spaces. 

The field intensities of Eqs. (18) and (19) vary 
in time with an angular frequency given by the 
coefficient of ¢ namely 


w= avy. (21) 


If we excite the coil with this frequency w, the 
wave density in the y direction becomes, with the 
use of Eqs. (14) and (15), 


a=w/vy=w(w/s)[ (en)! /v9 ]=(w/s)(w/v). (22) 


From this we can easily determine the wave- 
length of the field intensities along the y axis as 


/ =2n a=(2r W)(S/W)UV. (23) 


Along the coil wires the wave-length is w/s 
times greater; namely, 


A=(2r/w)v. (24) 


These wave-lengths are given, as usual, by the 
quotient of velocity and frequency. 

The spatial attenuation of Eq. (20) transverse 
to the wires now becomes, using Eqs. (22) 
and (23), 


B=(u/s)(w/v) =(u/w)(2r/dX). (25) 


Thus we see that the extension of the field in 
the z direction, given by 1/8, is directly propor- 
tional to the relative pitch s/u of the coil, and 
inversely proportional to the exciting frequency 
w. Hence for high frequency and small slope of 
the turns, the electromagnetic field concentrates 
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more and more in the neighborhood of the coil 
layer. If, for example, as in many power trans- 
formers, the slope of the winding is s/u=1, 1000 
and the velocity of the waves is v=1.5-10!° 
cm/sec., the attenuation space constant, with 
an exciting frequency 10* c.p.s., is 

1 1 1.5-10" 


= = 240 cm. 
8 1000 27-10' 


With 10° c.p.s. this constant reduces to 24 cm 
and with 10° c.p.s. to merely 2.4 cm. 

Figure 6 shows the general picture of the 
electric and magnetic lines of force, according 
to Eqs. (18) and (19). In addition to these 
components in the y and gz directions, there 
exist also circumferential components in the x 
direction, given by Eqs. (7) and (8), which are, 
however, negligible as long as the pitch s is only 
a small fraction of the circumference uw. 

For determining the constants JJ and EF in 
Eqs. (18) and (19), we insert the field intensities 
in Eq. (16). Then we obtain for the magnetic 
amplitudes, using Eqs. (6) and (20), 


ITo= —(w/u)H, (26) 
and for the electric amplitudes 

E.=(w/u)F. (27) 
Furthermore, an interdependence of the mag- 
netic and electric amplitudes can be derived by 
inserting the field components in Eqs. (9) or 
(10). If we take E, from Eq. (19) and H, from 


Iq. (18) and introduce them into the last Eq. 
(9), we obtain 


(€/vo")v,He=(s/u)iT, (28) 
and if then we use Eq. (15), this simplifies to 
TT, = (1/v9)(u/w)(e/p) Es. (29) 


All these relations refer to the internal as well 
as to the external insulating space. 


5. BOUNDARY CONDITIONS WITH RESPECT TO 
VOLTAGE AND CURRENT 

In order to coordinate all the field amplitudes 

with voltage and current in the winding, we have 

to consider the boundary conditions at the 

surface of the coil. The line integral of the 

electric field between coil and ground electrodes 
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in Fig. 2(c) gives the voltage E at the point 
considered. Thus we have over a distance |, 
using the last field component of Eq. (19), and 
confining ourselves to the amplitude of the 
spatial sine wave 

b 


E={ E..dz=(E2/B) sinh Bd. (30) 


0 


This is valid for the internal dielectric distance b. 
Since a similar expression yields the voltage E 
for the external distance d, only with opposite 
sign, the amplitudes in both spaces are de- 
termined by 


Fo;=(B8/sinh Bb) E 


(31) 
E..= —(8/sinh 6d)E. | 


These maximum electric intensities thus can 
be computed easily by 6, Eq. (25), if the exciting 
voltage is given by amplitude & and frequency w. 

The magnetic intensities at the interior and 
at the exterior surfaces of the coil layer are 
interlinked by the conditions that the compo- 
nents /7, perpendicular to the layer must be 
equal, while the difference of the axial compo- 
nents I, is determined by the current within 
the layer. From the first condition we have for 
z;=b and —z,=d, using the last Eq. (18), 


ITs; sinh Bb = —TI/», sinh Bd = Hp, (32) 


where JI) is merely an abbreviation for these 
maximum values. From the second condition 
there follows, using the first Eq. (18) 


IT,; cosh Bb—H,, cosh Bd = 4A... (33) 


Herein we denote by A, the amplitude of the 
linear current density in the circumferential 
direction of the coil, measured perpendicular 
to the field intensities //,. Since the real current 
I within the turns flows at an angle to the x 
direction, we see from Fig. 7 that the x compo- 
nent of the linear current density A is 


A,=(u/w)A =(u/w)(NI/k)=NI/h, (34) 


where N is the number of turns in a coil of the 
axial length h, while k is the projection of this 
length perpendicular to the wires. 

Considering Eq. (32), we can express both 
amplitudes 77, by the maximum value // of 
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the transverse magnetic field through the coil 
layer, and using Eq. (26) we obtain the compo- 
nents /7; as 


IT\;= —(u/w)(7o/sinh Bb) | 


(35) 
H,,=(u/w)(H)/sinh Bd). 


If we introduce these values into the left-hand 
side of Eq. (33) and the value of Eq. (34) into 
its right-hand side, there results 


4a NI/h= —(ull,/w)(ctnh Bb+ctnh Bd), (36) 


from which H/,) can be computed if the current is 
given by amplitude J and frequency w. 

We consider now the interdependence of 
magnetic intensity /7 and electric intensity E£, 
given by Eq. (29), and insert in (29) the values 
of Eeqs. (31) and (35). The internal as well as the 
external fields give the same relation 


— Fy =(B/v0)(€/p)tE (37) 


between the maximum magnetic intensity J/, 
and the maximum coil voltage £, both in the z 
direction and thus perpendicular to the wire 
layer. Inserting Eq. (37) in (36) in order to 
eliminate /7), we obtain a definite interdepend- 
ence between current J in the turns and voltage 
FE to ground at any point of the coil. The surge 
impedance Z of the coil as ratio of voltage and 
current is determined therefore by 

w N 1/8 


Z=E/IT=4rvo(p/e)! 
uh ctnh Bb+ctnh Bd 


(X10-%in Q). (38) 
In order to transfer the electromagnetic units of 
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Z into the practical unit ohm we have to multiply 


by 10~°, as indicated in parentheses. 


6. SURGE IMPEDANCE OF TRANSFORMER COILS 


In a discussion of the surge impedance, as 
given by Eq. (38), the first group of factors on 
the right-hand side and the last fraction term 
can be considered separately. The group contains 
merely a numerical term 47, the universal ve- 
locity of light vo, the quotient yw ¢ of the insulating 
space, and a few winding data of the coil. These 
are: the ratio w, u of turn length to circumference, 
which for transformers is always equal 1, and the 
number of turns N per length / of the coil 
which is a characteristic value of the winding, 
varying from 100 to 1 turns/cm for small-size to 
large-size power transformers. Thus into this 
group of factors there enter only very few data 
of the transformer; namely, the electromagnetic 
constants of the insulating space and the number 
of turns per axial unit length. 

The second part of Eq. (38), given by the 
fraction 


1/8 
ctnh 6b+ctnh Bd 


g= 


(39) 


determines the equivalent ground distance of 
the winding. We see this immediately if at first 
small values of the arguments 8) and Bd are 
specified, which means either relatively slow 
oscillations, see 8 in Eq. (25), or small distances 
to ground. Then the ctnh’s are equal to the 
reciprocals of their arguments and hence 


1/8 1 


Zo 


=— = . (40) 
1/6b+1/8d 1/b+1/d 


In this case the equivalent ground distance go is 
given by the actual interior and exterior ground 
distances in parallel connection, and thus the 
surge impedance for low frequency or narrow 
insulating spaces becomes 
N b 
Zo =4nvo(p/e)! 
h 1+5b/d 


(X10-%in @). (41) 


We see that this value is determined by only 
two more numerical data of the transformer 
than those just mentioned, namely the insulating 
distances 6 and d, and that in the range of 
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small 6b and Bd, the surge impedance is inde- 
pendent of the frequency and of the many other 
design data of a transformer. 

If only one insulating distance, perhaps the 
interior one b, is small and the exterior one d is 
large, measured in terms of 8) and 8d as men- 
tioned above, the ctnh of 8d approaches the 
value 1 and therefore the equivalent ground 
distance becomes 


1/8 b 
g= = =b. (42) 
1 Bb — 1 1 + Bb 


Thus the influence of one remote ground elec- 
trode disappears, as already indicated by Eq. 
(40). 

A few experiments may confirm the validity 
of formula (41). A coil consisting of a single- 
layer winding of flat copper wire 1.27 mm in 
cross section and 59 m long was wound on an 
insulating tube of 36.3 cm diameter with N=52 
turns over an axial length of h=65 cm. Concen- 
trically in the interior was a closed cylindrical 
steel sheet representing a core or a secondary 
winding at a distance from coil to core b=5.4 
cm. All this was imbedded in oil and there was 
no external metal tank. From the shape of the 
switching-on voltage curves, measured with a 
cathode-ray oscillograph, in comparison with 
those of a known smooth line, the traveling 
velocity along the winding could be determined 
as v=155 m/ysec., and the surge impedance as 
Z=920 ohms, approximately. 

With an average dielectric constant «=3.5 
and a permeability 7=1, our formulae (15) and 
(41) give 


v = 300/(3.5)'=160 m/usec. 
Z =120r(1/3.5)3(52/65)5.4=870 ohms, 


which is in good agreement with the measure- 
ments. 

Another coil, also with a single-layer winding, 
but with air insulation, consisting of N=4144 
turns of round wire of 0.254 mm diameter and 
3380 m length, was wound on an insulating 
cylinder of 26.0 cm diameter and h=116.5 cm 
length, surrounding a closed copper tube of 
24.38 cm diameter at an interior distance 
b=0.81 cm. The resultant dielectric constant 
was measured as «=1.34 and the permeability 
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was u=1. The natural oscillations of the coil, 
with open as well as with closed ends, were 
measured by cathode-ray oscillograms which 
yielded a wave velocity v=260 m/usec. Varia- 
tion of an ohmic resistance, connected between 
the distant end and ground, until no wave re- 
flection occurred on the cathode-ray oscillogram, 
vielded a surge impedance Z=9020 ohms. Our 
formulae (15) and (41) give 


v= 300/(1.34)!=259 m /usec. 
Z =120r(1/1.34)'(4144 /116.5)0.81 =9350 ohms 


again in close agreement with the experiments. 

Hence we see that for coils of different sizes 
with a small as well as with a large number of 
turns per cm length, our theoretical derivations 
are verified experimentally. 

For relatively high frequency and large value 
of 8, or for large distances } and d, the values of 
the ctnh’s in Eq. (39) approach 1 and thus we 
obtain as equivalent ground distance the 
asymptotic value 


g.=1/2B=(s/u)(v/2w). (43) 


In this range the surge impedance with use of 
Eq. (15) and noting that s=h, N, is therefore, 


Ln = (2r/w) (097 /€)(w/u*). (44) 


It varies inversely as the frequency and in 
addition to this and to universal constants is 
determined for tightly wound coils by the di- 
electric constant and the circumference of the 
coil. 

If we draw a curve for the equivalent ground 
distance g, dependent on 8, we obtain a shape as 
in Fig. 8, which also shows the dependence of 
the surge impedance Z on frequency w. The 
lower and the upper asymptotic ranges are 
separated by a boundary frequency w’ which is 
determined by the intersection of the asymptotes 
of the two ranges just mentioned. The boundary 
frequency, therefore, according to Eqs. (40) 
and (43), is 


w’ =(s/u)(v/2)(1/b+1/d). (45) 


For a power transformer of average size let us 
take v=150 m/yusec., s/u=1/1000, b=2 cm and 
d=20 cm. Then we obtain 


1 1.5-10%s1 1 
w’ = ( + )=4.12-10" C.p.s. 
1000 2 2 20 


“- - 


bo 
NM 
oe | 





This is much higher than the fundamental 
natural frequency of such a transformer, which 
is of the order of 10‘ c.p.s. The surge impedance 
of transformers of usual design remains constant 
therefore for frequencies of the order of the 
natural frequency up to a first group of higher 
harmonics and decreases only for frequencies of 
the order of very high harmonics. 

Another way of expressing the boundary 
between the two ranges in Fig. 8 is to compute 
the corresponding wave-length. By Eq. (23) 
this boundary wave-length is 


N = (4nu/w)[b/(1+56/d) |. (46) 


Thus we see, that only for waves which are in the 
axial direction shorter than 47 times the result- 
ant insulating distance does the surge impedance 
decrease with frequency. For b=2 cm and d= 20 
cm 


N’ =4n-1-2/1.1=22.8 cm, 


which is small as compared with a coil length of 
about 4 = 100 cm for such a transformer, forming 
natural oscillations of \=200 or 400 cm axial 
length dependent on the end conditions. In this 
example the boundary between the two ranges 
of constant and hyperbolically decreasing surge 
impedance is near the 9th or the 17th harmonic 
oscillation. For most transformers, with coils 
which are long as compared with the insulating 
distance, the use of the lower range Eq. (41) for 
the surge impedance is therefore sufficiently 
accurate. 

However, if we remove the ground electrodes 
on both sides of the coil to greater distances, the 
boundary frequency decreases and therefore the 
upper approximation of the surge impedance in 
Iq. (44) extends to lower frequencies. The 
reason for the variation of surge impedance with 
frequency now occurring is, as seen in Fig. 6, 
that not all the electric lines of force end at the 
ground electrodes, but that many of them pass 
directly between the plus and minus sections of 
voltage and current along the coil. Thus we have 
solved here rigorously the problem of the mutual 
influence of neighboring parts of the coil in the 
case of harmonic waves of voltage and current. 
We realize that this effect has no influence on 


the velocity of traveling waves, which, according 
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Fic. 8. Dependence of the equivalent ground distance, 
and of the surge impedance of a coil on the frequency of 
exciting oscillations. 


to Eqs. (14) and (15), is independent of the 
frequency. 

There still remains a discussion of the factor 
w/u in the surge impedance, Eq. (38), which 
is the ratio of turn length to circumference. For 
tightly-wound coils wu is always 1, but for 
loosely wound coils we can express it, combined 
with the turn density V/h, and using Eq. (6), as 


wN/uh=w/su=(1/u?+1/s?)?. (47) 


Thus both pitch and circumference of the coil 
influence the surge impedance. The deviation 
from the tight-coil formula becomes greatest 
when we consider a stranded concentric cable in 
which the wires are wound in spirals. Since such 
cables always have lengths which are very large 
as compared with the insulating distances, the 
boundary wave-length \’ is comparatively very 
small, and so the surge impedance becomes, 
with Eq. (40) for the lower range of frequencies, 


1/u?+1/s*)? 
Z =4Anvo(p/€)? 
(1/6+1/d) 


(K10-%in @). (48) 


This expression agrees with the usual relation 
for cables with small insulating distances, except 
for the influence of the twist of the conductors, 
which is accurately considered here. 

If we draw out axially a helical winding so 
that it finally becomes a straight conductor, the 
ratio of turn length w to circumference u gradu- 
ally increases towards ~. Since Eq. (17) has 
shown that the z dimension of any potential 
solution must then be lengthened increasingly 
up to “, we see from Fig. 6 that the lines of 
force straighten out radially more and more 
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until for a linear conductor they become straight 
lines perpendicular to the axis. This gives the 
transition to the well-known picture of lines of 
force for waves on smooth wires, where the 
electric and magnetic forces lie strictly in radial 
planes. 

For many applications it is useful to derive 
rigorous expressions for self-inductance / and 
ground capacitance c of the coil. Since wave 
velocity and surge impedance are determined 
by these magnitudes as 


v=1/(Ic)', Z=(l/c)', (49) 
we obtain the self-inductance 
l=Z/v=4np(w u)(N/h)g (50) 
and the ground capacitance 
c=1/Zv=(€/4v9") (u/w)(h/ N)/g, (51) 


both per cm wire length of the coil. For low 
frequencies and constant equivalent ground dis- 
tance go, as in Eq. (40), these last equations 
coincide with expressions which can easily be 
derived from the dimensions of transformer coils 
by elementary considerations. 

kor very high frequencies, however, by use of 
Iq. (43), we obtain 


l= (27 /w)prw/ u* 


C=(w/2r)u*/ pv*w. 


The self-inductance now decreases with fre- 
quency, caused by the rapid spatial attenuation 
of the magnetic field perpendicular to the winding 
area, as shown by Eq. (25). The capacitance 
increases with frequency, caused by the diminish- 
ing distances between the plus- and minus- 
sections of voltage along the coil. Hereby more 
and more electric lines of force close themselves 
directly to the neighboring parts of the coil 
without touching the ground electrodes as 
represented in the exterior space of Fig. 6. 
Finally, for very high frequencies or for remote 
ground electrodes, far above the criterion (45), 
the ground does not play any part whatever. 
Thus we have found in Eq. (52) the electro- 
magnetic characteristics for a long single coil, 
not too small in diameter, situated in empty 
space without ground. 
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7. CONCLUSIONS 


We realize from all the preceding considera- 
tions that, according to the rigorous Maxwellian 
equations, the wire velocity of waves in windings 
under the conditions stated always is identical 
with the fundamental velocity of light in the 
insulating space. 

Under conditions prevailing for the most part 
in actual transformers the surge impedance of 
the winding is, under the conditions stated, inde- 
pendent of the frequency of oscillations and 
dependent on only three characteristic data of 
the winding, namely turn-density, insulating 
distances and magneto-dielectric ratio of the 
insulating space. For extremely high frequencies 
or for very large insulating spaces on both sides 
of the coil, the surge impedance decreases 
gradually with frequency. Furthermore, the 
relative pitch of the turns, if large, enters the 
result and determines in every case the attenua- 
tion of the electromagnetic field perpendicular 
to the axis of the coil. 

As an important result of the rigorous theory 
we have found that, under the stated conditions, 
tightly wound coils perform, up to a boundary 
frequency, along their wires with respect to 
traveling waves exactly like smooth lines. In 
the direction of their axis, however, they perform 
nearly statically. Perpendicular to this axis the 
spatial field of harmonic oscillations appears, 
above that boundary frequency as exponentially 
attenuated. 


SPECIAL SYMBOLS 


b=internal insulating distance 
d=external insulating distance 


equivalent ground distance 
h=axial length of the coil 

= pitch of the helical coil 

u =circumference of the coil 


“ 


v=velocity of waves 


~ 
: 
\| 


wire length of a turn. 
= linear current density 


| 


N=number of turns of the winding 


Z =surge impedance 

a=wave density in space 
8=attenuation constant 

\= wave-length in axial direction 
w=angular frequency. 





Steady-State Solutions of Electromagnetic Field Problems 


I. Forced Oscillations of a Cylindrical Conductor 
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rhe object of this investigation is a study of the current distribution in or on the surface of a 
conductor and its associated field under the influence of a localized e.m.f. Steady-state solutions 


of the field equations are found for conductors of simple geometric form. The results clarify 


many electromagnetic problems involving localized sources, especially in the u-h-f region, 


for which ordinary circuit theory fails to give a satisfactory quantitative explanation. Part I 


treats the problem of a straight cylindrical conductor and shows the relation of the principal and 
complementary waves to the nature of the exciting field. A driving point impedance is calcu- 
lated for the case of an external field applied over a vanishingly short section of conductor. The 
driving point impedance is infinite for a conductor of infinite length and perfect conductivity. 
Likewise the case of a conductor of finite length bounded at either end by an infinite, perfectly 
conducting plane is discussed. This problem bears a direct relation to that of a hollow pipe 


excited\by a linear antenna. 


HE object of this investigation is a study of 

the current distribution in or on the surface 

of a conductor and its associated field under the 

influence of a localized e.m.f. Solutions of electro- 

magnetic field problems are ordinarily con- 

structed from wave functions satisfying a homo- 

geneous set of equations, which in rationalized 
m.k.s. units can be written 


VY XE+iwuH = 0, (I) 
VYXH—(e+iwe)E=0. (II) 


It is assumed here that all quantities contain the 
time in the form of a factor exp (iwt). The field 
intensity E is measured in volts/meter, H in 
ampere-turns/meter, and conductivity @ in 
mhos/meter. In free space wy=4r K1077 henry 
meter, €9= 10~*/ 36” farad/meter. When the usual 
boundary conditions are applied in an appropri- 
ate coordinate system, (I) and (II) reduce to a 
set of homogeneous algebraic equations whose 
determinant must vanish. The complex roots of 
this determinantal equation fix the frequencies 
and damping of the allowed or free modes of 
oscillation. The amplitudes, however, are com- 
pletely arbitrary and will be determined by the 
nature of the exciting source, to which (1) and 
(Il) give no clue. 

Exact steady-state solutions of electromagnetic 
problems have been found in a number of simple 
cases, such as that of a plane wave incident upon 
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a cylinder or sphere. The equations relating the 
coefficients of the various particular solutions are 
now inhomogeneous; hence the amplitudes of the 
excited modes are uniquely determined and the 
system oscillates with the frequency of the 
impressed field. In such problems the diffracting 
body is entirely enveloped in the incident wave. 
Quite as numerous as these problems of purely 
optical interest are those in which the impressed 
field is applied as a local e.m.f., confined es- 
sentially to a small domain of the conductor. The 
applied field intensity will be designated by the 
vector E’ measured in volts/meter. E’ is a 
specified function of position and time but is not 
necessarily analytic. It represents the intensity 
of any external force, of whatever origin, acting 
on the free charges of the conductor. It may be 
introduced by a mechanical or chemical agency, 
by the motion of an external magnet, or by a 
shielded transmission line. In the latter case, 
which is the common and practical one, it is 
obvious that the entire problem is modified to 
some extent by the presence of the line. It is 
reasonable to assume, however, that in most 
cases the effect of the leads on the current 
distribution in the principal conductor or upon 
the field at a distance is small and that the 
perturbation can be taken into account by a 
second approximation. 


It will be assumed in the following that a 
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conducting body is embedded in a homogeneous 
dielectric medium and that the external e.m.f. is 
applied within or on the surface of the conductor 
only. The grounds for the latter assumption are, 
perhaps, open to some question, since an applied 
field at the surface of a conductor must inevitably 
extend into the neighboring dielectric and thereby 
effect its polarization ; but the effect on the total 
field of such additional dielectric displacement 
currents is small relative to that of the con- 
duction current, and its neglect seems entirely 
warranted. Hence the field within the dielectric is 
governed by Eqs. (I) and (II). If the conductor 
is metallic the displacement current term is 
negligible, and since the current density J is 
proportional to the resultant field intensity, 
J=o(E+E’), we have to satisfy in this domain 
the inhomogeneous system 


VY XE+ie~H =0, (1) 
VY XH-—cE=cE’. (Il-a) 


Analytic continuation of the solutions from one 
domain to the other is insured by the boundary 
conditions which call for a continuous transition 
of the tangential components of E and H. Note, 
however, that the transition of E’ across a bound- 


ary is arbitrary and in general discontinuous. 
n X< (E.—E,) =0, n<(H.—H,)=0, (IIT) 


where n is a unit normal vector drawn from 
medium (1) to medium (2). The normai will 
always be drawn outward from a metallic surface 
into the dielectric. 

There is, of course, nothing new in_ these 
equations and they are discussed in most of the 
older texts on electromagnetic theory.' The pur- 
pose of the present work is to show that in many 
cases they are easily integrated and that the 
solutions can be of greater practical interest than 


those giving simply the transient oscillations. As 


‘a first example we consider the propagation of 


waves along a single, infinite conductor of circular 
cross section. The classical treatments of this 
problem are due to Rayleigh and to Sommerfeld. 


M. Abraham, Theorie der Elektrizitat (Teubner, seventh 
edition, 1923); O. Heaviside, Electrical Papers (Macmillan, 
1892). 
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It was shown by Sommerfeld that the system 
admitted a principal wave propagated with 
negligible damping and a velocity approaching 
that of light, together with an infinite series of 
complementary waves associated with damping 
factors so great as to result in immediate extinc- 
tion. In dielectric wires or hollow pipes, on the 
other hand, no principal mode occurs and for 
sufficiently high frequencies the complementary 
waves are propagated with relatively little 
attenuation. These are free modes of oscillation 
whose amplitudes are determined by the initial 
conditions and which are responsible for reso- 
nance phenomena. A ‘“‘surface impedance” was 
also defined by Rayleigh and used as a basis for 
the discussion of the alternating-current resist- 
ance of linear conductors. It differs from the 
driving-point impedance to be observed at the 
point of application of an external e.m.f. 

We shall assume that E’ is parallel to the axis 
of the conductor and is uniform over the cross 
section. In Fig. 1 the axis of the conductor is 
shown coinciding with the z axis of the coordi- 
nate system and the impressed e.m.f. is intro- 
duced in the neighborhood of z=0. Particular 
solutions of (1) and (11) in cylindrical coordinates 
r, @, 2 appropriate for the dielectric domain 
r>a can be taken directly from Sommerfeld’s 
discussion.” 


Thus for r>a 


E,=AH, (por) 
E,=(iu/ps)AH, (por) (e°"""*) — @D 
H,, = (iwes/p2)AH, (por) 


p2=(Ro?—Uu*)*, ko*® =w ope, (2) 


where JI, (por) is a Hankel function of th 
second kind insuring proper behavior of the field 
at infinity. 


Inside the conductor 
E,’=E,'=0, E,'’=E'(s), (3) 
where E’ is a prescribed function of z and f. 


2 A. Sommerfeld, Ann. d. Physik 67, 233 (1899); also i 
Riemann-Weber, Differentialgleichungen der Physik, Vol. 
Vieweg, 1935) eighth edition. 
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. Section of infinite conductor of circular cross section 
a embedded in dielectric medium &, ps. 


Equations (1) and (Il-a) now reduce to 


OE, OF, . 
+10 TT =), 
az or 
oll ; 
+o,k,=0, (4) 
Of 
10 
ril ) o\l ok’. 
kliminating #£, and /7, one obtains 
10 Ok, Er, . 
r + ~1wpyo lt 
ry or or Os" 
’ o?k’ : 
lwuyo lt — - (5S) 
ae? 
if now we place 
F=E,+E’ (6) 


and note that the field is finite on the axis, it is 
apparent that for r<a particular solutions of (5) 
must be constructed from the functions 


F=BJ,(pyrei*' 
k,? = 


’ /) 


pi _ (k;° _ u- 


- LW 10 | (8) 


where Jo(pir) is a Bessel function of the first 
kind. 

The impressed field £’ is a prescribed function 
of z. For the moment we shall assume only that 
E’(s) and its first derivative are piecewise con- 
tinuous and that the integral 


Ie’ \dz 


exists. Under these circumstances the impressed 
field can be represented by the Fourier integral 


E! =G(z)ei*! 


oF s 


(e'*! 2n) | du | da G(a)e’“{*-*), (9) 
x 


. L . 


The coefficients A and B of the particular 
solutions (1) and (7) are now considered to be 
functions of a 


parameter a. Upon integrating 


with respect to u and a@ the field at any point is 
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represented as a superposition of cylindrical 


waves. Thus for r>a 


— 


e wt os oe 2 
E,= du | da A(a)Ho (Nor) e ” 


: (10) 
e't a = A(a) 9 
H,= 1WE>o du da IT; (per) es" . 
OF San Oe Pe 
while for r<a 
e wt a7 es 
E,=- du | 
saan Paes 
X dal Bla) Jo(pir) —G(a@) e : 
11) 
e wt os es Bla) 
H.=e; du | da Ti(pirye(a 
2n , x e x Pr 


These equations are next subjected to the 


boundary conditions (III) at the surface r=a. 


There results the algebraic svstem 


B J«(pya)—A Hy (p22) =G(a), 


(12) 
0} : IWE2 (2) 
B J\(pya)— A Hy, (poa)=0 
Pi Pe 
from which we obtain 
x Gla) 1 
B . 
o,a J\(x) Z(u) 
x Jo(x) y Ho (y) 
Z(u)= —- 13) 
oa J\(x) tweoa IH (y) 
a * 
x=piad, yvV=p2a. 


In the absence of a driving force G, the natural 
modes of propagation are given by the roots of 
Z(u) =0, exactly the condition found by Sommer- 
feld. This condition has an interesting physical 
interpretation. By definition the wave impedance 
of a cylindrical wave in the radial direction is 

-E./H,. The rule for the algebraic sign is that of 
the Poynting vector. The ratio is positive when 
the components of E and H follow in cylindrical 
order.* It is clear that apart from a sign the first 
term of Z(u) is the radial impedance of the 
cylindrical field within the conductor, while the 
second term expresses this impedance in the 
external medium. Z() is the difference of the two 


impedances at the boundary. The free oscillations 
S. A. Schelkunoff, Bell Sys. Tech. J. 17, 17-48 (1938) 
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of any system of bodies embedded in a homo- 
geneous medium follow from the condition of 
equality, or “‘matching,”” of the normal com- 
ponents of wave impedance across the bounding 
surfaces, expressed in the present instance by the 
vanishing of Z(u). 

Let us now calculate the current passing 
any cross section of the conductor. Since 
J.=0,;(E.4+-E’), we find for the current J at any 
point Z 


7 


g(t) 
I =ae (2x)! | ——e~ “dy, (14) 
J. 28) 
where 
7 
g(u)=(1 (2n)!) | G(uye'*da (15) 
e L 


is the Fourier transform of the applied field G(s). 

The zeros of the impedance function Z(1) 
occur at the values of u corresponding to the 
propagation constants of the principal and com- 
plementary natural modes, and all are complex 
provided oa; is finite. Z(u) is symmetric in u and 
the roots lie in the first and third quadrants, with 
none falling on the real axis. It follows, moreover, 
from the asymptotic behavior of the Bessel and 
Hankel functions that as ux, Z(u)—< as 
(1/o,+1/tweo)u. Hence the path of integration 
along the real axis may be closed by a semicircle 
of very large radius in the upper or lower half- 
plane according as 2 is greater or less than zero. 
This is illustrated in Figs. 2(a) and 2(b), the 
dots indicating schematically the location of the 
Zeros. 

We assume the transform g(u) to be ana- 
lytic everywhere within the closed contour and 
properly bounded at infinity. Let ~, designate the 
roots of Z(u) =0. Then 


I= +a(2z)'e'*'2ri > Dy, (16) 


where b,, is the residue at the pole “, and the plus 
or minus sign holds as z is less or greater than 
zero. 

To determine the residues we need only the 
first term in the Taylor expansions of Z(1) about 
the poles. 

dZ 
Z(u)=(u—U,) - To, (17) 
du u=Un 
and hence 


weet g(u,je'“" 
IT =a(2r)*e'*'2mi > , 18) 


n (dZ/du)u 
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Thus the source excites all modes and these 
travel along the line with the propagation factors 
determined by Sommerfeld’s analysis. The ampli- 
tudes are uniquely determined by the form of the 
exciting function g(u). It will be recalled that if 
the radius of the conductor and its conductivity 
are not too small, the propagation factor of the 
principal wave is given approximately by 


Uy ke = w(€ou2)' = 27 X. (19) 


The effect of a finite conductivity is to introduce 
a very small imaginary part and a corresponding 
attenuation. The attenuation of the comple- 
mentary modes, on the other hand, is exceeding], 
large and so for all practical purposes (18) reduces 
to 


a(R») 
' . S é 

[ =a(2r)?2x1— 
(dZ /dut) uxke 


e'*! 22 (20) 


To calculate the denominator in (20) note that 


p20, kii\>>ke, xk a>1. (21) 
If we write kj}=a,;—7f;, it follows from the 
asymptotic representations of the Bessel func- 
tions that 


_  Jo(x) ; ; —— 
lim = lim cot [a(a,;—i8;—7/4a) |=1, | 
Boxe J\(x) By 


Nm 
Nm 


while from the first term of the series expansions 
of the Hankel functions in the neighborhood of 
y=0 it is easily shown that 

















Hy (y)/Hy (y)-y In (yyi/2) as yO, (23) 
where y =1.781. From this it follows that 
dZ /du—2ia(ps/e2)? In (yyi/2) as uke, (24) 
z<O z>O 
a b) 
ic. 2. Closed contours of integration in « plane. 
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and the current at any point on the line is 


:; g(R) 
[ =22(7e2/2y2) ee 05-888 8 
In (yyt/2) 


nN 
uN 


In case the external dielectric is air, the first 
factor reduces to 2r(€9/uo)'=1/60. 

For the sake of example consider the case of an 
applied field whose form is represented by the 
“impulse function.” 


‘ Eo é . (2d 
G(z) = : (26) 
(2r)? d 
with the property 
G(z)dz= Eo, (27) 


e x 


where Eo is a constant and d a parameter. Its 
Fourier transform is 
Eo 


g(k») = gf te 8, (28) 
(27)! 


By definition the applied voltage is 
V= | E’dz= Exe, (29) 
The ratio of this voltage to the current at the 
point <=0 depends in general on the parameter 
d, which measures effectively the length of the 
segment to which the driving force is actually 
applied. In the present instance 


V, 1(0) =120 e@ */2 In (yy2/2). (30) 


In the limit as d—0 the region of application 
reduces to a point and the ratio represents the 
true driving-point impedance of the line. 


Z;,=lim V/I(0)=120 In (yyt/2) ohms, (31) 


i--0 


assuming the external medium to be air. Practi- 
cally, it is only necessary that the region of 
applied field be small relative to the wave-length 
to justify the concept of a voltage V applied at a 
definite driving point z=0. In this limit the 
impedance is independent of the exact functional 
form of the applied field. This is easily verified for 
the case of a “‘rectangular’”’ distribution defined 
by 


G(z)=Ey) when z| <d, 
‘ , (32) 


G(z)=0 when z| >d. 


Li) 
w 
_ 


The ratio of V to the current at the midpoint is 
now 


V/I(0) =120 (dko/sin dke) In (yyt/2), (33) 


which reduces to (31) in the limit d=0. 

One will note that u is exactly equal to k2 and y 
equals zero only in the case of a perfectly con- 
ducting wire. In these circumstances the driving- 
point impedance of the system is infinite, a result 
entirely comprehensible from the fact that the 
wire is of infinite length. Sommerfeld’s solution 
has shown how currents, established in an 
unspecified manner, are propagated. We now see 
how these currents are related to the source and 
it appears that a finite current in a perfect con- 
ductor of infinite length cannot be established by 
a finite voltage. If the conductivity of the wire is 
finite, the input impedance is large but bounded. 
A current enters the wire and the inflowing 
energy is eventually dissipated in heat. 

It is of some interest to consider the case of a 
wire of finite length, although it has no direct 
bearing on the free antenna problem to be dis- 
cussed later. The procedure followed above will 
give us, by a slight modification, the driving- 
point impedance at the center of a linear, vertical 
antenna bounded at either end by horizontal, 
perfectly conducting planes. The problem is thus 
related indirectly to the establishment of waves 
in hollow pipes. 

The applied field is now represented by a 
Fourier series in place of an integral. Let the 
transverse, perfectly conducting planes be located 
at z=l and z= —l. The tangential component E 
must vanish over these surfaces and hence in 
place of (1) when r>a we construct a solution of 
the form 


© nT ’ 7 
E,=e'' > A, Il, (per) sin —(z+l), 
n 0 2p. 2 
a » nT 
E,=e'*' > A, Hy (per) cos —(s+)), 
n=0 ? 
(34) 
2 An 2 nr 


IT, =iwe.e'*' > 
7 0 Pe 


4nr? n*r*\! Qt n\\?\? 
no (EY -M(-(2)): 
? 4/? nN 4/ 


where \ is the wave-length in the diclectric. The 


Hy, (per) cos 31 (s+), 
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applied field is expanded in the series 


‘ 


. ~ 7 nT - 
E’ =e‘t > C, cos Thelen (35) 


so that the field within the wire is represented by 


nT 


E.=e'*" ¥ [BuJo(pir) — Cu] cos > (+0), 


(36) 
» B, nr 
IT,=o\e'*' > —Ji (pir) cos —(s+l). 
n=O Pi 2/ 
The boundary conditions at r=a lead to 
* Ge 7 ‘ 
B,= (37) 


a;a Ji(x) yo 


in which Z, is formally identical with (13) but 
with uw restricted to the discrete set of values 
nx /2/. The current at any point on the line is now 


a oe nT 
I =2rae'“' >> cos (s+/). (38) 
; - 


n=0 


Assume next the rectangular distribution of 
applied field defined in (33). E’ is then a sym- 
metric function of z, so that all odd coefficients of 
(36) are zero, while for the even coefficients one 
obtains 


ie=%. 2. ¢+-}: (FR) 


The input impedance at the central point of the 
wire is again found by taking the ratio of the 
applied voltage to the current at z=0 and passing 
to the limit as z—0. 


l | 2 1 
Z.= 1 +2 yy 3 (40) 
Tra |; Ls. 


or for the admittance 


Tra 4 
y= (142 Van), 41) 
l n=1 
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where Y,=1/Z,. Let us assume for the moment 
that the conductivity of the radiating cylinder is 
infinite. Then Y,, reduces to 

¥ 1 i 

Y2,= —- (€2/ 2)? 
(1—8,7)? 


H, ((2ra/d)(1—B,”)* ] 
x , (42) 
H® [ (2ma/d)(1—B,,7)? ] 

where 8, =nd/2l. If B,>1, the ratio of Hankel 
functions is a pure imaginary. Hence if \> 21, all 
admittance terms with the exceptions of Vo have 
susceptance components but no conductance. 
The complex power associated with these terms 
is purely reactive. If now AX is decreased, a critical 
wave-length is reached at which VY» assumes a 
nonvanishing conductance. Further decrease of \ 
introduces successively in step-wise fashion the 
conductances of higher modes, and with each 
entry there is a corresponding increase in real 
power input representing outward radiation of a 
higher order wave. It will be noted in passing 
that the sum of these terms fails, in the present 
instance, to approach a limit, due to the as- 
sumption of a voltage concentrated at a mathe- 
matical point rather than over a small but finite 
segment of the cylinder. This matter will be 
referred to later. 

The term Yo has a nonvanishing conductance 
at all wave-lengths and represents the principal 
wave. Its field is strictly two dimensional, there 
being no variation along the z axis. The com- 
ponent £, is zero everywhere and energy flows 
radially outward from the cylinder. Such a mode 
is possible because of the absence of any further 
boundaries transverse to the conducting planes. 
It will be recalled that the principal wave is 
completely suppressed in hollow wave guides of 
finite cross section. 

In case the conductivity of the cylinder is 
finite the term xJo(x) /J;(x) in (13) contributes a 
small conductance to each admittance com- 
ponent at all wave-lengths, accounting for the 


dissipation of a certain amount of energy in heat. 
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Che method discussed in Part I for the determination of steady-state solutions of electromag 


netic problems is applied to a finitely conducting sphere as the simplest illustration involving a 


body of finite dimensions. The external e.m.f. is applied across an infinitesimal strip at the 


equator. The induced field appears as the sum of an infinite number of partials or modes. As the 


frequency is increased these modes become successively dominant. There are a series of reso- 


nance points. The real part of the driving-point admittance rises step-wise with increasing 


frequency and the behavior of the system is notably different from that of a normal circuit 


Expressions are given for the power dissipated in heat and in radiation and it is shown that the 


complex power input at the driving point is equal to the integral of the complex Poynting vector 


over the surface of the sphere. 


Hk preceding paper' has dealt with the 
forced oscillations of a cylindrical con- 
ductor under the influence of a local applied 
e.m.f. In the case of a cylinder of infinite length 
there is no radiation loss and the real part of the 
driving-point impedance is due solely to the finite 
conductivity of the material. In the case of a 
vertical cylindrical conductor bounded at either 
end by horizontal, perfectly conducting planes of 
infinite dimensions, the radiation is mainly two 
dimensional. In both cases the conductors extend 
to infinity and the significance of the results is in 
some respects obscure. In this present paper, we 
consider certain forced oscillations of a metallic 
body of simple geometric form and finite dimen- 
sions——a sphere—and discuss in more detail the 
resonance phenomena, the radiation impedance 
to be measured at the driving point, and the re- 
lations between radiation power, conduction loss 
and power input. The foundations are laid for a 
discussion of the generalized case of a prolate 
spheroid of arbitrary eccentricity.’ 
If E’ is the intensity of the applied field, one 
has at any point in a homogeneous medium 


VY XE+iwyuH = 0, (1) 
VXH—(o+iwe)E=cE’. (2) 
Any influence of E’ upon the polarization of di- 


1J. A. Stratton and L. J. Chu, J. App. Phys. 12, 230 
(1941). 


2L. J. Chu and J. A. Stratton, J. App. Phys. 12, 241 
(1941). 
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electric material has been neglected, since it is 
understood that the external field is to be applied 
directly to conductors, with the result that 
polarization currents will be masked by neighbor- 
ing conduction currents. From (1) and (2) it 
follows that E and H satisfy 


VXV KE+F°E=iwycE’, (3) 
—-VXVXH+FH=-—cVXE’, (4) 


where k? = w*eu —twyue. 

Consider a metallic sphere of finite conduc- 
tivity embedded in a dielectric of zero con- 
ductivity, with o;>wee It will be assumed that 
E’ is applied uniformly over the equatorial plane 
and in a direction normal to this plane. This 
corresponds to the assumption made in I. 
Strictly, such a field does not exist, since it is 
irrotational. Due to the high conductivity of the 
metal, however, the total field, which is the 
resultant of the applied and induced fields, is 
confined to a thin laver near the surface, and the 
resultant current distribution is governed essen- 
tially by the value of E’ within this layer alone. 
Its thickness is small .relative to the wave- 
length, so that E’ may be assumed uniform 
over the layer cross section and consequently 
v XE’ =0. To this approximation the field equa- 
tions are homogeneous in the dependent variables 
H and E+E’, within the conductor as well as 
within the external dielectric. Assuming, there- 
fore, that E’ is applied in the z direction and is 
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symmetrical about the z axis, one finds for the 


spherical components of the induced field : 
In the dielectric, medium 2, R>a, 


_ (2 1 
H,=e'*' > A,h, (k2R) P,(cos 6), 
n=1 


— giwt . " 
Er=- : > Anh, (RoR) 
lweoR sin @n=1 ' 
(5) 
, l l 
X {cos 6 P,(cos 0) — P,,.:(cos 8) |, 
Eo=i(u2/e2)'e**! S A Lh ns(ReR) 
n=1 
(2 l 
—(n/koR)h,, (keR) \P,(cos 8), 
while in the metal sphere, medium 1, R<a, 
>. ! 
H,=e'*' Z B,, jn(RiR) P,(cos 6), 
n=1 
— 
Er=—— > Ban j,(RiR) 
o1R sin 6 n=! 
(0) 


1 1 
<I cos 9 P, P, ' Ep’, 


Ee (Why 101) ” si > B,| } 


1 


(n/kiR)7,(RiR) JP ,1— Ee’, 


l ki R) 


where : 
ky = | 1Wj10 1) P Ro =wl€op2)?, 


h.. (RoR) =(m/2koR)!H,,.,( RoR), 
i (kiR) =(2/2kiR)'J,.\(RiR). 
P,,(cos 0) = —(d/d0)P,(cos 6). 


The applied field E’ is expanded in a series of 
spherical harmonics as 


5° C,(R) P,(cos @). (7) 


— 


E,’ = {Wy 1a) *e we 


the coefficients C,,(.R) being functions of R alone. 
The conditions on the continuity of the tan- 
gential components of E and H at the boundary 
R=a then lead to 


wus \? Cte Jnl kia) 
B,=-—{ - —, A,=- B,,, (8) 
10} i,(k,a) Z, h® (Roa) 


Lwpe| hy, —1( Rea) n 


Z,= - 
ko 9 Roa 
h, (Rea) 
1p) ; = kya) n | 
’ : Q 
ky Jn( Ria) kia 
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The quantity | 4 expresses the difference be- 
tween the radial wave impedances Ey///,4 inside 
and outside the sphere, and free oscillation of 
the sphere occurs when these impedances are 
matched, or Z,=0, an equation whose roots 
determine the characteristic frequencies of the 
natural modes.* Because of the extremely large 
value of o in metals one has, even for very small 
values of radius a, the condition k,a\>1, so 
that 


(topr/Ri){ jn—a( Ria) / jn(Rkoa) —n/kya | 


S 


= — (twys/o1)’, 


a very small quantity, and hence 


Z nt (po/e2) Th, (Roa) /h, (koa) —n/kea | 


+ (lwp a\)?. (10) 


Now koa=2na/\, where X is the wave-length 
measured in the dielectric. It is therefore a 
measure of the curvature of the sphere. When 
kywi— «x, the terms within the square brackets 
approach —1, and Z,, reduces to 


oe/€2)? + (twui/ oi)’, 


the difference between the intrinsic impedances 
of dielectric and conductor. One will recall that 
the wave in the dielectric propagates away from 
the center of the sphere while that in the con- 
ductor supplying the conduction loss travels 
toward the center. Looking outward from the 
center, the intrinsic impedance of the conductor 
is —(twu;/a;)?. When koa<n, the terms within 
the square brackets are essentially real, and 


be n koa (koa)?* 
Z =i( ) | + + 
€o / koa = 2n—1 2n—1)*(2n—3) 
2(kea)® 
+ + 
In —1)*(2n—3)(2n—5 | 


1 : 
-( ) . (11) 
0} 


The @ component of current density in the 
metal is according to (2). 


0 


1 
Jo=a(Eot+ Ey’) = RIT.), (12) 


R AR 


P. Debye, Ann. d. Physik 30, 57 (1909). 
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Radiation conductance of first few modes as 
functions of 2wa,X. 


Fic, 1 


and hence the total current crossing the surface 
#= constant ts 
olga . 
I= | | JoR sin 0 dR do 
=2xa sin 6 Ha! rea. 13) 


This determines the current distribution on the 
sphere. In particular, the current crossing the 
equatorial plane §6=90° is 

es (al) 


“(wu io;)* > P,*(0). 14 
1 Z 


I == 2ra a 


In passing, one sees that whatever the value of 
E’ in the interior of the sphere, the resultant 
held E+E’, the 
density, to the high 
conductivity except within a thin layer at the 
This skin effect 
spherical Bessel functions of complex argument. 
Thus the assumption that V XE’=0 is justified 
so long as E’ is effective only within the thin 


and therefore also current 


vanish everywhere due 


surface. is governed by the 


layer at the surface. 
In virtue of the assumption that E’ is irrota- 
tional, one may write 


(15) 
where V is the potential, or voltage with respect 
to some reference plane. Consequently 


R)(aV/00), Ep’=—O0V/0R. (16) 
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hic. 2. Radiation susceptance ol first few modes as 


functions of 2ra_X. 


This voltage is applied at the equator across a 
small strip or band of width aAé@ on the surface. 
On the surface of the sphere Er’ =0 and hence 


vv  w—Aé r+A0 
E,’ = — tor < @: ; > 
aXdé 2 2 
17 
r—A0 , 
; | 
EF,’ =0 for 6 
r+AA 
) 


In virtue of the orthogonal properties of the 


P..‘(cos 6) the coefficients C,(a) are found to be 


lo1\*? (2n+1) 
C,(a)=—- ra(O) Ve, (38) 
wuis 2an(n+1 


where 
n—1 . 
P,'(0) =1"—!2! n'/\( Ht 
) 


when v7 is odd, and P,'(0)=0 when n is even. 


Upon introducing this into (14) one is led to : 
2n+1 [P,'(0) } 
Ib=arV > (19) 
n(n+1) } a yi 


while for the input admittance at @=90° one 


obtains 
Ty 2n+1 [P,'(0) }? 


Y,=—=r) 20) 
V’ n n(n+1) Ye 


a ae 


JOURNAL OF APPLIED PHYSICS 








hh 


eee 


The input impedance Z; is the reciprocal of Yj. 
The summations are over odd values of n. 

This expression for the input admittance ob- 
viously suggests for the equivalent circuit of the 
sphere an infinite system of parallel, resonating 
networks, the admittance Y, of each branch 
being represented by one of the modes of the 
sphere, or 


2n+1 [P,'(0) }? 
Y.=f —, (21) 
n(n+1) Pe 


The real and imaginary parts are respectively 
conductance g, and susceptance b,. If one 
assumes the conductivity of the sphere to be 
infinite, the admittance Y, reduces to a pure 
conductance as kea— x, and contains only sus- 
ceptance when koa <n. For intermediate values 
of koa the admittance is complex. 

In Figs. 1 and 2 the real.and imaginary parts 
of the admittance Y,, are plotted as functions of 
2ra/X. The conductivity of the sphere is assumed 
infinite. Since there is then no loss in the metal, 
the ordinates represent radiation conductance 
and radiation susceptance. The conductance 
associated with each mode or harmonic term 
rises to a peak at resonance, then drops off 
gradually to a constant value. The total con- 
ductance rises step-wise with increasing a/X as 
shown in Fig. 4. Thus as the applied frequency 
is increased the harmonic terms add their contri- 
butions successively, each new term giving rise 
to a corresponding increase in the radiation. 
Inversely, for a fixed wave-length an increase in 
the radius of the sphere increases its effectiveness 
as a radiator. 

From Fig. 2 it appears that the radiation sus- 
ceptance of each successive harmonic rises to a 
maximum near resonance and then falls rapidly 
to zero as the wave-length decreases. The sus- 
ceptance is always positive, showing that the 
admittance is capacitive. It is evident from the 
figure that the sum of the susceptances fails to 
approach a limit. The ratio Y,/¥,—-2 in facet 
approaches unity as n— ~«. This is a consequence 
of our assumption that the impressed voltage is 
applied across a surface of discontinuity at 
§=n/2, a singularity implying an infinite value 
of field intensity. Actually the voltage will be 
applied across a strip or gap of small but finite 
width, with a corresponding decrease in the 
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magnitudes of the higher order harmonics. The 
general form of the conductance and susceptance 
curves will not be materially altered, but their 
sum then converges, as it obviously must on 
physical grounds. In this case the susceptance of 
each harmonic is small except in the neighbor- 
hood of its own resonance point. The total input 
susceptance will vary only slightly with the ratio 
a/X as compared with the total input conduc- 
tance. Since the susceptance of each harmonic 
term is always positive, under no condition will 
the input admittance, or the input impedance 
become a real quantity. 

It is apparent that the resonance phenomena 
associated with a radiating sphere differ funda- 
mentally from those encountered in conventional 
circuit theory. In this connection it will prove 
interesting in a following paper to observe the 
transition of resonance curves from the case of a 
linear oscillator represented by a prolate spheroid 
of unit eccentricity to that of a sphere of zero 
eccentricity. 

Figure 3 illustrates the amplitudes, or absolute 
values, of the first few admittances. 

The power input into the system is equal to 
one-half the product of the applied voltage and 
the conjugate of the current at the driving point. 
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Upon rationalizing Z,* in the denominator one 
obtains 
- 2n+1 [| P,'(0) 
P, Vi2>— 
2 n(n+1) z 
| fue h,,-\( Rea) n Lap} 
X42 + . tas 
L \e. | koa ao, J | 
h, (k x1) 


The real part of P; must equal the sum of the 
power radiated from the system and that dissi- 
pated in the conductor. The radiated power P, is 
given by the integral of the radial component of 
the Poynting vector over a sphere of radius 
concentric with the 


than a and 


greater con- 
ductor. 
P.=4| de | E,I1,*R? sin 6 0 

(“) rR? ) 2n+1 [ P.'(0) 

l 

€. 2a* ™ n(n+1) : 
h,, (kea)Z, 
l 9 i ) 
Kh» (RoR) hy 1 ReR) h, (RoR) }. (24 
koR 
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Since 


pa l 2 
h, (koa) *=h,, (koa) h, (koa) 
on the surface of the sphere (R=a), we have 


2n+1 | P,'(0) 


— ° 
n(n+1) pa 





=o 


h,, (Rea) nN 
x - 25 
Roa 
h, (Rea 
The power dissipated as heat in the metallic 


conductor is evidently 
P,=P,—P,=(14+1)—(wp/201)? V 
? 


2n+1 
XL 
n(n-+1) 


[P,,'(0) }? 


Zz ,i|* 


(26 


One will note that as the conductivity becomes 
infinite, not only does the real part of the input 
power P; equal that of the radiated power P,, 
but also the imaginary parts are equal to each 


other, the wattless component of P, being ’ 
measured at the surface of the conductor. The 

real part of P,, from which one finds immediately J 
the radiation resistance is s 


Rel P J=<(es uo)? Vi? 
2n+1 [P,(0 
[koa h 


Watts, 27 


— 7 
nin+1) koa) 


the prime in the denominator indicating differ- 
entiation with respect to the argument koa. It 


can be shown that 
VY *=2P,;/V 28) : 


That is, when a unit r.m.s. voltage is applied to 


. . . 
the system, the power input is equal to the 
conjugate of the input admittance, as illustrated 
in Figs. 1, 2, and 4. 
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The method discussed in Parts I and II for the determination of steady-state solutions of 
electromagnetic problems is applied to a prolate conducting spheroid. An external e.m.f. is 


applied across the central section of the spheroid and the resultant field is represented as an 


infinite sum of wave modes expressed in terms of spheroidal functions. To each mode there 


corresponds a definite resonant frequency. As the applied frequency passes through the reso- 


nance point of a given mode, the corresponding term in the driving-point admittance changes 


from capacitive through a pure conductance to inductive. The sharpness of resonance increases 


with the eccentricity. Curves are drawn showing the behavior of the real and imaginary parts of 


the admittance for several values of eccentricity. The solution goes over, in the limit of zero 


eccentricity, to that discussed in II. When the eccentricity approaches unity one obtains a linear 


antenna of finite length and the driving-point impedance at the resonance frequency of the first 


mode is shown to be about 72 ohms. 


N Part I of this investigation’ we have dis- 

cussed the generation of waves in an infinite 
straight conductor, while Part IT? dealt with the 
forced oscillations of a sphere and the relation of 
radiation conductance and susceptance to the 
nature of the applied e.m.f. Although these 
examples give a clear insight into the general 
problem of forced oscillations, they are otherwise 
of no great technical importance. By the same 
method, however, we shall treat now the forced 
oscillations of a prolate spheroid, a problem of 
somewhat greater mathematical complexity but 
having a direct bearing on the operation of linear 
antennas, since a prolate spheroid reduces to a 
line of finite length as the eccentricity approaches 
unity. Furthermore, since the case of the sphere 
is obtained by allowing the eccentricity to be- 
come zero, one may follow the transition from the 
linear to a spherical configuration and the effect 
of increasing the central cross section on the 
radiation impedance. 

The free oscillations of a prolate spheroid were 
first studied in detail by Maclaurin® and by 


Abraham,‘ and much work has been reported 


tJ. A. Stratton and L. J. Chu, J. App. Phys. 12, 230 
1941). 

2]. A. Stratton and L. J. Chu, J. App. Phys. 12, 236 
1941). 

*R.C. Maclaurin, Camb. Phil. Trans. 17, 33 (1898). 
*M.. Abraham, Ann. d. Physik 66, 435 (1898); 2, 32 
1900). 
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since then on the orthogonal solutions of the 
wave equation in spheroidal coordinates. Several 
spheroidal functions have recently been investi- 
gated by Stratton® in an effort to establish rela- 
tions between the angular and radial functions 
defined as series of Gegenbauer and _ Bessel 
functions, respectively. This work has been ex- 
tended by the present authors and the results 
reduced to a form suitable for practical appli- 
cation.® The formulas obtained are the basis of 
the present paper. The free oscillations of a 
spheroid have lately been discussed by Page and 
Adams,’ who consider also a problem of forced 
oscillations. The driving field in their case is a 
homogeneous plane wave of length very much 
greater than the smaller dimension of the sphe- 
roid. A closer approach to the problem considered 
here has been made by Ryder,* based on the work 
of Page and Adams. The primary advantage of 
the functions which we shall employ in the 
following lies in the fact that they have been 
appropriately defined for every part of the 
domain of the independent variable and_ all 
analytic connections are specified. Convergence 
difficulties are thus largely circumvented and the 

J. A. Stratton, Proc. Nat. Acad. Sci. 21, 51—56 (1935); 
21, 317-321 (1935). 

®L. J. Chu and J. A. Stratton, J. Math. Phys. To be 
published. 


7 L. Page and N. I, Adams, Phys. Rev. 53, 819 (1938). 
*R. M. Ryder, thesis, Yale University (1940). 
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numerical discussion need not be confined to 


eccentricities near unity or zero. 

Let the foci of a family of confocal spheroids be 
located on the z axis at the points +f. A right- 
hand system of spheroidal coordinates £, n, @ is 
defined in terms of the cylindrical system r, ¢, 2 
by 

r=fl(t?—1)(1—7n’*) ]}', o=¢, z=fétn, (1 
where §21, —1=y=1, 0O=@=2rm. They are re- 
lated also to the spherical coordinates R, cos 0, ¢ 
by 


R= f(t?+n*—1)}, 
cos 6=én/(t2?+n?-—1)', d=. (2) 
As §— the spheroid goes over to a sphere such 


that 
R—ft, cos 0—-n. (3) 


Therefore & is often called the radial, and 7 the 
angular variable. If ig, i,, i, are unit vectors in the 
direction of £, n, ¢, respectively, then 

Véi=is/h:, Vn=i,/h,, Vo=1i6/hs, (4) 


and the metrical coefficients are 


ws 


he =fl(?—1)(1— 9?) Jj}. 


We shall consider a prolate spheroid of infinite 
conductivity embedded in a homogeneous perfect 
dielectric medium. It will be assumed that the 
applied electric field intensity-is circularly sym- 
metric about the major axis and linearly polarized 
in a direction parallel to this axis. All components 
of the field vary with the time as exp (tw). 
Then the Maxwell equations for the field in the 
external dielectric (€2, 42; o2=0) are 


1 0 
iweok: (h I] f 
hAedn 
—1 0 
lweok, (hall.), (6) 


hah: dé 
1 [o(h:E;) d(h,E,) 
tweet 


heh, On og 


Upon eliminating E; and E, one obtains for //, 
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the equation 


0? fl, Oo f/Tl, 
eR (te)sa-o 8") 
ae?\ h, On? \ h. 

0 fH, 0 fll. 
+4é ( —4yn ( ) 

at\h On\ h 


I] 
+ f*ko* £-— n°) G 7) 
h 


where ko= w(€o10)?. Let 
H.,=h, UE) Vine**. (8) 
I-quation (7) is separable and leads to 
(1—#2)U" —4£U' +(b,— frre?) U=0 (9) 


and an identical equation for V(m), where }; is one 
of a set of separation constants chosen such that 
I], remains finite at the poles n= +1. It appears 
that both U and V satisfy the general equation 


(1—2°)W”" —2(a+1)2W’'+(b-—c?s?)W=0 (10) 


for the special case a=1, c=fke. The angular 
function JV (») will be represented by the 
spheroidal function of the first kind Se;, ;(fke, 7) 
defined as 


1 ‘ P : ant 
Sei, (fk, n) =(n?—-1)7! DO! ant" 'T ysl), (11) 


where 

| P ° n! 1 

T n+2(n) = (9? —-1)? Pisiln), (12) 

(n+2)! 
. 

the P,.; are associated Legendre functions as 
defined by Hobson.® and the prime over the 
summation sign indicates that ” has the values 
n=(, 2, 4, 


? 


3, 5-++ if lis odd. The function 


if / is an even integer, and n=1, 


1 
Sey.1.— (9? — 1)? Piailn) 


as fke—0, i.e., as the spheroid degenerates into 
a sphere. At large distances from the center the 
held must reduce to a wave traveling radially) 
outward. To insure the proper behavior as > ~ , 
we choose for U() the radial spheroidal function 
of the fourth kind. 


Rey (fhe, § . 


=— : ; DY’ a hnsr( fhe), (13) 
fkot Do’ a,a"—' > 


- a m" 


Str 


’E. W. Hobson, Theory of Spherical and Ellipsoidal 
Harmonics (Cambridge University Press, 1931). 
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where /t,,; is the spherical Hankel function 
defined by 


o£) =(14 2fkot)!Hy.s o( fkeé), (14) 


and #1,,32 is the ordinary, or cylindrical, 
Hankel function of half-order. This construction 
. a ' af 

of the angular and radial functions Se; and 


4 » ae 
Re,,, employs the same set of coefficients a, and 
simplifies the computation accordingly. As 
fkot— « we find that 


1 


~— ‘ 
( fké)? 


to 
Paad 
t 


> 


4 : 
Re, i(fke, £) > 


while at the singular points &= +1 the function 
has poles of the first order. Thus the series 
- - . 4 . 
representation (13) of the function Re,,; diverges 
on and within a circle in the & plane of unit 
. . . i 

radius about the origin. However Re;,; may be 

4 “rr . - 
replaced by Se,,;, a different representation of 
one and the same function which converges in 
the neighborhood of £=+1. The function 

a" ‘ ‘ , , 

Se;.; is a linear combination of the two inde- 

. . 1 ‘ 2 . 2 . 
pendent solutions Se;,,; and Se;,;, where Se;,; is 

‘ : 1 . 7 
obtained by replacing the P,,; in Se;,; by the 
: 1 ; 

associated Legendre functions Q,,,; of the second 
eS l ,; ; 
kind.® But, Q,.; contains a factor that becomes 
infinite when 7 is a negative integer <—2. Con- 
sequently the series does not break off on the 
lower side at n=0 or n=1 as in the case of 


1 . ' a 
Se;,: the summation of the Se, series extends 


from negative infinity to positive infinity, over all 


even or all odd values of n as | is even or odd. 
The complete steady-state solution is, there- 


fore, 


Hs=>- Ah, U(E) V(ne'! 
= f{ (€2-—1)(1— 7?) Jie'! 


Xd A Re:*(E)Sei'(n), (15) 


1 


_  $ pe—2_ 
E;= : e* 
lwes \ £° — n° 


r d 
ye AiRei'(§)- [(1—7?)Se,"() ], (16) 
0 a7 
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: | 1—7* \3 
E. = - ete! 


1Weé» &- - n° 
a : 
X 2, ArSer'(n)- [(£2—1)Re,*(£) ], (17) 
0 d g 


where Re; i(fke, —€) and Se; i(fko, n) have been 
abbreviated to Re,*(é) and Se,'(n), respectively. 
The coefficients A; are to be determined by the 
distribution of applied e.m.f. 

We now apply to a perfectly conducting sphe- 
roid of eccentricity 1/£) an electric field E’ that is 
independent of ¢. The nature of this field need 
not be otherwise restricted. It may be the radi- 
ation field of a distant source whose wave-length 
is large relative to the dimensions of the spheroid. 
In such a case (15)-(17) above constitute the 
reflected wave.’ It may also be the field appearing 
near the terminals of a transmission line. In this 
case, a small section of the metallic spheroid is to 
be replaced by dielectric and the inner terminal 
faces of the metal are connected to the line as 
shown in Fig. 1. It is clear that the current from 
one terminal of the split spheroid will pass into 
one conductor of the transmission line, through 
the oscillator, and back to the other terminal of 
the spheroid. If the system is electrically balanced 
the current passing out of the one terminal must 
equal that entering the other, and the voltage or 
potential difference across the gap must equal 
that across the terminals of the line. The trans- 
mission line, therefore, serves the purpose of 
providing a metallic link to the split spheroid and 
applying a voltage across the gap. Neglecting 
distortion of the field caused by the presence of 
the line, we replace the line and source by a 
conductor short-circuiting the gap to insure 
continuity of current, and assume that a voltage 
equaling that across the terminals of the trans- 
mission line is applied over this short-circuiting 
link. The altered or equivalent system thus con- 
tains a complete spheroid with an external field 
intensity applied symmetrically over the surface 
of a short segment. The relation between current 
and voltage across this inserted segment is 
identical with that at the terminals of the line. 
This argument applies equally well to the case of 
a cylindrical conductor and that of a sphere, 
although in the latter case the cross-sectional 
dimensions are of the order of the wave-length 
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and the line must be designed such as to produce 
a circularly symmetric field. 

Since the conductivity has been assumed infi- 
nite the current within the spheroid must vanish. 
This implies that the sum of the applied field E’ 
and the induced field E is zero. Since the transi- 
tion of tangential components of E across a 
boundary is continuous, it follows that 


E,+E,' =0 at &=£o, (18) 
where E, is the field in the dielectric and given by 


(17). Let the applied field on the surface of the 
spheroid be expanded in the series 


1 1—7°* 
Ey = —- 


1weéo| Eo" — 7° 





> B,Se;'(n). 19) 


from the boundary condition (18) one obtains 


B 
A -, (20) 
(d/d&o)| (£9? — 1) Rex*(£o) | 
Togive a numerical illustration we shall assume 
a concentrated voltage V applied across an 
infinitesimal segment of length /A,Ayn at the 
equator »=0. Then 
E.’ = V/h,An when 9» <An/2, 


21 
E, =0 when » >An 


Nm 


From (19) and (20), and the orthogonal 
properties of the Se,'(m), we find 


B (wee V f)(Se(0) N a (2?) 
lweo V Se,'(0) 
1 (23 


f  Ni(d/d&o)[(Eo2—1) Re,*(Eo) ] 


where the normalizing factor is 


N 1 —n*)[ Se,'(m) |*dn 


Ae” 
22>’ 


(n+1)(n+2)(2n+3) 


B, and A, are zero for odd values of / since 


Se (0) =0 if 1 is odd. 


7 anemic, * SIE Fic. 1. Sphe 


roidal antenna 
transmission line 
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The field in the dielectric is now completely 
determined. At large distances from the spheroid 
the asymptotic expression for Re,;'(é) leads to 


IT g=(es uo) (ALS — (1 fEk2*)e wi—kesé 


Jt 


xX > 1'A,Se)"(n), (25) 


in which /=0, 2, 4--- 


center to the point of observation approaches fé. 


The distance from the 


The radiation pattern is determined by the sum 
of angular functions. 

It is of greater interest to learn how the power 
input of the system is distributed among the 
various “‘space harmonics,’ or modes, charac- 
terized by the index /, and the relation between 
the applied voltage and the current of the system. 
The surface current density J, on the spheroid is 
equal in magnitude to //, at the surface, or 
J,=Hg,\ ::,. The total current crossing a parallel 
of latitude on the spheroid in the y direction, i.e., 
from south to north, is 


J ,h.do= 2r fl &o7-—1)1- n*) }H, 
(26) 


2a f?(&o?—1)(1—7?) > A. Re;*(Eo)Se,'(n). 


Across the equator, where the voltage is applied, 
the current is 


tkof 
Inh=- V’ 
60 
_ (£9? — 1) Res*( Eo) [| Se:'(0) | 
xs —— amperes, (27) 
v (d déo)[ (Eo?—-1 ) Re,*( 9) IN, 


since 2mwe2= 27 €2/ 2) wl €2u2)'=k2/60 in ration- 
alized m.k.s. units. If one assumes the current 
uniform over the strip 4,An at the equator the 
input admittance is 


In —tkef 
V 60 
; &o°>—1)Re '(Eo)| Se '(Q)) 
x> ——mbho. (28) 
» (d d&q)| (£o7—1) Rer*( Eo) JN 
The reciprocal of this quantity is the input 
admittance Z;=1/ Y;. 
The expression for ); suggests at once, as in 
the case of the sphere, an equivalent system 
composed of an infinite number of circuits con- 
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Fic. 2. Components of input admittance as functions of L/X. 


nected in parallel to whose common terminals the 
voltage V is applied. Each component circuit 
represents one of the infinite sequence of wave 
modes. Consistent with the spherical case we 
shall define the admittance of the (/+1)th mode 


as 
tkof (&)2?—1)Re,*(£o){ Se,'(0) 
Vs 29) 
60 (d/dé&o){ (£o2?—1) Re/*( Eo) JN 


S 
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The total input admittance Y; is obtained by 
summing the ),,,; over even values of /. The real 
part of an admittance component Yj,; is the 
conductance g;.1, while the imaginary part is the 
susceptance b;,). 

In Fig. 2 the components of input admittance 
are plotted versus L/d for a series of values of the 
ratio L/D, where L is the length of the major 
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axis, D the length of the minor axis of the metal 
spheroid, and X is the wave-length measured in 
free space. The resonance point of the first mode 
falls within the range from 0.2 to 0.8 for L/X, the 
conductance g, of this mode reaching its peak 
value in the neighborhood of L/AX=0.5. In the 
same region the susceptance ); changes from 
positive to negative. The greater the departure of 
L/X from resonance, the smaller are the con- 
ductance and susceptance of the first mode. 
Within the same range the conductances of the 
higher modes are insignificant and the suscept- 
ances increase almost linearly with L/X. As the 
range of L/X is extended, the third, fifth and 
higher modes successively exhibit similar reso- 
nance phenomena in the vicinity of L/A1.5, 
2.5---, and then one by one drop out. 

A decrease in the eccentricity from unity to 
zero changes the shape of the spheroid from that 
of a thin wire of finite length to that of a sphere. 
2—1)}, 


where £) is the reciprocal of the eccentricity. At 


Eo /( 


Numerically the ratio L/D is equal to 


c 
s 


large eccentricities, say for L/D=7070 (Fig. 
2(a)), the conductance g; and susceptance }, of 
the first mode are sharply resonant. As the ratio 
L/D is diminished, the curves broaden at the 
peak while at the same time the susceptances of 
higher modes are augmented in magnitude. When 
L/D=7.05, (Fig. 2(d)), the conductance g,; tends 
asymptotically to a limiting value comparable in 
magnitude with the resonance value as L/A\ 
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Fic. 3. Input impedance and its real and imaginary parts 
as functions of LX. 


increases. This step-like increase in conductance 
is most striking in the limiting case of a sphere, 
L/D=1, discussed in Il. Summarizing briefly, it 
appears that for a thin wire each wave mode 
dominates only at or near its resonant frequency, 
while for a sphere the equivalent parallel branch 
circuit corresponding to each mode acts essentially 
as a Capacitance when the frequency is below the 
resonant value and as a resistance when the 
frequency is above this value. 

It should be remarked here that the present 
example is not wholly satisfactory, since the 
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series of input admittance terms fails to con 
verge. For large values of / the ratio of consecutive 
terms VY ,,:/¥,-, is always less than one but 
approaches unity as a limit as /— x. The same 
situation was encountered in the spherical case. 
This difficulty is no fault of the mathematical 
formulation of the problem but results from our 
assumption that the voltage is applied across a 
segment of vanishing length, implying in turn an 
infinite field intensity and infinite current density 
at the point of application. Actually there is 
always a finite separation between the terminals, 
and the length of this segment evidently bears an 
important relationship to the performance of the 
spheroid. A review of the properties of the 
angular spheroidal functions and the process of 
computing the A, shows that an increase in the 
length of the segment across which V is applied 
will result in a considerable decrease in the 
magnitude of the coefficients A, for large values 
of 1, and an inappreciable change for the lower 
modes. Thus the curves plotted in Fig. 2 are not 
far wrong if we ignore the contributions of higher 
modes. On this basis the real and imaginary 
parts and absolute magnitude of the input 
impedance have been plotted in Fig. 3. The real 
part, or radiation resistance, is zero at L/A=0 
and reaches the value 72 ohms when L/X is 
slightly less than one-half. This is in accord with 
the result obtained by calculating the radiation 
from the retarded potentials, assuming a sinu- 
soidal current distribution. The fact that this 
value occurs at a point slightly below L/A=0.5 
shows the effect of a finite thickness of wire and 
of a current distribution that deviates slightly 
from the sinusoidal. 

The radiation resistance attains a maximum 
somewhat below L/A=1, and passes through 
subsequent minima near L/A=1.5, 2.5, ete., the 
resonance points of the third, fifth and higher 
modes. The behavior of the reactance curves is 
illustrated by Fig. 3(b), showing that they pass 
through zero near resonance, and again at higher 
values of L/X, as one sees for the case L/D =7.05. 
The recurrence of zero reactance points is to be 
expected, for as the third mode starts to resonate 
at L/X=1.5, the reactance has to pass from 
negative value to positive. As shown in the curve, 
the spheroid of L/D=7 has a very small re- 


actance over a rather large range of L/X. For 
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large values of L/D there will be a sudden drop of 
reactance from large positive values to large 
negative values somewhere beyond the range of 
L/X shown in Fig. 3(b). The input impedance 
curves (Fig. 3(c)) indicate that sharp resonance is 
associated with thin spheroids, and that the 
impedance at resonance is almost independent of 
the ratio L/D unless L/D is near unity. 

Figure 4 shows the current distributions along 
a very thin spheroid. In this case the magnitudes 
of the third and higher modes are insignificant as 
demonstrated by Fig. 2(a) and 2(b). The abscissa 
n is directly proportional to the distance along the 
spheroid. At resonance L/A20.5 the distribution 
is sinusoidal, and the curve is sharper or flatter as 
L/ is less or larger than one-half. 

The power input into the spheroid, assuming 
that the voltage is applied across a very small 
distance, is one-half the product of the applied 
voltage and the conjugate of the current in the 
region where the voltage is applied. Thus 


kof 


P;=4VI,* =1—|V 
120 


s | Se(O) |7( £97 — 1) Re:*( Eo) 
Nid /dé, )| (£)7— 1) Re,*( Eo) | 


where P; is the input power, J»* indicates the 
conjugate of the quantity J) and Re/*(é ) is the 


radial spheroidal function of the third kind, 
which is the conjugate of Re;'(£)). Obviously 


P;=3)\V)\?Y;,*, (31) 


and needs no further comment. The real part of 
P; should equal the sum of the power dissipated 
in the metal and that radiated into space. There 


is no dissipation since the metal is assumed to 
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Fic. 4. Current distribution as function of 7 for 
thin spheroidal antenna. 





have infinite conductivity. The power radiated 
into space is equal to the integral of the Poynting 
the 
metallic spheroid. If we choose for this surface 


vector over a closed surface containing 
a confocal spheroid of eccentricity 1/£ then the 


radiated power P, is 


“de 


v/f Oo anf 


d nh hE 


Ro f P A | Se, (0) 
i Viz>’ 
120 \ 
& 1) Re,*(&)(d/dé)(E 1) Re,*(&) 
x $2) 


d/dé& )((é 1)Re;*(&o) |? 


It can be shown that the real part of P, is equal 





to that of P; by rationalizing the denominator 
of P; and evaluating the Wronskian of 


0 
(£?—1) Re,*(&)- -L(€*— 1) Re,*(&) |. 
O§ 


Particularly at the surface of the conductor, 


< 


&= ££) and 


P = /P.. 33) 


Thus we establish the theorem that the complex 
power input into the spheroid is equal to the 
complex radiated power at the surface of the 
spheroid. 

We shall discuss later in another place the 
relation of these results to engineering formulas 


in current use for practical antenna design. 
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Much experimental work has been done on the creep of metals and the long-time creep test is 


the standard means of rating resistance to slow deformation of metals under tensile stress. In 


engineering practice, frequently, cases are found in which the total strain and not the stress 


remains constant. In what has become known as relaxation, an initially pure elastic strain 


gradually is transformed into a permanent strain while the length of the bar remains unchanged 


Several years ago some new apparatus was designed and built primarily to duplicate the above 


conditions in connection with the study of 


bolt 


materials for steam turbines. The unusual 


wccuracy and exceptionally good test results obtained from the automatic relaxation machine 


during more than two years of operation indicate that the type of test may well become a means 


of comparison as standard as a creep test and at a great reduction of testing time. 


I. INTRODUCTION 


NDUSTRY has become considerably inter- 

ested, during the last fifteen years, in what 
happens to metals under stress at high tempera- 
lhe the i 


efhciency of steam turbines are partially set by 
the the 


tures. effective limits to increase in 


temperature at which steels may 


stressed without their “creeping” bevond 
allowable amount for a period of years. Creep is 
the term used to describe the slow plastic de- 
formation taking place in a material due to 
applied stress. This phenomenon is hardly notice- 
able in steels at room temperatures but becomes 
of serious magnitude in most of even the best 


alloy steels at the maximum operating tempera- 
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be 
the 


tures of steam turbines (near 500°C). The creep 
test, which measures the deflection at constant 
stress, has become the standard method of 
rating resistance to slow deformations of metals 
under tensile stress. It is, however, a long-time 
test because the time normally required for 
accomplishing one is from three to six months. 
In high temperature apparatus, frequently, 
parts are not subjected to constant load but to 
total the 


portant example of this. Attempts to use creep 


constant strain: bolts are most im- 
data in estimating the relaxation which occurs 
unsuccessful. The relaxation 
built 


urgently needed information. The machine has 


have so far been 


machine was _ therefore to supply this 
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proved itself more than successful since it filled 





the direct need for design purposes, and now it 
can be shown that the results have further 
significance. The relaxation test is a short-time 
test since it normally requires only one week. 
Krom the large number of test results now 
available we conclude’ that one can obtain 
qualitative information in terms of creep—and 
possibly even quantitative information. More 
tests will be required before we can ascertain 
how accurately it will be possible to estimate the 
creep curve. 

The relaxation phenomenon is exhibited when 
a material is stretched to, and subsequently held 
at, a given elongation. Creep then acts to 
elongate the bar and the release of load causing 
elastic shortening balances this effect. The total 
strain remains constant. 


e+e’ =ce=0,/E (1 
e’=0/ E (Hooke’s law) 

ldo de’ 

Edt dt | 


+ 
\| 
NR 


é’’=plastic strain, ¢«’=elastic strain, e=total 
strain, o=stress, ¢;=initial stress, = modulus 
of elasticity. The plastic elongation can be 


measured by the decrease in stress which is of a 











magnitude that can conveniently be observed. 





In our own laboratory we previously experi- ; 
; ; : ° Fic. 2. Automatic relaxation machine. 
mented with a smaller relaxation machine which 


O’= 40,200 LB./ IN2 | 
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Fic. 1. Diagrammatic sketch of automatic 
relaxation machine. biG. 3. Stress-strain diagram of relaxation test 
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Fic. 4 
gave good results for copper at low tempera- 
tures.' In recent years relaxation tests have been 
used by W. Barr and W. E. Bardgett,? N. L. 
Mochel,’ E. L. Robinson,‘ and C. C. Davenport, 
but the number of test results available is com- 
paratively small. Our present automatic relaxa- 
tion machine was developed and built at th 


Westinghouse Research Laboratories.® 


: | Boyd, Proc ASI 
?W. Barr and W. E. 

122, 285 (1932) 
N. I ° Mochel, lrans 


M. 37, Part II (1937) 
Bardgett, Proc. Inst. Mech. Eng. 


A.S.M.E. 59, 453 (1937). 


‘E. L. Robinson, Trans. A.S.M.E. 59, 451 (1937) 
C. C. Davenport, J. App. Mech. 5, 55 (1938). 
*A. Nadai and J. Boyd, Proc. Fifth International 


Congress of Applied Mechanics, Cambridge, Massa 
chusetts, 1938, p. 245. Special acknowledgments are due to 
Mr. A. D. Irwin, for design of the machine. 
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Relaxation record of steel at 500°C. 


Il. RELAXATION MACHINE 


To duplicate relaxation conditions, load must 
that the 


shortening balances the permanent elongation. 


be continuously removed so elastic 
This ideal type of test is approached by having 
very many small steps of unloading. The required 
elements of the machine, therefore, are a spring 
loading mechanism and an accurate but sturdy 
extensometer. A diagram of the machine is 
1. See, also, Fig. 2. 


specimen (/7) is loaded through a 


shown in Fig. 

The tensile 
lever arm (B) by springs (C) whose extension is 
D). The 
springs are loaded and unloaded by a 
(Ji. The 


measured by a follow-up mechanism 
worm 


gear drive motor driven follow-up 
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device is directly coupled with the recorder so 
that the stress, which is directly proportional to 
the spring elongation, can simply be calculated 
from the chart. 

The (K) is attached to the 
specimen at the gage length of 6 inches. The 


extensometer 


steel extensometer arms extend down through 
the bottom of the furnace and at (G) the small 
relative motion is mechanically amplified. An 
electric contact with a micrometer measuring 
device accurately determines the total strain. 
After the specimen is loaded this contact operates 
the unloading motion for a strain of approxi- 
mately 2 10~-* inch per inch. 

The electric furnace (A) heats the specimen 
and controls the temperature at any desired 
level below 750°C. The differential expansion of a 
quartz rod and a nickel tube operating another 
contact device is utilized to interrupt part of the 
heating current. In this way the specimen tem- 
perature varies no more than +1°C. 

On all of the control contacts a 6-volt potential 
is applied which, when contacts are closed, is 
amplified by a thyratron tube to operate a relay. 
This small voltage eliminates sparking and there- 
fore increases accuracy and decreases wear. These 
automatic devices have given no trouble and 
require attention only several times a year. 

Figure 3 shows a loading curve and indicates 
that the small strains required to operate a 
loading motor are less than 1/100 of the total 
strain. 

Figure 4 is a picture of a continuous 1100-hr. 


load-time record of a relaxation test on an alloy 
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Fic. 5. Relaxation of Cr-Mo-W steel at 500°C. 
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steel at 500°C. This unusually long test was used 
to check both the reliability of the machine and 
the accuracy of the extrapolation. At points B 
and C the pen was shifted so that the line would 
the paper. At B the 


reduced to } the initial scale because the stress 


stay on time scale was 
was changing more slowly. At E a slight dis- 
turbance can be noticed; this was caused by a 
change in room temperature, which could not 
always be controlled. The disturbance was of no 
consequence since the stress was reduced less 
than 1 percent and for only a short time. For 
perfect test results air temperature and circula- 
tion must not be disturbed near the machine. 
III. ResuLts 

Since one of the main purposes of the con- 
machine 


struction of the automatic relaxation 





was to study mechanical properties of high 
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kic. 7. Relaxation of Cr-Mo-W steel (2) at 500° 
Remanent stress vs. initial stress. 
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Fic. 8. Relaxation values at 1000 hr. Initial stress 
a, = 40,000 Ib. /in2 


temperature bolting materials, a great amount of 
data has been gathered on this subject. It was 
found that, usually, after one week the stress 
decreased very slowly and so the usual time of 
testing became about 150 hours. 

Figures 5 and 6 show typical plots of test data. 
It is convenient to use the semilog plot because 
the logarithmic time scale fits both the rapid 
initial the subsequent 
changes. Curve FE of Fig. 5 has an inflection 
point at approximately one-half of the initial 
stress. This observation and the assumption 
that eventually the 
asymptotic to the zero stress axis have been the 
guide for extrapolation of the relaxation data. 
Figure 6 is one of the extrapolated series of 
curves in which curve C is a replot of the record 
shown on Fig. 4. The original curves of this 
material were drawn and extrapolated before the 
check test was run. The check test not only 


relaxation and slow 


relaxation curve becomes 


agrees well with the original test whose initial 
stress was 40,000 Ib. sq. in. but it also coincides 
with the extrapolation. In most cases when check 
tests were run starting with the same initial 
stress there was agreement between the observa- 
tions. When more than one relaxation machine 
becomes available, other long tests will be run 


to check even further extrapolations; however 


those to 10,000 hours are probably fairly ac- 
curate. 

A cross plot of Fig. 6, taking time as the 
parameter, is shown in Fig. 7. This is a useful 
plot for the engineer as it simplifies interpolation. 
It happens that a parabolic function o,—¢,x 


=¢(or)* (where o;=the initial and op =the final 
stress) expresses the curves closely for all of the 


high strength materials. The parabolas are drawn 
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through the origin O making them tangent to. 
line OA (Fig. 7). Theoretical considerations of 
possible creep laws show, however, that this 
assumption may not be true, for the indication 
is that straight lines fan out from the origin O at 
the low stresses. For test data on very weak 
materials it was found that the parabola does 
not fit. Although we realize that this is an em- 
pirical statement and we do not know, therefore, 
the more exact shapes of the curves of Fig. 7, 
we do know that for strong materials such curves 
seem probable from the observations. 

In order to compare the relaxation of steels at 
various temperatures, a comparison was made of 
the relaxed values at an initial stress of 40,000 
lb./sq. in. and at a time of 1000 hr. Figure 8 
shows readily which is the strongest material 
because of its resistance to creep. Materials may 
well be compared by means of a “‘creep critical 
temperature range’’ which would be the range of 
temperature near the point of maximum slope 
of this type of curve. In creep tests there is also 
a range of temperature at which the creep rate 
for a given stress increases very rapidly, but 
since there are no creep test data at various tem- 
peratures on the above materials it is impossible 
to show that are 


identical. The significance of this critical tem- 


these temperature ranges 
perature can be appreciated in that it is known 
from several other types of tests that the alloy 
K-42-B has excellent high temperature proper- 
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Fic. 9, Relaxation of K-42-B at 1000°F. 
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ties, and a glance at the critical temperature 
comparison brings out the order of increasing 
creep resistance of the steels which were com- 
pared. 

Since relaxation 


under certain conditions is 


essentially an integrated creep phenomenon 


under variable stress, it would seem that we 
should be able to correlate the two types of test 
results. Unfortunately, no complete series of both 
tvpes of test results on the same material are 
available, but an attempt at least to correlate 
tests which were available for K-42-B 
at 538°C (1000°F) will be described. Figure 9 is 
the 


in., and the parabolic 


several 
the plot of the relaxation test at initial 
stress of 40,900 Ib. ‘sq. 
curves used to interpolate the information re- 
quired. In Fig. 10 the interpolated relaxation 
that a 20,000 


is covered. The assumptions for the 


curves are drawn so range at 
Ib. ‘sq. in. 
correlation are that, given the same material at 
zero time, at equal times and stresses the plastic 
strain rates are also equal. As mentioned before 
(see Eq. (2)), the permanent strain e’’ is a func- 
tion of the change of stress in the relaxation 
curve, and a creep curve ¢«’’=f(t) at a constant 
stress can be plotted using an equivalent stress 
o= Ke’ as the coordinate instead of the strain. 
A graphical integration is made of the relaxation 
curve slopes at the stress of 20,000 Ib. sq. in. 
and a creep curve is calculated. The slopes at 
points A to F are redrawn near the creep curve 
to show that this is the basis of the correlation. 
Both calculated and tested curves are plotted on 
the coordinates and_ their 


normal strain-time 


close correlation can be noted. Two other 
attempts at this type of correlation were equally 
successful. Although in no case was there suffi- 
cient information, all evidence so far available 
points in the direction that this correlation can 
give good results. 


The relaxation test is useful beyond the direct 
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Fic. 10. Correlation diagrams. 


application of information regarding the behavior 
of metals at constant total strains. It establishes 
the relative ability of a material to withstand 
stress at high temperatures, particularly in the 
applications of bolting materials. It is also 
valuable in that it may be used to estimate creep 
values in certain industrially important ranges of 
stress and temperature for high strength alloys. 
K:ffectively, this is a much needed short time test 
to predict mechanical properties of metals at 


high temperatures. 
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Engineering Defense Training Program 


Che following is a list of the engineering defense training 
courses which have been authorized to December 30, 1940. 
Courses are all of college grade. The cost of instruction 
will be paid by the Federal Government, but students are 
required to provide their own textbooks and supplies. 


Alabama 
Alabama Polytechnic Institute, Auburn, Alabama: 
Engineering, (2) 
4) Industrial Relations, 
7) Military 
duction Engineering, 


1) Construction 
Engineering Economy, 
5) Machine Design, (6) Marine Engineering, 


Naval Architecture, (9) Pro- 


Engineering Drawing, (3) 
Highway Engineering, (8) 
10) Surveying, (11) Tool Engineering 
University of Alabama, University, Alabama: (1) 
tor Ordnance Inspectors, (2 
Time and Motion Study, (3 


Correlation Course 
Industrial Organization and Management, 

Industrial Practice, ¢ 
Practice, Electric 


ast Iron Pipe and 


Fittings Industry, (4) Industrial Power Industry, 


5) Industrial Practice, Steel Industry, (6) Industrial Safety and Safety 
Materials, (8) Mecha- 


9) Physical Metallurgy and 


\dministration, (7) Mechanics and Strength of 


nisms and Elementary Machine Design, 
Metallurgy of Iron and Steel, (10) Tool and Die Design 
Arkansas 


University of Arkansas, Fayetteville, Arkansas: (1) 


Drawing 


Engineering 


California 
California Institute of Technology, Pasadena, California: (1) Pro 


duction Engineering 2) Production Supervision, (3) Radio and 
Electronics. 
University of Southern California, Los Angeles, California: (1) Ma- 


terials Inspection and Testing (Unit A), 
Unit B), (3) 


2) Materials Inspection and 
Testing Materials Inspection and Testing (Unit ¢ 


Colorado 
Colorado School of Mines, Golden, Colorado: (1) Production 
Major Nonferrous Metals, (2) Theory and Application of Welding 
Colorado State College of Agricultural and Mechanical Arts, Fort 
Collins, Colorado: (1 Machine 
(3) Materials Inspection and Testing 
University of Colorado, Boulder, Colorado: (1) 
2) Airplane Stress Analysis, (3 
Machine 
Production Engineering 


Engineering Drawing, (2) Design, 
Advanced Electroni 
Instruments, Engineering Drawing and 
Material Testing and In- 


Technology of Heavy Chem- 


Machine Design, (4 Design, (5 
spection, (6) 
icals Industry, (8) Theory of Elasticity 

University of Denver, Denver, Colorado: (1 


and Other Defense Materials 


Testing of Chemicals 


Connecticul 


Yale University, New Haven, Connecticut: (1) Applied Mechanical 


Design, (2) Engineering Drawing, (3) Inspection of Materials, (4) Met- 
allurgy, Production Control, (6) Welding Supervision 
Delaware 

University of Delaware, Newark, Delaware: (1) Motion and Time 

\nalysis 
District of Columbia 

George Washington University, Washington, D. C.: (1) Concret 
Inspection and Testing, (2) Construction Materials, Properties, Tests 
and Specifications, (3) Map Making and Interpretation, (4) Radio 
Communication, (5) Ship Construction, (6) Structural Steel Design, 


7) Surface and Aerial Navigation 


Howard University, Washington, D. C.: (1) Engineering Drawing 


Florida 


University of Florida, Gainesville, Florida: (1) Engineering Drawing 


Georgia 
Georgia School of Technology, Atlanta, Georgia: (1 
Materials, (2 Metal Inspection, (4 


( onstruction 
Engineering Drawing, (3 Textile 


Testing 
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Idaho 
University of Idahe, Moscow, Idaho: (1 
neering Drawing. 


Design Drafting, (2) Engi- 
Illinois 
Bradley Polytechnic Institute, 
Drawing. 
Illinois Institute of Technology, Chicago, Illinois: (1) 
Tools 
Drafting and Elementary Design, (5) 


Peoria, Illinois: (1 Engineering 

Design of 
Shelters, (2 Equipment, 
Industrial 
Machine Design, (8) Metal 
10) Strength of Materials and 
Structural Designing and Drafting, (12 
Testing Methods, (14 
Tool and Fixture Design, (16 


Bombproot Design of and Testing 
3) Diesel Engines, (4 
Management, (6) Inspection Methods, (7) 
9) Production Planning Routing, 


Steel Design, (11 


lurgy, 
Testing and 


Inspection of Concrete, (13 Time and Motion 


Study, (15 Welding Engineering 


Indiana 


Purdue University, Lafayette, Indiana: (1) Aeronautical Engineering, 


2) Applications of Electronics, (3) Engineering Drafting, (4) Explosives 


Inspection, (5) Machine Design, (6) Materials Inspection, (7) Materials 


Testing and Inspection (2 courses), (8) Metallurgy, (9) Production 
Engineering courses), (10) Production Supervision (2 courses 
11) Tool Design (2 courses). 

Rose Polytechnic Institute, Terre Haute, Indiana: (1) Electric 
Communication, (2) Engineering Drawing and Machine Design, 
3) Illuminating Engineering, (4) Materials Inspection and Testing, 
5) Production Engineering. 

University of Notre Dame, South Bend, Indiana: (1) Chemical 


Materials, (2) Production 
Tool and Die Design 


Analysis of Metallurgical 
3) Physical Metallurgy, (4 


Engineering, 


lowa 
lowa State College, Ames, Iowa: (1) Materials Inspection and 1 
ing, (2) Tool Engineering 
Kansas 
Kansas State College, Manhattan, Kansas: (1) Aeronautical Engi- 
neering, (2) Engineering Drawing, 3) Explosives, } Materials 


5) Tool Engineering. 
University of Kansas, Lawrence, Kansas: (1) Aeronautical Engi 
Machine Design, (4) Materials 


Production Engineering and Supervision. 


Inspection and Testing, 


neering, (2) Engineering Drawing, (3 


Inspection and Testing, (5 


Kentuc ky 
University of Louisville, Louisville, Kentucky: (1) Aerial Photo 
graphic Mapping, (2) Applied Mechanics (Statics), (3) Elementary 
Unified Mathematics, (4) Ferrous Physical Metallurgy 5) Non 
ferrous Physical Metallurgy, (6) Strength of Materials 


Loutstana 
Louisiana Polytechnic Institute, Ruston, Louisiana: (| 
Drawing 
Southwestern Louisiana Institute, Lafayette, Louisiana: (1) Engi 
neering Drawing 


Engineering 


Tulane University, New Orleans, Louisiana: (1) Engineering Drawing 
? courses 2) Machine Design 
Maryland 
Johns Hopkins University, Baltimore, Maryland: (1) Electricity 
Applied to Aircraft, (2) Industrial Organization and Management 
3) Materials Inspection and Testing, (4) Mechanical Design for Radi 
Engineers, (5) Plastics, (6) Production Fundamentals Rad 


8) Time and Motion Study 

University of Maryland, College Park, Maryland: (1) Acronautical 
\ircrait Tool 
Radio Engineering, (5) Radio Testing and Inspecti 


Engineering, 


Dratting and Design, (2) Inspection, (3) Aircraft Engi 


neering, (4 


Massachusetts 
Massachusetts Institute of Technology, Cambridge, Massachusetts: 
1) Aeronautical Engineering, (2) Aircraft Engines, (3) 
Applied 
Coordination and Expediting, (6 


Applications of 


Metallography, (4 Mathematics, (5) Defense Production 
ome ha 
Exterior Ballistics, (9) In 
Marine 


Inspection 13) 


Dynamic Analysis of Serv 


nisms, (7) Engineering Fundamentals, (8 
strumentation and Vibration Measurement, (10 
11) Naval Architecture, (12) Ordnance 


Testing, (14) Vibrations and Their 


Engineering, 
Textile 
Applications. 
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Northeastern University, Boston, Massachusetts: 1 Applied 
Mechanics, (2) Elementary Airplane Structures, (3) Engineering 
Drawing, (4) Higher Mathematics, (5) Instrumentation, (6) Machine 


Design (2 courses), (7) Materials Inspection and Testing (2 courses 


8) Methods Engineering, (9) Production Engineering and Supervisiotr 
19) Production Supervision 

Tufts College, Medford, Massachusetts: (1 
2) Chemistry in National Defense, (3) 
Construction, (4) Engineering Drawing, (5 
Instructor Training, (6 


Applied Mathematics, 
Concrete Inspection and 
Foreman Conferences and 


Metallurgy, 
Theory and 


Machine Design, (7 
8) Production Engineering, (9 
Practice of Welding. 
Worcester Polytechnic Institute, Worcester, Massachusetts: (1) Ele 
ments of Machine Design, (2) Elements of Tool Engineering, (3) Ma 


Physical 
Radio Engineering, (10 


terials Testing and Inspection. 


Michigan 

Detroit Institute of Technology, Detroit, Michigan: (1) 
Design, (2) Tool and Die Design. 

Michigan College of Mining and Technology, Houghton, Michigan: 
1) Engineering Drawing, 
and Testing. 

University of Detroit, Detroit, Michigan: (1) Aeronautical 
neering, (2) Drawing, Descriptive Geometry, and Shop 
3) Materials Inspection and Testing, (4) 
5) Production Supervision. 

University of Michigan, Ann Arbor, Michigan: (1 
2) Materials Inspection and Testing, (3) 
Supervision, (4) Tool Engineering. 

Wayne University, Detroit, Michigan: (1 
2) Engineering Drawing I, (3) Engineering Drawing II 
C,eometry), 


Machine 


2) Machine Design, (3) Materials Inspection 
Engi- 
lathematics, 
Production Engineering 
Machine Design, 
Production Engineering and 


Electronics Engineering I 
Descriptive 
4) Engineering Drawing III (Freehand Drawing), (5) In- 
dustrial Organization and Factory Management, (6 
7) Material Testing, (8) 


and Measurements, (10 


Machine Design, 
Metallurgy, (9) Shop Processes, Mathematics 


Telephone Communication 


Minnesota 
University of Minnescta, Minneapolis, Minnesota: (1 
Machine Design, (3 
+) Optical Engineering 


Engineering 


Drawing, (2 Metallurgy and Metallography, 


Mississippi 
Mississippi State College, State College, Mississippi: (1 
Mater 


Engineering 
Drawing, (2 ials Inspection and Testing, (3) Tool Engineering 
Missouri 
Washington University, St. Louis, Missouri: (1 
2) Material Production 


4) Production Supervision, (5) Tool Engineering 


Machine Design, 


Inspection and Testing, (3 Engineering 


Nevada 


University of Nevada, Reno, Nevada: (1) Civil Engineering Drawing 


New Hampshire 


University cf New Hampshire, Durham, New Hampshire: (1) Engi- 
neering Drawing 


New Jersey 


Newark College of Engineering, Newark, New Jersey: (1 
Machine Design 3 


ing 


neering Drawing, (2 Materials Inspection and 


Testing 

Rutgers University, New Brunswick, New Jersey: (1) Machine 
Design 

Stevens Institute of Technology, Hoboken, New Jersey: (1) Elements 


of Engineering Drawing, (2) Introduction to Engineering 


New Mexico 
New Mexico College of Agricultural and Mechanical Arts, State 
College, New Mexico: (1) Engineering Drawing, (2) Machine Design 
3) Materials Inspection and Testing. 
New York 
Cornell University, Ithaca, New York: (1) Advanced Stress Analysis 


and Elastic Stability, (2) Electronics, (3 
Metallurgy, (4 


Materials Testing and Physical 


Mechanics and Elementary Aircraft Design. 
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New York University, New York, New York: (1 
ment Protection, (2) Airc: 
Advanced Aircraft 


\erial Bombard 


ift Stress Analysis Procedure, (3) Elements o 


Engine Design 4) Food Inspection, (5) Funda 
mentals of Aircraft Fitting Design, (6) Gaging and Inspection Methods, 
Heat Treatment of Aluminum Alloys, (8) Introductory Aerody 


namics, (9) Marine Engineering, (10) Materials Testing with Labo- 
ratory, (11) Metallurgy and Metallography, 
13) Production Control, (14) Time and Operation Study. 

Pratt Institute, Brooklyn, New York: (1 


12) Military Sanitation, 


Engineering Drafting and 


Practice, (2) Production Engineering 

Union College, Schenectady, New York: (1) Chemistry of Metal 
2) Design of Highly Stressed Structures and Apparatus, (3) Drafting 
ind Drafting Room Practice, (4) Electron Tubes and Their Applica- 


tions, (5) Elementary Electrical Theory and Practice, (6 
Mechanics and Strength of Materials, (7 


Manufacturing Methods, (8) 


Elementary 
Engineering Materials and 
Engineering Problems and Their Mathe- 
9) Industrial Applications of Optics, (10) Industrial 
11) Internal Combustion Engines, (12 


matical Solution, 
Chemistry, Methods Improve- 
ment by Motion Study, (13) Principles and Practice of Radio Communi- 
cation, (14) Production and Cost Control, (15) Tool Design 


North Carolina 
Agricultu:al and Technical College of North Carolina, Greensboro, 
North Carolina: (1) Materials Inspection and Testing. 
Duke University, Durham, North Carolina: 
Engineering, (2) Machine Design. 
North Carolina State College, Raleigh, North Carolina: (1 
Inspection, (2) Chemical 


1) Communication 


Aircratt 
Diesel 
Electronics, 
Machine Design, (8) Ma- 
Production Engineering, (10) Pr 


Testing and Inspection, (3) Engi 


neering, 4) Engineering Drawing, (5) Experimental 
6) Instrument Men and Topographers, (7 
terials Inspection and Testing, (9 


duction Supervision 


North Dakoia 


University of North Dakota, University, North Dakota: (1) Engi 
neering Drawing, (2) Materials Inspection and Testing 
Ohio 
Case School of Applied Science, Cleveland, Ohio: (1) Electrical 
Maintenance and Control, (2) Elements of Motion Economy, (3) Engi- 
neering Surveying, (4) Explosives, (5) Fundamentals of Machine 


Design, (6) Industrial Instrumentation, (7) Inspector and Metallurgica! 


Technician Training, (8) Internal Combustion Engine Testing, (9) Ma 
chine Detailing. 
University of Cincinnati, Cincinnati, Ohio: (1 


and Testing, 


Materials Inspection 


2) Production Supervisor. 


Fenn College, Cleveland, Ohio: (1) Engineering Drawing, (2) Ma 
terial Inspection and Testing, (3) Nonmetal Materia! Testing, (4) Pro 
duction Engineering, (5) Production Supervision, (6) Stress Analysis 


7) Structural Detailing, (8) Tool Engineering 
Ohio State University, Columbus, Ohio: (1 
munications 
Machinery 


6) Inspection of Materials of Chemical Industry 


Air Navigation, (2) Com 
Plant 
Explosive s 
A.S. T. M. methods 


Inspection of Construction Materials, (8 


Laboratory Technique, 3) Design of Chemical 


and Equipment, (4) Engineering Drawing, (5 


and Chemicals, (7 Inspection 


of Instruments, (9) Inspection and Testing of Small Electrical Powe: 
Units, (10) Inspection and Testing of Small Internal Combustion 
Engine Power Units, (11) Inspection of Steel and Cast Iron Products 
12) Inspection of Welded Products, (13) Machine Design, (14) Physica 
Metallurgy, (15) Production Engineering, (16) Production Supervision, 


17) Tool Engineering, (18) Welding Engineering Design. 


Oklahoma 


Oklahoma Agricultural and Mechanical College, Stillwater, Okla- 
homa: (1 Machine Design, (3) Materials 


Preduction Engineering, (5) Tool Engi- 


Engineering Drawing, (2 
Inspection and Testing, (4 
neering. 

University of Oklahoma, 
Engineering, (2 


\eronautical 


Oklahoma: (1 


Engineering Drawing (2 courses). 


Norman, 


Oregon 


Oregon State College, Corvallis, Oregon: (1 
neering 


\eronautical Eng 








Penns ylvanta Texas 


Bucknell University, Lewisburg, Pennsylvania: (1) Engineering Agricultural and Mechanical College of Texas, College Station, 
Drawing, (2) Materials Inspection and Testing, 3 Productiot Texas: (1) Aeronautical Engineering, (2) Aircraft Inspection, (3) Camp 
Engineering. Sanitation, (4) Engineering Drawing, (5) Material Inspection and 
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N a paper by Raymond B. Block of the same title 
appearing in this journal,'! experiments on rolling a 
steel ball filled 


liquid were reported and it was found that the data for 
different liquids did not fit either Stokes’ or Oseen’s law 


down an inclined tube with a viscous 


for the viscous drag on a sphere. The usual assumption 
made in fluid mechanics is that there is no relative motion 
of the liquid at the surface of the sphere, but he concluded 
that if tangential slipping is postulated, the data can be 
used to calculate a coefficient of slip which varies froin 
liquid to liquid. It is suggested following Brillié? that this 
slip should be greater with polar compounds. 

The hypothesis of slip has intrigued investigators since 
the time of Helmholtz, and many experiments have been 
conducted in order to prove its existence, with no positive 
results.4 Hence, if the conclusions of Block’s paper are 
justified they have very far-reaching significance, among 
other things the fact that any measurement of viscosity 
would always depend on the apparatus and the liquid 
under investigation. 

The criticisms of the work under discussion fall into 
1) whether Stokes’ 
law should apply in the present experiment even if there 


three general categories: or Oseen’s 


were no slip, (2) whether the effect of slip is also present 
in the measurement of viscosity and should be allowed for 
in reducing the data, and (3) whether glycerin itself is a 
polar compound. 

Considering these objections in turn, the first relates to 
the applicability of Stokes’ or Oseen’s law in the present 
case. Both these equations, as the author points out, 
apply to the uniform translation of a sphere through an 
infinite amount of liquid and not to combined translation 
and rotation in a relatively small tube. The effect of the 
walls of a cylindrical tube, assuming the sphere sym- 
metrically placed in the tube has been considered by 
Ladenburg,® by Faxen,’ and by Francis.* For the apparatus 
used by Block the correction factor for the walls will be 
quite large and in that case will depend on the speed and 
the viscosity. Hence, it should be taken into account when 
The 


imposed rotation of the sphere complicates the problem 


comparing different liquids. matter of the super- 
and no rigorous treatment of it has vet been given. It 
would have to be shown that a linear relation did exist 
between the drag, F, and the quantity, wed, where yu is the 
absolute viscosity, d the diameter of the sphere and v its 
velocity. Even assuming for the moment that this is so, 
a more serious objection lies in the fact that both these 
equations that the 
“‘Reynolds’ number”’ pvd/u, where p is the density of the 
liquid, must be much less than unity. 


bee ‘n 


were derived on the assumption 


Many experiments® 


have performed to test this relationship. These 
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include dropping Rose metal spheres in several vegetable 
oils and steel spheres in water and glycerol. A summary" 
of the results shows that Stokes’ law is precisely obeved if 
the Reynolds’ number is less than 0.5 and that Oseen’s 
law is obeyed for values of the Reynolds’ number up to 
0.8. Above these points the experimental values begin to 
deviate from the theoretical curves as would be expected. 
Block shows three plots in which the variable IH», which is 
the Reynolds’ number as defined here, appears. For the 
experiments using pure glycerin Il. ranges from 0.012 to 
0.25, for the solution containing 10 percent methyl alcohol 
in glycerin it ranges from 0.63 to 2.0, and for the other 
solutions used IT, is even higher. It can be seen that except 
in the case of pure glycerin neither of the theoretical 
equations are even supposed to apply to the present 
experiment, and hence any deviations therefrom cannot 
be taken as evidence of slip. The true explanation probably 
lies in the partial or complete neglect of the inertial terms 
in the equation of motion.5 

The second criticism of the conclusions drawn in the 
paper under discussion is that in any investigation of the 
existence of slip by this method account should be taken 
of its effect on the measurement of viscosity by means of 
the modified Ostwald capillary viscosimeter. If we assume 
a coefficient of slip, 8, then Poiseuille’s equation for the 
discharge rate, Q, of a liquid through a tube of radius r 
and length / under a difference of pressure p, becomes: 


O=(nr'ip 8u)(1+4u/8r). 


Taking the values of wu and 8 given in Table IV of 
Block’s paper and assuming a radius of the capillary of the 
viscosimeter of 0.2 (1+4u/6r) 
approximately equal to 4. (The end effects will modify 


this result only slightly.) This means that a liquid which 


cm, the factor becomes 


slips on the surface will flow through the viscosimetet 
four times as fast as a similar liquid of the same absolute 
viscosity which does not show slip. This result seems highly 
unlikely, but if the object of the present work is to prove 
the existence of this phenomenon then either this factor 
should be taken into account in the formulae or else it 
should be definitely established that none of the liquids 
used exhibit slip on the glass walls of the viscosimeter. 
In this latter connection Whetham'® found no observable 
deviations from Poiseuille’s equation for the flow of water 
through either silvered or unsilvered glass tubes. 

Lastly, it has been well established that glycerol is a 
polar compound!*~!® and the value of its dipole moment!” 
has been given as 2.20 X 10~'§ as compared with 1.68 x 107" 
for methyl alcohol and zero for a nonpolar substance. It 
appears then that one would expect both substances to 
behave alike even on the hypothesis suggested by Block 


i) 
uw 
~ 





undoubtedly 


and that the differences observed by him are 


due to the great disparity in the viscosities with its re- 
sultant effect on the inertial terms in the equations of 
motion 

It should also be pointed out that Arnold® and Allen"™ 
were able to verify Stokes’ law for the low range of 
Reynolds’ number using colza oil, linseed oil, aqueous 
alcohol and water, all of which are or contain polar 
substances. 

If there is still any question as to the non-existence of slip 
in spite of the many experiments that have been performed 
to find it, a few of which have been cited above, then the 
critical experiment should be «somewhat as follows: the 
sphere should be allowed to fall freely in a cylinder filled 
with the liquid and the diameter of the cylinder should be 
at least thirty times that of the sphere. Two liquids of the 
same viscosity should be used so that any differences ob 
served could not be ascribed to viscosity. It is known that 
a pure mineral oil such as Nujol is nonpolar?’ and this 
could be compared with some vegetable oil such as castor 
oil or a solution of stearic acid in caprylic acid which is 
known to be polar. It would be preferable to arrange the 
constants of the system so that the Reynolds’ number is 


small, in order that the theoretical equations could be 





applied to the data. Such an experiment should answer 
unequivocally what Block’s work fails to do. 

The author is indebted to Professor J. C. Hunsaker for 
suggesting this comment and to Professor A. M. Gaudin 
for helpful discussion. 
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New York 

12-21 American Institute of Chemical Engineers, Chi 
cago, Illinois 

19-23 American Society for Metals, Los Angeles, Cali 
fornia 

27-29 American Society for Refrigerating Engineers, Cin 
cinnati, Ohio 

June 

1- 6 Society of Automotive Engineers, White Sulphur 
Springs, West Virginia 

2— 6 American Medical Association, Cleveland, Ohio 

12-20 American Society of Mechanical Engineers—Semi- 
Annual, Kansas City, Missouri 

16-19 American Society of Heating and Ventilating 
Engineers, San Francisco, California 

16-20 American Institute of Electrical Engineers (Sum 
mer Convention), Toronto, Ontario 

18-20 American Physical Society, Pacific Coast Meeting, 
Pasadena, California 
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20-21 American Physical Society, Providence, Rhode 
Island 

23-27 American Association for the Advancement of 
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and more accuracy from science. With 
proper treatment your Eppley precision 
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authoritative sources the emf of cadmium 
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In this advanced text and reference book the author places primary emphasis on dynamic rather than static field theory, 
and the greater part of the book is devoted to the propagation of plane, cylindrical, and spherical waves, the theory of 
radiation, and representative boundary-value problems. The m.k.s. system is used throughout. 
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